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APPENDIX  A 


THE  WILLIAM:;  AFK  AIK  QlJAl.llY  SU'DY 

1.  INTRODUCTION 

This  appendix  begins  with  a  brief  do  script  ion  .<t  r  ij  *  •  Williams  A  F  B  experi¬ 
ment  and  concludes  with  the  presentation  of  stal  is!  icaL  summaries  of  the  data 
acquired  during  the  13-month  monitoring  period  beginning  June  197b.  A  con¬ 
siderably  more  extensive  review  of  th<"  experimental  operations  is  contained 
in  the  paper  of  Sheas ley  et  al.:-4~l  a  report  issued  jointly  bv  Nor  thru; 
Services,  Inc.  and  the  USEPA  Environmental  Monitoring  and  Support  Laboratory . 

2.  MONITORING  LOCATIONS  AND  STATION  SITING  RATIONALEA-2> 

The  placement  of  the  five  ambient  air  quality  monitoring  trailers  :  s 
shown  in  Fig.  A-l .*  The  rationale  for  selecting  these  trailer  locations  i- 
provided  in  Table  A-l.* 

3.  AIR  QUALITY  MONITORING  STATION  DESCRIPTION^-2 

Each  monitoring  station  consists  of  a  mobile  enclosure  (trailer)  approxi¬ 
mately  2.4  m  wide,  4.3  m  long,  and  3.1  m  from  ground  to  roof.  The  following 
monitoring  instruments  are  rack-mounted  in  each  trailer: 

1.  Heckman  Model  6800  Air  Quality  Gas  Ch  roma  tog  rapli  for 
analysis  of  carbon  monoxide,  methane,  and  total  hydro¬ 
carbons.  The  Beckman  6800  has  an  automatic  calibration 
c  a pali  i  1  ity  . 

2.  Monitor  Labs  Model  8440-R/FR/DA  Dual-Chamber  NC/N0X 
Analyzer  with  teflon  filter,  flowmeter,  molybdenum  Nt^-to- 
NO  converter  and  0.05  ppm  full  scale  lowest  range  option. 

3.  R.M.  Young  Co.  Model  35003  Gill  Propeller  Vane  for  wind 
speed  and  direction  measurements. 

4.  MR1  Integrating  Nephelomoter  that  will  be  used  as  a  measure 
of  light  scattering  due  to  particulates. 

Each  trailer  also  has  an  Esterline  Angus  LI  1-028  strip  chart  1  a be r a : o r v 
recorder  for  calibration  and  back-up  purposes. 

The  intake  system  for  the  above  air  quality  instruments  i  n<  ludes  a 
glass  ballast  chamber  (not  in  system  at  present  )  to  provide  a  sample  air 
residence  time  (ratio  of  volume  to  flow  rate)  of  at  least  twice  tie  data 
system  interrogation  rate.  The  chamber  will  allow  the  averaging  o!  tt.iusieut 
high  concentr.it  ions  of  pollutants  that  might  otbeiwise  pass  unobserved.  In 
addition,  the  Heckman  6HU0  has  a  separate  mixing  chamber  to  allow  for  the 

♦Figures  and  tables  appear  consecutively  at  the  end  ot  this  appendix. 


five-minute  instantaneous  sampling  rate.  A  separate  sampling  train  is 
provided  for  the  MRI  integrating  nephelometer  so  that  the  aerosol  size 
distribution  will  not  be  distorted  prior  to  entry  into  the  nephelometer.  A 
"T"  was  installed  in  the  air  quality  sampling  train  so  that,  il  necessary,  air 
could  be  sampled  from  heights  other  than  the  fixed  entry  port  several  teet 
above  the  roof  of  the  trailer. 


4.  PRIMARY  DATA  FLOW 


a.  Aerometric  and  Meteorological  Data 

During  the  period  June  1,  1976,  to  June  3U,  1977,  the  U.S.  KPA/Las  Vegas 
collected  reliable  data  at  Williams  AFB  on  five  pollutants  (CO,  CH/4,  THC ,  NO, 
N0„)  and  b-scatt  plus  wind  speed  and  direction  from  the  five-station  network 
described  briefly  above.  The  raw  data,  collected  at  one-mi nutr  intervals,  was 
corrected  for  span  and  zero  drift  before  computation  of  the  hourly  average 
quan  t it ies  required  for  the  AQAM  evaluation  effort.  The  procedure  for  produc¬ 
ing  the  hourly-average  tape  ( DS  III)  is  indicated  schematically  in  Fig.  A-2 
and  described  in  detail  in  the  report  of  Sheesley  et  al  ,A~1 

Additional  meteorological  data,  taken  routinely  by  the  U.S.  Air  Force  at 
Williams  AFB,  was  recorded  onto  coding  sheets,  keypunched,  and  stored  on 
magnetic  tape.  This  data  set  also  includes  mixing-depth  measurements  taken  by 
the  EPA  with  an  acoustic  sounder.  For  logistical  reasons,  this  data  set  was 
not  merged  (as  erroneously  suggested  in  Fig.  A-2)  with  the  aerometric  d.ita. 
The  format  and  information  structure  of  the  meteorological  data  tape,  along 
with  details  of  the  acoustic  sounder  measurement  program,  are  also  given  by 
Sheesley  et  al  and  further  described  in  Appendix  B. 

b .  Emissions  Data 

Figure  A-2  also  shows  the  essential  steps  taken  in  the  creation  of  the 
source  inventory  data  set  ( DS  I)  and  the  detailed  aircraft  activity  data  set 
(OS  If).  The  source  inventory  data  set  was  created  by  running  the  AQAM  Source 
Emissions  Inventory  Code  on  emissions  data  for  Williams  AFB,  recentLy  compiled 
and  coded  by  SRI  and  checked  for  consistency  by  the  Air  Quality  Research 
Division  of  USAF/CEEDO.  Modifications  to  the  emissions  inventory  are  de¬ 
scribed  in  Appendix  C.  The  aircraft  activity  data  set  (DS  11)  contains 
hour-by-hour  tallies  of  aircraft  operations  (i.e.,  takeoffs,  landings,  touch 
and  goes,  aborts,  and  transient  aircraft  arrivals  and  departures)  broken  down 
according  to  aircraft  type  and  runway  number.  This  data  set  serves  as  input 
to  that  version  of  AQAM  (i.e.,  AQAM  II),  incorporat i ng  the  actual  hourly 
aircraft  emissions  as  opposed  to  the  version  (AQAM  1)  that  uses  the  fractional 
apportionment  of  total  yearly  aircraft  operations  dictated  by  the  AQAM 
source  Emissions  Inventory  Code. 

Generation  of  data  set  DS  II  was  accomplished  by  processing  information 
or.  each  individual  aircraft  operation.  This  information,  extracted  from  IISAF 
Flier. t  Forms  ATC  355,  ATC  90,  and  TW  124  by  USAF  student  pilots  and  CEf-.lKi 
:  ecr  i i  ar.s ,  was  -.od<  -I  and  rteypun'  fieri  according  '•>  the  format  sp.-.  if  fU  it 


I  :il<  1 1-  A  2  .  1 1  n  I  or  t  'in. 1 1  i  1  y  ,  all  l  In-  data  rt-jii  r  <■<!  to  |<r  i.dus  •-  U'.,  !  1  was  not 
i  on  t  a  i  nod  in  available  IISAF  flight  larms  ,  fin-si-  <l<  <  i  t  i  one  i  <-s  o  r  i  g  i  ua  1  I  / 
pointed  out  by  Hum  .in  an<l  bun phy ^  required  tic-  i  mp  I  --.ni-n  t  a  r  i  on  u ;  a  t «-  «/ 
s  tmpl  i  tying  assumptions;  none*  of  which  ant  presently  thought  to  sign  if  it  an:  ly 
alter  the  space-time  distribution  ol  aircraft  ‘“missions  .  These  d*_- 1  i  c  u-:v.  ie  > 
and  approximations  art-  discussed  in  Appendix  D. 


5.  THK  AFROMhTKIC  DATA 


a.  Frequency  Distributions 

As  a  first  step,  preceding  more  detailed  analyses,  on<-  generally  e/ami 
the  frequency  distributions  of  the  measured  quantities.  Vi.,ual  inspection  a: 
such  distributions  quickly  gives  one  an  estimate  ol  the  mean,  median  ,  star, car', 
deviation,  range,  and  statistical  nature  (e.g.,  normal  ,  lognormal.)  trie 

measurables  and  occasionally  reveals  systematic  problems  in  the  data.  rcr 
example:  an  uncharacteristic  bimodal  snape  in  the  station  2,  October  Co  dis¬ 
tribution  led  to  the  discovery  of  an  instrument  calibration  problem  covering  a 
two-week  period.  Though  tl-e  pollutant  distributions  for  each  station  and 
month  were  examined  during  the  diagnostic  phase,  the  distributions  presented 
in  Figs.  A-3  to  A-8  combine  all  stations  and  months,  as  the  subtle  differences 
between  stations  and  months  are  shown  more  clearly  in  the  statistical  tabula¬ 
tions.  Distributions  for  the  "calculated"  pollutants  Nth  (  =  N0X  -  NOJ  and 
NM1IC  (2  THC  -  CH^)  are  given  in  Figs.  A-9  and  A-10,  respectively. 

The  meteorological  variables,  wind  speed  (Fig.  A-1I)  and  wind  direction 
(Fig.  A-I2)  were  also  investigated  in  this  manner.  Of  particular  interest  is 
the  greatly  enhanced  probability  of  winds  from  the  ESF  and  West.  Subsequent 
analysis  of  this  phenomenon  revealed  a  strong  correlation  between  wind  direc¬ 
tion  and  time  of  day.  This  strong  diurnal  dependence  of  the  wind  direction 
(see  Fig.  A-12)  is  characteristic  of  the  mountain-valley  How  found  in  '  in- 
Valley  of  the  Sun,  where  Williams  AFB  is  located. 

Use  lulness  of  the  frequency  distribution  extends  beyond  the  0  i  ;o-:nc '  i  0 
phase  to  the  analysis  phase.  For  example,  a  detailed  tit  of  the  NO  and  N0y 
frequency  distributions,  assuming  only  a  positive  concentration  disir ibut ion 
convoluted  with  a  (latissian  instrumental  resolution  function,  should  disclose 
the  overall  accuracy,  o  ,  of  the  N()/NOx  instrument  near  threshold.  Visual 
examination  ot  the  "nog.it  ive  concentration"  tail  of  thi  s-'  distributions 
suggests  a  o  ot  a  few  ppb . 

Instrument  threshold  and  repeatability  noise  might  also  be  indicated  on 
frequency  distributions  ot  the  data.  If  a  significant  problem,  the  repeat¬ 
ability  noise  resolution  function  could  be  folded  into  tin-  experimental  data 
to  produce  frequency  ideograms. 


b .  St  at  i st  ic a  1  Summa r i os 

Belore  presenting  a  statistical  summary  ol  the  ai-mm.-ti  n  d.il.i,  010  tirct 
must  di-tine  the  range  ot  data  values  that  are  a  c  1  ••  p  I  a  h  1  e  or  pliy.i  .1  i  i 
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allowable  given  the  instrument  range  setting,  instrument  threshold  and  accu¬ 
racy,  as  well  as  "data  windows"  used  in  the  program  that  generated  the  hourly 
average  data  tape  ( DS  III).  While  some  discussion  regarding  the  selection  of 
the  "data  windows"  is  given  by  Sheesley  et  al.,A~l  selection  of  threshold 
concentrations  C-p  (such  that  C^j^  <  C  <  C-p  are  set  equal  to  C-p)  was  based  on  a 
combination  of  conversat  ions  with  the  experimenters,  instrument  manufacturer's 
specifications,  and  a  strong  desire  to  avoid  taking  the  logarithm  (if  a  number 
(necessary  in  the  computation  of  geometric  means)  too  close  to  zero.  The 
minimum,  maximum,  and  threshold  concentrations  for  measured,  as  well  as 
computed,  pollutants  are  given  in  Table  A-3 . 

Statistical  properties  of  the  entire  13-month  data  sample  are  given  in 
Table  A-4 ,  and  for  each  of  the  13  months  individually  in  Tables  A- 3  to  A-17. 
In  these  tables,  the  geometric  mean  and  geometric  standard  deviation  should  he 
considered  meaningless  when  the  arithmetic  mean  is  negative. 

c.  Cumulative  Frequency  Distributions 

Cumulative  frequency  distributions  for  the  13-month  data  sample  are 
presented  (Figs.  A-13  to  A-20)  in  1 ogpr obab i 1 i t y  format  for  each  of  the 
pollutants.  The  interesting  characteristic  of  the  logprobabil ity  plot  is  that 
a  lognormal  distribution  transforms  into  a  straight  line  that  intercepts  the 
titty  percentile  point  at  the  geometric  mean  and  has  a  slope  proportional  to 
the  geometric  standard  deviation. 

These  figures  relate  to  equivalent  tables  (A-18  to  A-25)  at  the  end  of 
this  appendix,  giving  the  concentrations  at  several  useful  cumulative  percen¬ 
tile  levels  for  each  pollutant  species. 
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Table  A-2.  Aircraft  Activity  Data  Format 


Header  Card  -  One  Per  Day 

75122SC  27S  327  129  14  S  9 
Col .  Content 

1-6  YY1MDD  =  Year,  month,  day 

C  =>  Total  aircraft  counts 

8-11  '  30 L 

12-15  '  30C 

I 

16-19  1  30R  Full  traffic  count  on  these  runways 

20-23  }  12R  Rigb1  justify  integers 

24-27  j  12C 

28-31  ;  12L 


Aircraft  Activity  Cards 


75122SN30LF-S 
7S122SI30CC-130 
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0620,0620,1410,1600 


Col,  Content 


1-6 


8-10 
11-20 
21-24  ^ 


YYWIDD 


N  =>  Normal  operations  -  times  are  local  takeoff  and  landing  times 
Slashes  separate  sorties. 

I  =>  Transient  inbound 


0  =>  Transient  outbound 


Times  are  in  GMT  (i.e. ,  Zulu) 


V  =>  Deviation  -  times  are  local  engine  start  and  takeoff  times. 
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A  =>  Non-weather  aborts  -  times  are  local  engine  start  times 


Runway 

Aircraft  type  -  left  justified 


26-50  ' 
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Table  A-13.  ^illU" 


Table  A- 15.  KIUI/- 


Table  A-18.  Cumulative  Frequency  Percentile  Concentrations 
(PPM)  for  Williams  AFB  Hourly  NO  Data: 

June  1976-June  1977 


Percent ile 

Stat ion  No  , 

1 

2 

3 

4 

5 

50. 

0.002 

0.003 

0.002 

0.003 

0.002 

80. 

0.004 

0.004 

0.004 

0.008 

0.005 

90. 

0.006 

0.006 

0.005 

0.012 

0.008 

95. 

0.008 

0.007 

0.007 

0.018 

0.010 

98. 

0.012 

0.011 

0.011 

0.028 

0.013 

99. 

0.016 

0.016 

0.015 

0.038 

0.016 

99.5 

0.023 

0.021 

0.020 

0.048 

0.022 

99.9 

0.052 

0.045 

0.037 

0.075 

0.036 

99  .95 

0.060 

0.050 

0.045 

0.096 

0.045 

99.99 

0.106 

0.115 

0.094 

0.131 

0.059 

Table  A-19. 

Cumulat ive 
(PPM)  for 
June  1976- 

Frequency  Percentile  Concent 
Williams  AFB  Hourly  N0X  Data. 
June  1977 

rat  ions 

Station  No. 

Percentile 

1 

2 

3 

4 

5 

50. 

0.008 

0.009 

0.009 

0.014 

0.008 

80. 

0.014 

0.016 

0.017 

0.027 

0.015 

90. 

0.020 

0.022 

0.024 

0.040 

0.023 

95. 

0.028 

0.030 

0.033 

0.055 

0.032 

98. 

0.039 

0.041 

0.048 

0.080 

0.047 

99. 

0.053 

0.052 

0.061 

0.100 

0.061 

99.5 

0.071 

0.060 

0.076 

0.120 

0.073 

99.9 

0.106 

0.094 

0.106 

0.158 

0.096 

99.9  5 

0.125 

0.107 

0.114 

0.170 

0.102 

99.99 

0.206 

0.228 

0.198 

0.259 

0.108 

I 
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Table  A-20. 

Cumu 1  a  t i ve 
(PPM)  for 
June  1976- 

Frequency  Percent  i 
Will  iams  AFB  Hot.i  1  y 
June  1977 

1  e  Concent 
Nn^  Data: 

ra  t i on s 

Slat  i on  Ne 

Percent  lie 

l 

2 

3 

5 

S 

50. 

0.005 

0.006 

0 . 00  7 

0.010 

u.oot 

O 

00 

0.011 

0.013 

0.01.4 

0.0  2  0 

0.012 

90. 

0.016 

0.019 

0.020 

0  .  cl  ‘  i 

0.020 

95. 

0.022 

0.026 

0.026 

0 . 04  4 

0.028 

98. 

0.034 

0.036 

0 . 04 1 

0.06) 

0 . 04  1 

99. 

0.047 

0.045 

0 . 05  3 

0 .0  74 

0.051 

99.5 

0.058 

0.052 

0.065 

0 . 084 

0.060 

99.9 

0.099 

0.077 

0.088 

0.108 

0.074 

99.95 

0.103 

0.082 

0.099 

0.128 

0.077 

99.99 

0.108 

0.113 

0.108 

0.14  7 

0.097 

Table  A-21 . 

Cumulat ive 
(PPM)  for 
June  1976- 

Frequency  Percent! 
Williams  AFB  Hourly 
June  1977 

le  Concentrations 

CO  Uat a: 

St  at  ion  No 

Pe t  cent  i  1  e 

l 

2 

3 

- 

5 

50. 

0.096 

0.124 

0.081 

0 . 2  m  ; 

0.1  5  i 

80. 

0.168 

0.221 

0 . 1  5m 

0.514 

0.24  4 

90. 

0.262 

0.318 

0.2  52 

0 . 7  58 

0.34  7 

95. 

0.386 

0.444 

0.400 

1  .02  7 

0 . 582 

98. 

0.6  37 

0 .646 

0.600 

1  .  .04 

0 . 

99. 

0.8  56 

0.832 

0.765 

1.77  5 

1.050 

99.5 

0.988 

1  .029 

0 . 9  6  H 

2.182 

1  .261 

99.9 

1  .553 

1  .548 

1  .454 

>  .  140 

1.810 

99.9  5 

1.810 

1  .667 

1  .  b84 

3.56  7 

1  .4  72 

99.99 

2.026 

2.471 

2.261 

1 . J 

2.172 

Table  A-22.  Cumulative  Frequency  Percentile  Concentrations 
(PPM)  for  Williams  AFB  Hourly  CH^  Data: 

June  1976-June  1977 


Percent  lie 

Station  No. 

1 

2 

3 

4 

5 

50. 

1 .615 

1 .603 

1 .  581 

1  .623 

1  .624 

80. 

1  675 

1  .686 

1.654 

1.735 

1  .697 

90. 

1.722 

1.743 

1.712 

1 .861 

1.779 

95. 

1.792 

1 .800 

1.780 

2.010 

1  .882 

98. 

1.891 

1 .893 

1 .890 

2.246 

2.082 

99. 

1.966 

1.966 

1  .966 

2.420 

2.237 

99.5 

2.071 

2.043 

2.069 

2.655 

2.415 

99.9 

2.735 

2.422 

2.392 

3.446 

3.120 

99.95 

3.135 

2.995 

2  .825 

3.689 

3.199 

99.99 

3.339 

3.100 

2.987 

3.980 

3.440 

Table  A-23. 

Cumulative  Frequency  Percentile  Concentrations 
(PPM)  for  Williams  AFB  Hourly  THC  Data: 

June  1976-June  1977 

Percentile 

Station  No. 

1 

2 

3 

4 

5 

50. 

1 .651 

1.704 

1.637 

1.775 

1  .685 

80. 

1.749 

1 .839 

1.760 

1.993 

1  .811 

90. 

l  .829 

1.917 

1 .856 

2.218 

1.938 

95. 

1 .951 

2.028 

1 .974 

2.495 

2.097 

98. 

2.131 

2.161 

2.119 

2.866 

2  .  309 

99. 

2.244 

2.290 

2.253 

3.175 

2  .  522 

99.5 

2.364 

2.369 

2.424 

3.550 

2.693 

99.9 

3.027 

3.265 

3.152 

4.836 

1 . 34  3 

99.95 

3.203 

3.390 

3.406 

5.35  5 

3.516 

99.99 

3.898 

4.328 

4.384 

6.238 

4.006 

99.99 


4  .  006 


Tab  1  e  A-24  .  Cumulative  Frequency  l'e  t  c  eiu  i  1  ■■  Convent,  rat  ions 
(PPM)  for  Williams  A Kb  Hmirlv  NMHC  Data: 

June  197 6- June  1977 


St  a!  i  on  \t.i . 

Percent i le 

1 

2 

3 

U 

50. 

0.048 

0.079 

0.064 

0.1)1 

0.038 

80. 

0.107 

0.176 

0.137 

0.2° ' 

0.116 

90. 

0.155 

0.229 

0 .  i  94 

J  .  <  •  (  1  l) 

11.168 

95. 

0.224 

0.288 

0 . 2  3  7 

0.32/ 

0.2  34 

98. 

0.322 

0.378 

0.331 

0.776 

0. 3  JO 

99. 

0.397 

0.435 

0.417 

)  .  00  7 

0.4  12 

99.5 

0.559 

0.515 

0 . 5  1  6 

1.297 

0. 3  30 

99.9 

0.990 

0.727 

1.222 

2 . 94 1 

0.842 

99.95 

1.152 

0  .842 

1.372 

3.075 

1  .233 

99.99 

1.554 

1.3  30 

2  .00? 

3.98  3 

!  .  7  !4 

Table  A-25. 

Cumulative  Frequency  Cero 
Coefficients  (x  1 0 n-*) 
Hourly  bg^x  Data:  June 

ent  ilt'  Sc  at  t  er  ing 
for  Willi  am  s  A  F  B 

1  9  76-June  1977 

St  at  ion  1 

No  . 

Pe  rcent i  le 

1 

2 

3 

<4 

') 

50. 

0.469 

0.49  5 

0.  51  ) 

0.4  76 

0 . 4  74 

80. 

0.  /  17 

0.809 

0.7  71 

0. 7  81 

0.7  77 

90. 

0 .945 

1  .051 

0.982 

1  .0'.:'  . 

0.948 

95. 

1.207 

1  .352 

1  .  24  3 

1  .  .'Mi 

1  .  24  9 

98. 

1.669 

1  .894 

1  .800 

1.71  4 

1  .  7  84 

99. 

2 .030 

2.333 

2.4  5  3 

2.117 

2.133 

99.5 

2.499 

3.136 

3 . 2  54 

2  ,69b 

2  . 6  06 

99  .9 

3.861 

3 . 9  74 

5.846 

i.  701 

1.466 

99.95 

4.408 

4 . 5  50 

5.962 

4 .19/ 

4.07  1 

99.99 

3.067 

5.828 

7.000 

S  .  6  hf} 

4.468 

APPENDIX  B 


METEOROLOGICAL  DATA:  PROCESSING  AND  ASSUMPTIONS 


1 .  INTRODUCTION 

The  prime  objectives  of  this  research  effort  were  to  evaluate  and 
assess  the  accuracy  of  AQAM,  using  only  the  meteorological  input  data  that 
is  routinely  collected  at  U.S.  Air  Force  bases.  Initial  time  constraints 
dictated  that  the  EPA  collect  and  process  this  meteorological  data  in  a  manner 
similar  to  ETAC;  however,  the  standard  ETAC  hourly  meteorological  data  base 
for  WilLiams  AFB  ultimately  was  used.  These  data,  in  the  DATSAV-SUKFACE 
format,  were  then  merged  together  with  the  EPA  provided  acoustic  sounder  data. 
Additional  EPA/NOAA  acquired  data  on  u,  v,  w  wind  velocities,  solar  radiation, 
and  vertical  temperature  gradient  exists  for  several  months  of  this  study  but 
iias  yet  to  be  integrated  with  the  hourly  average  meteorological  data  base. 

The  method  for  determining  Turner  stability  class  from  the  USAF/ETAC 
data  is  discussed  in  this  appendix.  In  addition,  the  acoustic  sounder  data 
coding  is  interpreted  in  the  context  of  determination  of  an  hourly  mixing 
depth.  The  mixing  depth  is  a  key  AQAM  model  input  provided  neither  directly 
uy  ETAC  or  indirectly  by  a  Short-Term  AQAM  calculation.  Finally,  a  model  lor 
mixing  depths  is  presented,  based  on  ETAC-measur ed ,  surface-level  meteorologi¬ 
cal  parameters. 

2.  METEOROLOGICAL  DATA 

Meteorological  parameters  used  to  determine  stability  class  by  the 
Turner®'*  method  are  wind  speed,  cloud  ceiling  height,  and  amount  of  sky 
cover.  All  parameters  were  supplied  by  USAF/ETAC  but  only  wind  speed  values 
could  be  used  directly.  A  comprehensive  table  of  codes  for  ceiling  heights 
provided  easy  interpretation  of  this  data  field  but  the  sky  cover  data  field, 
also  presented  in  code,  was  interpretable  only  with  the  aid  of  ETAC  personnel 
iGapt.  1 .  Clark,  private  communication,  February  197V). 

Sky  cover  denotes  the  amount  (to  the  nearest  tenth)  of  sky  covered  by 
clouds  or  hidden  by  surface-based  obscuring  phenomena.  In  order  to  completely 
resolve  the  sky  cover  data  field,  three  steps  were  Laken:  tl)  a  description 

of  the  coded  field  was  examined;  (2)  it  was  then  compared  to  the  Turner 
requirement  (see  Table  B-l);  and  (13)  an  estimated  sky  cover  (tenths)  within 
bounds  of  the  requirement  was  assigned.  For  example,  the  ETAC  given  code  tor 
skv  cover  is  -8,  indicating  "broken"  conditions  having  a  sky  cover  of  O.b-U.9, 
inclusive;  but  the  Turner  method  requires  the  range  O.b-1.0,  inclusive,  to 
define  the  same  broken  condition;  therefore,  wherever  a  code  -8  is  supplied, 
the  value  O.b  is  assigned  as  the  sky  rover.  In  the  case  of  "partly  obscured" 
conditions,  we  indicate  only  that  the  amount  of  cover  is  less  than  complete 
though  the  possibilities  span  the  interval  0. 1-0.9.  Table  B-2  gives  the 
frequency  of  observed  sky  cover  codes  for  the  13-month  study  period.  The 
following  table  is  used  to  illustrate  the  assignment  procedure  for  the  six 
different  sky  cover  codes  supplied  among  13  months  ol  hourlv  data. 


ETAC  Code 
(NS) 

Descr i pt i on 

fur  ner 

Kec|u  i  remen  t 

Sky  cove r 

0 

CLEAR.  The  state  oi  tin- 
sky  when  it  is  either 
cloudless  or  the  skv  covet 
is  less  t han  0.1. 

CoV  0.4 

u.  1 

-3 

SCATTERED.  A  skv  cover  of 

0.1  through  0.5  at  and 
between  the  level  ol  a  lavei 
aloft  . 

COV  -  0.5  (day) 

CoV  v  0  .*4  (  n  igti  t  ) 

0 . 4 

-8 

BROKEN.  A  sky  cover  ot  0 . 6 
through  0.9  at  and  below  the 
level  of  a  laver  aloft. 

D.r-  <i’V  1  .u 

(> .  6 

9 

OVERCAST.  A  sky  cover  ot 

1.0  (ten  tenths)  at  and 
below  the  level  of  a  layer 
aloft . 

COV  =  1.0 

l  .0 

8 

OBSCURED.  The  condition  when 
sky  is  completely  hidden  bv 
surface-based  obscuring 
phenomena;  e.g.,  fog,  smoke, 
precipitation  forms,  etc. 

COV  =  1.0 

1  .0 

10 

PARTLY  OBSCURED.  The  con¬ 
dition  when  0.1  or  more  but 
not  all  of  the  sky  is  hidden 
by  surface-based  obscuring 
phenomena . 

COV  <  1.0 

.9 

3.  STABILITY  CLASS  INDEX  AS  DETERMINED  BY  EPA  METHOD 

An  examination  of  the  Single  Source  (CRSTER)  ModelH'2  revealed  that 
one  of  seven  stability  classes  is  determined  from  meteorological  data  for 
each  hour  by  its  preprocessor.  The  first  six  of  these  categories  (i  -4) 
correspond  to  i’asquill's  classifications  (A-F).,}-i  The  seventh  category 
corresponds  to  the  "dashes"  in  I’asquill's  original  classification  and  repre¬ 
sents  the  existence  ol  a  strong,  ground-based  nocturnal  temperature  inversion 
and  nonde  finable  wind  flow  conditions.  Hie  CKS'IEK  model  preprocessor  re¬ 
stricts  changes  in  stability  to  one  class  per  hour  a.-,  show'll  in  Fig.  B-l.  This 
type  of  restriction  allows  undesirable  daytime  stability  classes  such  as  5,  t>, 
and  7.  Figure  8-2  shows  the  result  of  unrestricted  hourlv  class  changes. 

Initially,  the  preprocessor  determines  the  hour  angle  of  the  sun  aid  the 
times  of  sunrise  and  sunset  from  the  day  number,  longitude,  and  time  zone  to 
permit  differentiation  of  daytime  and  nighttime'  cases  hv  l he  Woolf^-i  method. 
For  daytime  cases,  tfie  appropriate  insolation  class  is  selected  by  means  of 
the  Turner  objective  method  using  cloud  cover,  ceiling  height,  and  solar 
elevation  as  indicators.  This  method  assigns  net  radial  mu  indices,  using  the 
criteria  shown  in  Table  B- J ,  for  cases  where  the  total  cloud  rover  '  5  / ID.  It 

the  cloud  cover  >  5/10,  but  less  than  10/10  (overcast),  the  insolation  class 
is  reduced  bv  one  category  when  the  piling  height  is  gi. niter  than  ]<  ,f0f.  it 
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(See  legend  following.) 


i 


ACOUSTIC  MIXING  DEPTH (M) 


LEGEND  FOR  THE  STABILITY  CLASS  INDEX  ANALYSIS 

TOP  GRAPH 

Acoustic  Sounder  determined  mixing  layer  height  (meters) 
Undefined  layers. 

Return  due  to  noise,  possible  turbulence,  rain,  etc. 
Daytime  noise  due  to  airbase  activity. 

Transition  period. 

Mechanical  turbulence  superimposed. 

An  expected  height  value  was  not  supplied. 

Inoperative  equipment. 

Windspeed  (knots) 

BOTTOM  GRAPH 

Stability  class  index  determined  by  EPA  method. 
Stability  class  index  determined  by  ANL  method. 

Cloud  amount  in  tenths  of  total  cloud  cover. 
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The  computed  times  of  sunrise  and  sunset  are 
indicated  by  lowered  numerals. 
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and  by  two  categor  i  es  for  ceilings  between  7,000  and  ift.OOO  tt.  For  c  »_•  1 1  i  nj»  s 
below  7000  ft  and  10/10  cloud  cover  (  i  .  e  .  ,  ovcrcasL ) ,  a  n*t  ladiation  of  0  is 
defined  and  neutral  stability  is  specified.  With  th<-  ox-  option  of  the  i0/10 
low  cloud  cases,  the  net  radiation  index  during  daytime  hours  is  never  rt-dtn.  •  a 
below  1,  or  "weak."  The  final  stability  category  is  selected  from  lable 
and  Turner's  insolation  classes. 


Table  B-l.  Turner  Requirement  Related  to  Skv  Covt ra 


The  net  radiation  index  used  with  wind  speed  to  able ;  r:  stability  class 
is  determined  by  the  following  procedure: 

1)  If  the  total  cloud  cover  is  10/10  and  the  ceiling  is 

less  than  7000  ft,  use  net  radiation  index  equal  to 

0  (whether  day  or  night). 

2)  For  nighttime  (between  sunset  and  sunrise): 

a)  If  total  cloud  cover  4/10,  use  net  radiation 
index  equal  to  -2. 

b)  If  total  cloud  cover  >  4/10,  use  net  radiation 
index  equal  to  -1. 

3 )  For  daytime  : 

a)  Determine  the  insolation  class  number  as  a 
function  of  solar  altitude  from  Table  A-2. 

b)  If  total  cloud  cover  <  5/10,  use  the  net  radiation 
index  in  Table  A- 1  corresponding  to  the  insolation 
class  number. 

c)  If  cloud  cover  >  5/10,  modify  the  insolation  class 
number  by  following  six  steps. 

1)  Ceiling  <  7000  ft,  subtract  2. 

2)  Ceiling  _>  7000  ft  but  <  16,000  ft,  subtract  i. 

3)  Total  cloud  cover  equal  10/10,  subtract  1. 

(This  will  only  apply  to  ceilings  _>  7000  It 
since  cases  with  10/10  coverage  below  7000  ft 
are  considered  in  item  1  above.) 

4)  If  insolation  class  number  has  not  been 
modified  by  steps  (1),  (2),  or  (3)  above,  assume 
modified  class  number  equal  to  insolation  class 
ruinbe r  . 

5)  If  modified  insolation  class  number  is  less  than 
1  ,  let  it  equal  1 . 

6)  Use  the  net  radiation  index  in  Table  A- 1  corre- 
ponding  to  the  modified  insolation  class  number. 


aTaken  from  Ref.  B-l. 
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BE 


j - ^ 


0 

-3 

-8 

8 

9 

10 

1976  June 

545 

138 

28 

9 

July 

54 

594 

87 

9 

August 

191 

512 

29 

8 

1 

3 

September 

122 

494 

78 

22 

4 

October 

300 

341 

85 

17 

1 

November 

464 

195 

35 

26 

December 

391 

232 

70 

51 

1977  January 

127 

267 

80 

28 

2 

February 

2  56 

300 

92 

18 

6 

Ma  rch 

330 

281 

82 

44 

7 

Apr  i  1 

273 

359 

64 

23 

1 

May 

258 

387 

84 

15 

June 

130 

457 

82 

24 

_2 

Total 

3441 

4557 

896 

294 

1 

27 

X  of  Total  ( 92 16 )a 

37.3 

49.4 

9.7 

3.2 

.01 

.29 

aExcludes  11  days  for 

which  pollutant 

observat ions 

are  missing,  19 

January  2,  12,  15,  19 

,  20,  21, 

22,  27 

,  28,  29, 

and 

1977 

June  30. 

Table  B-3  .  Insolation  Classes  as  a  Function  of  Solar 
Altitude  for  Cloud  Cover  <  5 / 1 0 a 


Solar  Elevation 
Angle  (a) 

Insolat  ion 
Class 

Net  Radiation 
Index 

0’  <  a  <  15° 

Weak 

1 

15°  <  a  <  35° 

Si  igh  t 

2 

35°  <  a  <  60° 

Moderate 

3 

60°  <  a 

Strong 

4 

aFor  >  5/10  cloud  cover,  see  Section  3. 


4.  STABILITY  CLASS  INDEX  AS  DETERMINED  BY  ANL  METHOD 

This  section  is  included  mainly  to  document  the  method  used  to  determine 
stability  class  indices  during  the  initial  phase  of  the  Air  Quality  Assessment 
Model  validation  effort. 

In  the  initial  phase  of  the  validation  effort,  pollutant  concentrat ions 
were  estimated  at  nine  receptor  locations.  Subsequent  data  analysis  on  the 
acoustic  sounder  codes  revealed  that  certain  discrepancies  existed  at  and 
around  the  time  of  sunset.  One  obvious  discrepancy  was  the  measurement  of  a 
finite  lid-height  occurring  several  hours  after  the  presumed  sunset  time; 
thus  a  possible  unst  able  condition  was  indicated  as  a  stable  condition. 
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Table  B-4 . 


Stability  Classification  Criteria 


Sur  face 
wind  speed 
(knots) 

Davt ime 

Insol  at  ion 

N i ght t ime 

__  5/10 
Cloud 

<  5/10 

C 1  oud 

Strong 

Moderate 

Si igh  t 

Weak 

Ov  (;  r  c  a  t 

<  1 

1 

1 

2 

3 

4 

h 

7 

“  2 

i 

2 

2 

3 

4 

6 

7 

3 

1 

2 

2 

3 

4 

6 

7 

4 

1 

2 

3 

4 

4 

5 

6 

5 

l 

2 

3 

4 

4 

5 

6 

6 

2 

2 

3 

4 

4 

C 

6 

7 

2 

2 

3 

4 

4 

4 

5 

8 

2 

3 

3 

4 

4 

4 

C 

J 

9 

2 

3 

3 

4 

4 

4 

5 

10 

3 

3 

4 

4 

4 

4 

5 

11 

3 

3 

4 

4 

4 

4 

4 

>  12 

3 

4 

4 

4 

4 

4 

4 

A  comparison  of  the  ANL  method  t  i  tin  i-.i'A  method  showed  that  both  are  the 
same  in  every  essential  except  that  the  ANL  method  assumed  the  fo 1  lowing: 

1.  Designed  for  urban  conditions,  only  the  first  five  sta¬ 
bility  categories  were  allowed. 

2.  The  times  of  sunrise  and  sunset  were  estimated  from  charts 
presented  in  the  Smithsonian  Meteorological  Iai-lis. 

3.  Fractional  wind  speeds  were  rounded  w!m:;  apt.  lie-.!  t' 

Turner's  Stability-Net  Radiation  table. 

Of  the  two  methods,  EPA's  method  was  selected  for  use  in  AOAM  h  •  >  •-  i 

its  more  accurate  determination  of  sunrise  and  sunset  times;  the  r.  ,n  for 
which  can  be  seen  in  Fig.  B-2  where  tne  lid  heigh;  as  by  tie-  a-,  oust  ic 

sounder  is  approximately  130  m  at  the  hour  of  sunset,  1 700  hr,  and  a  ^1  ightly 
unstable  atmospheric  condition  (class  3)  is  correctly  ditermimd  r,\  tin-  EPA 
method.  On  the  other  hand,  a  slightly  stable  atmospher  ic  condit  in,  '■  <•  1  iss  5) 
is  incorrectly  determined  by  the  ANL  method. 

5.  MIXING  DEPTH,  DIRECT  MEASUREMENT  BY  ACOUSTIC  SOUNDER 

The  Acoustic  Sounder  information  is  presented  in  code  as  n  six-digit 

number  of  the  form  ABmDEn,  where  AB  and  UE  are  average  height  values  (fens 
of  meters)  (rum  which  echoes  were  received,  "in"  is  a  code  used  to  explain  a 
special  local  effect,  and  "n"  is  a  code  used  to  designate  a  specific  atmo¬ 
spheric  phenomena.  In  all  cases  of  inoperative  equipment  the  code  is  999994. 

The  average  height  values,  AB  aid  DE,  fall  within  the  20-500  m  vertical 
range  of  the  recording  chart  paper. 


The  special  effects,  m,  cause  graying  of  the  sound  pulse  as  indicated  by 
the  following  four  values: 

1.  Superimposed  return  due  to  noise,  turbulence,  or  prec i pi ta t ion . 

2.  Background  acoustical  noise  due  to  airbase  activity. 

3.  Uncertain  acoustical  return  due  to  unstable-stab le  atmospheric 
trans  it  ion . 

4.  Apparent  mechanical  turbulence  due  to  high  wind  speed. 

The  specific  atmospheric  phenomena,  n,  are  conditions  applied  to  the 
Williams  Air  Force  Base  study.  The  values  of  "n"  (1-10)  are  described  in 

Table  B-5,  which  gives  the  AFB  Acoustic  Sounder  Classification  Codes. 

i'he  table  below  shows  the  assignment  of  mixing  depth  values  based 
on  these  codes.  They  are  used  only  in  the  development  of  a  mixing  depth 
algorithm  (Section  6). 


Acoustic  Sounder  Class 

Mixing  Depth 

(meters ) 

1 

Top  of  radiation  inversion 

0.0 

2 

Drainage  winds 

0.0 

3 

Layered  return 

a.  Top  of  radiation  inversion 

0.0 

b.  Top  of  stratified  echoes 

0.0 

c.  Surface  based 

0.0 

d.  Undefined  layers 

9990. 

e.  Transition  period 

9993. 

f.  Mechanical  turbulence 

ABO.  or 

9994 

i  f 

mi  ssing 

4 

Inversion  base  forced  aloft 

ABO.  or 

9998 

if 

missing 

5 

Subsidence  inversion 

ABO.  or 

9998 

if 

missing 

6 

Frontal  inversion 

ABO.  or 

9998 

if 

missing 

7 

Daytime  noise 

9992  . 

8 

Return  due  to  noise 

9991  . 

9 

Top  of  mixing  layer 

a.  Top  known 

DEO.  or 

9998 

if 

miss i ng 

b.  Top  unknown 

9998. 

c.  Transition  period 

9993. 

d.  Mechanical  turbulence 

DEO.  or 

9998 

i  f 

miss  ing 

10 

Instrument  inoperative 

9999. 

Figures  B-3  and  B-4  present  frequency 

d  i  str ibut ions  of 

acoustic 

sounder 

measurements  based  on  the  hour  of  occurrence  and  the  month  of  occurrence, 
respectively.  These  two  figures  supply  the  bases  for  the  amount  of  present 
data  deemed  valid  using  the  above  interpretation  of  the  Cl  ass i f icat ion  Codes. 


6.  DEVELOPMENT  OF  A  MIXING  DEPTH  ALGORITHM 

For  dayt  imp  hours,  the  acoustic  sounder  provided  mixing  depths  lor  94b 
hours  of  the  13 -month  experiment.  'Ilie  paucity  of  data  is  due  (  •>  t  tie  rapid 
daytime  development  'I  flu-  mixed  layer  arid  t  lie  limited  ranging  ia|i/ihilil  tea 
(  i  .  e  .  ,  -  400  m)  ol  tin  uiuu'ili'  noureler  .  'lit « -  dala  are  lutnjiaied  in  Fig.  It-’) 


Fig.  B-3.  Frequency  Pistribut ion  of  .Voustlc  Sounder  Measurements  bv  Hour 
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ACOUSTIC  SOUNDER  (METERS) 


Table  B-5.  Acoustic  Sounder  Classification  Codes 
for  Williams  AFB  Study 


1.  Top  of  radiation  inversion:  ABO  001 
AB  is  top  in  10s  of  meters 

This  type  of  inversion  is  typical  of  night  time  desert  locations. 

It  is  ground  based  around  sunrise  hours  but  aloft  during  late  evening 
hours;  stable  to  very  stable  atmospheric  conditions.  Much  of  the  turbu¬ 
lent  mixing  occurs  only  in  a  shallow  layer  near  the  ground.  The  thick¬ 
ness  of  this  shallow  layer  varies  over  a  very  wide  range,  from  meters  to 
hundreds  of  meters  making  it  impracticable  to  define  the  layer  height 
from  this  code. 

2.  Drainage  winds:  with  or  without  stratified  layers  above:  ABO  DE2^ 

AB  is  top  of  Drainage  in  10s  of  meters 

DE  is  top  of  stratified  layers  in  10s  of  meters 

a.  AB:  top  of  drainage  ABO  C02 

b.  DE :  top  of  stratified  echoes  ABO  DE2 

c.  top  of  stratified  echoes  unknown  ABO  002 

Drainage  winds  are  defined  for  the  dense  air  that  flows  downhill  in 
mountain-valley  terrain.  They  were  observed,  specifically  during  winter 
months,  to  change  directionally  during  their  westward  onset  at  approxi¬ 
mately  2100  hr,  to  calm  winds,  and  then  eastward  after  2200  hr.(*' 

There  may  or  may  not  be  stratified  layers  above  the  winds  but  the  layer 
below  the  winds  is  assumed  to  be  stable.  The  mixed  layer  as  a  whole  may 
span  several  hundred  meters  in  height  making  it  impracticable  to  define 
the  layer  height  from  this  code. 


Layered  return:  use  when  there  are  layers  in 

addition 

to  radiation 

inversion:  ABO  DE3 

a . 

AB:  top  of  radiation  inversion 

ABO 

003 

b. 

DE:  top  of  stratified  echoes 

ABO 

DE3 

c  . 

surface  based  —  no  radiation  inversion 

000 

DE3 

d. 

undefined  layers 

000 

003 

or 

ABO 

003 

if  radiation 

inversion  is 

determinab le 

e . 

transition  (unstable  to  stable) 

AB  is  top  of  transition  layer 

AB3 

003 

£. 

Mechancial  turbulence  superimposed 

AB  is  top  of  echoes 

AB4 

003 

A  layered  return  is  produced  when  sounding  echoes  are  stratified 
indicating  the  presence  of  multiple  layers.  In  addition  to  the  radiation 
inversion  as  defined  in  1  above,  the  formation  of  multiple  layers  is 
strongly  influenced  by  the  onset  of  drainage  winds. '*)  The  mixed  layer, 
during  these  nighttime  turbulence  episodes,  is  defined  by  category  (f) 
and  will  extend  to  the  top  of  the  echoes. 
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Table  B-b.  (Cont'd) 


Invers ion 

base 

forced  aloft 

by  surface  heating: 

ABU 

D£4_ 

a . 

AB: 

base  of  echo 

ABO 

IM-.4 

b. 

DE: 

top  of 

strat i tied 

echoes 

AilO 

DE4 

c  . 

top  of 

strat ified 

echoes  unknown 

ABO 

004 

This  type  of  inversion  is  due  to  morning  solar  heating  and  is 
limited  to  one  or  two  hours  at  Williams  Air  Force  Base.  The  mixed  layer 
in  this  case  will  extend  to  the  base  of  the  echo. 

5.  Subsidence  inversion:  ABO  DEb 


a  . 

AB: 

base 

of 

echo 

ABO 

U'Jb 

b. 

DE: 

top 

o  f 

stratified 

echoes 

ABO 

DE5 

c  . 

top 

o  f 

strat i f  ied 

echoes  unknown 

ABO 

00  b 

A  subsidence  inversion  is  usually  caused  by  synoptic  stale  high 
pressure  over  the  region.  Since  this  condition  did  not  occur  throughout 
13  months  of  records,  the  assignment  of  mixing  depths  posed  no  problem. 
However,  mixing  depths  would  have  been  treated  in  4  above. 

6.  Frontal  inversion:  ABO  DE6 


a  . 

AB: 

base  of 

echo 

A  BO 

DE6 

b. 

DE: 

top  of 

strat i f ied 

echoes 

ABC) 

DK6 

c  . 

top  of 

strat ified 

echoes  unknown 

ABO 

00b 

Frontal  inversion  (same  as  4  and  5  above  except  for  cause). 

7.  Daytime  noise:  0CE2  000 

Daytime  noise  is  an  undetermined  graying  of  the  return  (background 
acoustical  noise)  due  to  air  base  activity. 

8.  Return  due  to  noise,  possible  turbulence,  rain,  etc.  00_1_  000 

9.  Top  of  mixing  layer  (thermal  spikes),  normal  daytime  situation  ABO  DK2 


a . 

DE :  top  known 

000 

DE9 

b. 

top  unknown 

000 

009 

c  . 

transition  period 

003 

1)E9 

DE  is  top  of  transition  layer  return  (NOTE: 

this 

is  similar  to 

case  3d  above;  however,  thermal  spikes  are  more  prominent.) 

d.  apparent  mechanical  turbulence  due  to  higher 

wind  speed  009  l)E9 

DE  is  top  of  echoes 

During  normal  daytime  situations,  the  mixing  depth  is  unobtainable  only 
if  the  top  is  unknown  (9b)  or  if  the  sounding  return  record  is  uncertain 
( 9c  )  . 

10.  Instrument  inoperative  or  unable  to  read  chart  994  999 

^*^Bowen,  J.,  F.PA  has  Vegas,  private  communication,  March  1979. 


with  the  theory  of  Nozaki,®"^  as  corrected  by  ETAC,  (J.B.  Hooten,  private 
communication,  October,  1978),  which  gives  the  mixing  depth  L  in  meters  as; 

121  ( 6-J )  (DPT)  J(0 . 087 )  (Uz  +  0.5) 

L  =  _  +  _______  . 

where : 

J  =  the  stability  class  index, 

DPT  =  the  dew  point  depression  (°C)  i.e.,  the  difference 

between  surface  temperature  and  dew  point  temperature, 

Uz  *  the  wind  speed  in  knots  at  anemometer  height 
(Z  =  373.4  cm)  , 

ZQ  =  5.0  cm  is  the  roughness  length  characteristic  of  the 
anemometer  exposure,  and 

f  =  0.0001  is  the  Coriolis  parameter. 

The  comparison  shows  that  the  theory  overpredicts  the  observed  mixing  depth 
substantially.  In  addition,  comparison  of  Figs.  B-6  and  B-7  show  the  inabil¬ 
ity  of  this  model  to  duplicate  the  observed  lid  height  time  dependence.  One 
is  cautioned,  however,  that  the  average  height  of  the  lid  during  afternoon 
hours  is  certainly  underpr ed ic ted  by  these  acoustic  sounder  data  as  the 
instrumental  sensitivity  cuts  off  at  500  m. 

A  substantially  better  fit  to  the  acoustic  sounder  data  has  been  obtained 
by  chisquared  optimization  of  the  equation  and  choosing  for  lid  height  the 
maximum  of  the  two  expressions  as  follows: 


t-t 


max 


C1  *  C2  <TC 


T  ) 

MIN 


c5 


C39 


C  41 
4  dt 


or 


H  *  C,  •  WSS 
o 

where : 

C'i  *  66.452  meters 

C2  »  6.9644 

C3  =  1.4402 

C4  *  -14.941 

C5  -  6.1866 

C6  =  8.4085 

t  *  hour  of  day  (0-23) 

T  *  temperature  (K)  at  hour  t, 

WSS  =■  wind  speed  (knots)  at  hour  t, 
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Fig.  B-8.  Mixing  Depth:  A  Comparison  of  Observations  versus 

a  Theoretical  Model 
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1( j  55  the  highest  temperature  attained  up  to  hour  t, 

^MiN  s  the  i°west  temperature  attained  up  to  hour  t, 
tMAX  =  the  h°ur  at  which  Tq  occurred. 

Results  of  this  comparison  are  presented  in  Fig.  B-8 .  Although  based  entirely 
on  surface  observations,  one  obtains  a  correlation  coefficient  of  0.69,  which 
is  comparable  to  results  obtained  using  vertical  temperature  sound ings .  **“6 
This  expression  for  mixing  depth  was  then  used  for  both  AQAM  I  and  II  predic¬ 
tions  throughout  the  13-month  period. 
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APPENDIX  C 


SOURCE  EMISSIONS  INVENTORY  PROGRAM:  INPUT  MODIFICATIONS 

1.  INTRODUCTION 

The  Source  Inventory  Program^-- 1  computes  the  total  annual  pollutant 

emission  for  several  source  categories  at  a  given  airbase,  thus  it  must  be  run 
successfully  before  the  Short-Term  Dispersion  Cocle^_2  can  £>e  usrd . 

The  required  input  to  the  Source  Inventory  Program  was  prepared  at: 
Williams  Airforce  Base,  Ariz.,  by  the  Stanford  Research  Institute  ( SRI ' 

The  input  consists  of  operational  information  related  to  aircraft  and  nonair¬ 
craft  source  activity  (see  Figs.  C-l  to  C-7),*  as  well  as  information  about 
the  airbase  environ  sources  (Fig.  C-8),  such  as  the  Phoenix  area.  Modifica¬ 
tions  to  some  of  the  airbase  operational  information  were  requested  by  the 
U.S.  Air  Force  and  involved  the  addition  of  aircraft  names,  engine  names, 
engine  iden t i f i c a t  i on  numbers,  engine  fuel  rates,  and  engine-pollutant- 
emission  rates.  These  items  were  modified  as  were  the  number  of  airbase 

aircraft  parking  areas  and  their  center  coordinates,  the  number  of  taxiway 
segments  and  their  end  point  coordinates,  and  the  annual  account  of  each 

aircraft  type  and  runway  usage.  In  addition,  data  related  to  aerospace  ground 
equipment  and  aircraft  refueling  operations  were  extended  to  include  three 
additional  aircraft  types.  Finally,  all  training-fire  point  sources  were 
deleted  from  the  input,  though  times  of  actual  training  fires  were  noted. 

The  requested  modifications  are  detailed  in  this  appendix. 


2.  NAMELIST  DATA,  REASSIGNED  PROGRAMMED  DATA 

The  Williams  Air  Force  Base  (WAFB)  AQAM  Source  Inventory  Program  input 
as  supplied  by  Stanford  Research  Institute  (SRI)  has  undergone  several  addi¬ 
tions  and  updates  per  Air  Force  requests. 

The  namelist  input  to  the  Source  Inventory  Program  consists  of  three 
namelist  group  names:  (1)  EGDATA ,  engine  data;  (2)  ACOATA ,  cycle  data; 
and  (3)  DSDATA,  temporal  distribution  data. 

These  data  are  considered  to  be  good  overall  averages  of  aircraft 
engine  emission  factors,  aircraft  landing  and  takeoff  parameters,  and  the 
temporal  distribution  of  aircraft  and  airbase  activities;  thus  they  are 
automatically  used  by  the  program  when  the  user  has  not  input  other  data. 

EGDATA  is  the  only  group  name  considered  in  the  update  process  for  which 
the  reassignments  involve  the  addition  of  aircraft  name's,  engine  names,  engine 
ID  numbers,  fuel  rates,  and  engine  pollutant  emission  data. 

Two  table's,  Tables  C- 1  and  C-2,*  are  presented  to  record  the'  change 
from  the  SRI  supplied  input  to  the  updated  input.  Note-  that  the  e*xtende‘d  list 


♦Figures  and  tables  appear  consecutively  at  the  end  of  this  appendix. 


of  aircraft  types  results  from  actual  observations  at  WAFB  and  from  guidelines 
submitted  by  the  Air  Force  dated  12  October  1978. 

3.  DATA  SET  4,  AIRBASE  AIRCRAFT  AND  RUNWAY  TOTALS 

The  information  coded  in  this  data  set  defines  the  total  number  ol 
aircraft  types,  runways,  parking  areas,  a  special-case  wind  condition,  and 
taxiway  segments.  The  following  table  shows  the  value  of  these  parameters  for 
tiie  SRI  and  updated  inventories. 


SRI  Update 

Aircraft  Types  4  6 
Runways  6  6 
Parking  Areas  4  3 
Special-Case  Wind  1  1 
Taxiway  Segments  16  24 


4.  DATA  SET  5,  AIRCRAFT  ACTIVITY 

Activity  for  aircraft  is  defined  as  the  total  annual  number  of  arrivals, 
departures,  and  touch-and-go  operations.  Each  arrival  and  each  departure  of 
an  aircraft  is  considered  as  an  operation  in  the  landing  and  takeoff  cycle.  A 
touch-and-go,  however,  is  a  complete  cycle  and  occurs  during  a  training  flight 
when  an  aircraft  approaches  and  lands  on  a  runway,  travels  down  the  runway  for 
several  seconds,  accelerates,  and  lifts  off. 

Estimates  of  activity  for  aircraft  were  made  from  actual  counts  at  WAFB 
taken  during  the  period  1  June  1976  to  30  7'nc  1977.  A  review  of  counts  of 
June  and  July  1976  shows  that  the  information  needed  to  determine  touch-and- 
go  operations  is  presented  as  "daily  estimated  averages  of  total  runway 
activity."  These  daily  estimated  averages  could  not  be  used  realistically  in 
the  summation  methods  required  for  an  annual  number  of  operations.  Therefore, 
the  activity  for  aircraft  is  redefined  to  reflect  twice  the  activity  informa¬ 
tion  presented  for  the  6  month  period  August  1976-January  1977. 

Tables  C-3  and  C-4  are  presented  to  record  the  change  from  the  SRI 
supplied  activity  to  the  updated  activity,  respectively.  Note  the  new 
TRANSIENT  aircraft  category  consists  of  the  F4,  C130H,  and  C140  aircraft 
types.  In  regard  to  touch-and-go  cycles,  an  account  is  made  for  both  normal 
and  transient  operations  by  F5,  T37,  and  T38  aircraft. 

5.  DATA  SET  6,  AIRCRAFT  PARKING  AREAS 

Information  in  this  data  set  must  describe  the  geometries  of  the  aircraft 
parking  areas.  The  center  coordinate  of  each  square  in  a  series  of  squares 
making  up  Che  parking  area  is  defined  along  with  the  length  of  a  side  of  each 
square.  The  model  assumes  the  square  is  situated  so  that  a  line  drawn  paral¬ 
lel  with  its  right  or  left  side  wil1  N  directed  north-south. 


A  6 


Tables  C-5  and  C-6  record  the  change  from  the  SKI  supplied  geometries  to 
the  updated  geometries.  Figure  C-9  shows  the  location  of  aircraft  parking 
areas  relative  to  taxiways,  runways,  and  ambient  air  monitors. 

6.  DATA  SET  7,  AIRCRAFT  TAXIWAY  PATH  SEGMENTS 

Each  taxiway  path  used  by  aircraft  at  WAFB  is  defined  as  a  series  of 
connected  straight  line  segments.  Information  in  this  data  set  describes  the 
geometries  of  these  straight  line  segments  and  assigns  to  each  an  identifying 
line  number.  The  line  number  is  used  for  defining  the  particular  segments 
that  will  be  used  to  make  up  a  complete  taxiwav  path. 

Tables  C-7  and  C-8  record  the  change  from  the-  f'K]  supplied  geome  t  r  i  <- and 
the  updated  geometries.  Figure  C-10  shows  the  nunibe; e  i  taxiway  segments 
relative  to  the  airbase  runway  configuration. 

7.  DATA  SET  8,  AIRCRAFT  RUNWAY  INFORMATION 

All  information  concerning  airbase  runways  is  defined  in  this  data  set. 
Annual  arrivals  and  departures  are  defined  for  each  runway  (see  Tables  C-9  and 
C-1C)  as  well  as  length  and  direction.  Its  length  is  the  physical  length  of 
the  runway  pavement;  its  direction  is  determined  by  its  orientation  in  rela¬ 
tion  to  true  north  (see  Tables  C— 1 1  and  C-12). 

Each  individual  aircraft  type  may  use  several  different  runways  for 
landing  and  takeoff  operations.  Therefore,  inbound  and  outbound  taxiway  paths 
are  defined  to/ from  each  aircraft  parking  area  to/ from  every  airbase  runway  by 
a  sequence  of  taxiway  segment  numbers. 

Table  C— 1 3  records  the  update  of  aircraft  movement  over  the  taxiwav 
segment  configuration  (see  Fig.  C-10). 

8.  AEROSPACE  GROUND  EQUIPMENT  EMISSIONS 

Aerospace  Cround  Equipment  (AGE)  consists  of  all  motorized  equipment 
except  refueling  tanks  used  to  support  incoming  and  outgoing  aircraft.  These 
support  vehicles  generally  consist  of  coolers,  power  generators,  heaters,  and 
hydrant i-  test  stands.  The  emissions  for  this  equipment  must  be  determined 
and  inpu'  directly  into  the  model.  The  model  assumes  that  all  AGE  activities 
nii  u  in  the  aircraft  parking  areas,  but  the  emissions  are  calculated  sepa- 
r •  .v  trim  those  emissions  resulting  from  aircraft  parking  activities.  Since 
•  •  '  it.  -'t  this  data  set  involved  only  the  replacement  of  the  single 

-*"■  V  a'  •  !■  r  v  with  K4 ,  C13QH,  and  C141  aircraft  types,  Table  '  14  is 

. t  i  'he  up*! a  t  ed  i  npn t  . 

•  ‘V.  SP  1 1  ..Aid. ,  e.Nlt  VhNIINi.  lOIAI.S 

•  •  !  'i'  i  :  e  ■  aid  vent  ng  i  r  i  i  i  !  t  i 
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fuel  drained  from  the  aircraft  fuel  lines.  Since  the  update  of  this  data  set 
involved  only  the  replacement  of  the  single  TRANSIENT  category  with  FA,  C130H, 
and  C1A1  aircraft  types,  Table  C— 1 5  is  presented  to  record  the  updated  input. 

10.  DATA  SET  13,  TRAINING  FIRE  POINT  SOURCES 

Training  fire  point  sources  are  defined  as  shallow  ground  level  sites  on 
the  airbase  that  are  filled  with  fuel  and  ignited  for  the  purpose  of  training 
airbase  personnel  in  the  art  of  fire  fighting. 

The  update  to  this  data  set  involved  the  deletion  of  all  sites  (a  total 
of  4),  designated  by  SRI,  from  the  inventory.  Thus,  training  fire  point 
sources  are  not  considered  in  the  WAFB  modeling  effort. 
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Fig.  C-l.  Point  Source  Locations  Relative  to  Tax i wavs 
Runways ,  and  Anb i en t  Air  ion itors 
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iDient  Air  Monitors 


Table  C-l.  Engine  Data  (SKI) 


(ACNAME) 

AIRCRAFT 

NAME 

AIRCRAFT 

ID 

( ECNAME) 
ENGINE 
NAME 

( IDACEF) 
ENGINE 

ID 

THRUST 

SETTING 

(EGFF) 

FUEL  RATE 
1000  LB  /UR 

Pollutant  Emission 
Pounds  per  1000  los 
(EGEMEC) 

Data  In 
of  Fuel 

CO 

HC 

NOX 

PM 

Transient 

50 

J79-C15 

i 

Idle 

1.131 

56.70 

1  !>.  70 

1.47 

0.50 

Normal 

2.720 

11  .-0 

1.33 

4.25 

2.12 

Military 

8.921 

2.28 

0.22 

8.94 

2.  36 

After  Br 

32.240 

4 . 00 

0.01 

3.11 

0.15 

F5 

13 

J85 

6 

Idle 

0.433 

180.00 

79.90 

l  .26 

0.013 

T38 

32 

Normal 

1.462 

4  ) .  30 

1.  3  7 

2.32 

0.017 

Military 

2.630 

29.  30 

•t.84 

2.b  8 

0.015 

After  Br 

8.323 

26.00 

0.07 

1.99 

0.008 

T37 

31 

.169 

14 

Idle 

0.231 

127.00 

19.50 

1.53 

0.  729 

Normal 

0.698 

49.10 

1.29 

2.67 

0.017 

Military 

1.093 

31 . 30 

0.30 

) .  oO 

0.020 

Table 

C-2. 

Engine 

Data  (Update) 

(ACNAME  3 
\IRCRAET 
NAME 

AIRCRAFT 

ID 

( ECNAME) 
ENGINE 
NAME 

(IDACEF) 

ENGINE 

ID 

THRUST 

SETTING 

(EGFF) 

Follutant  Emission 
Pounds  per  1000  lbs 
( EG  DIF  C) 

Data  ir. 
of  Fuel 

1000  LB  /HR 

CO 

HC 

NOX 

TM 

F4 

12 

J79-G15 

1 

Idle 

1.130 

57.00 

12.00 

2.50 

0.50 

Normal 

3.300 

9.40 

1.10 

4 . 8U 

1  .  do 

Military 

8.929 

2.20 

0.20 

8.90 

2  .  20 

After  Br 

32.240 

4.00 

0.01 

i.  1  1 

0 .  I  5 

(1141 

14 

TF33-P3 

4 

Idle 

0.900 

84  .00 

10  7.  0 

1.8 

U .  2  4 

Normal 

3.797 

6.  30 

2.6 

3.8 

0.99 

Military 

7.4  36 

1  .  ?() 

O.b 

10.0 

1.73 

E5 

1  3 

J85 

6 

Idle 

0.469 

1 78.  U 

JO.O 

1.30 

0.00  4 

T  38 

32 

Normal 

1.000 

7  3.6 

6.40 

i  .80 

0.00  7 

Military 

3.527 

29.(3 

0.80 

7 . 60 

0.018 

After  Br 

8.32  3 

26.0 

.07 

1.99 

0 . 008 

T37 

31 

J69 

14 

Idle 

0.231 

129.0 

19.0 

1  .  5 

0 .  3  r> 

Normal 

0.288 

107.0 

11.1 

1  .  7 

0.28 

Military 

1.095 

32.0 

0.5 

; .  b 

0.0  7 

C130H 

44 

T56-A15 

21 

Idle 

0.493 

18.1 

15.1 

2.4 

U.  48 

No rmal 

0.827 

8.5 

3.4 

7 

0.47 

Military 

2.392 

1 . 6 

0.2 

11.7 

0.7] 
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Table 

C-3.  Aircraft  Activity  (SRI) 

Aircraft 

Name 

Annual  Number  of 

Arrivals 

Departures 

T/G  Cycles 

1. 

F5 

8391 

8391 

1682 

2. 

T37 

39662 

39662 

63357 

3. 

T38 

77281 

77281 

55963 

4. 

TRAN 

8112 

8112 

0 

Table  C-4.  Aircraft  Activity  (Update) 


Aircraft 

Name 

Annual  Number  of 

Arrivals 

Departures 

T/G  Cycles 

1. 

F5 

6566 

6666 

15493 

2. 

T37 

25394 

25904 

34614 

3. 

T38 

27180 

27686 

64125 

4. 

F4 

1310 

1320 

0 

5. 

C130H 

66 

62 

0 

6. 

C140 

20 

22 

0 

Table  C-5. 

Aircraft 

Parking  Area  Geometries  (SRI) 

Center 

Coordinate  (UTM) 

Area 

Square 

X 

Y 

Length  (km) 

1 

1 

438.10 

3685.57 

0.10 

2 

1 

437.70 

3686.10 

0.20 

2 

437.90 

3685.90 

0.20 

3 

438.00 

3685. 70 

0.20 

3 

1 

437.95 

3684.70 

0.20 

2 

438.10 

1684.40 

0.20 

3 

438.11 

1685.05 

0 . 1 4 

4 

1 

4  $8.  i  5 

Ifi8  3  .  3 

0.  If 

> 

438.1  ' 

i  h8 .  r> 

0.1, 

8ft 
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Table  C-6 .  Aircraft  Parking  Area  Geometries  (Update) 


Area 

Square 

Center  Coordinat 

X 

e  (UTM) 

Y 

Length  (km) 

1 

1 

437. 41 

36  86. 10 

0.13 

2 

437.54 

3685.97 

0.13 

3 

437.67 

3685.84 

0.13 

2 

1 

437.83 

3685.50 

0.  ]  3 

2 

437.83 

3685. 30 

0.  13 

3 

437.83 

3685.10 

0.13 

3 

1 

437.71 

3684.81 

0.13 

2 

437.58 

3684.68 

0.13 

Table 

C-7.  Taxiway 

Segment  Geometries  (SRI) 

Line 

No. 

Ground  Level 
of  One  End 

Coordinates 
of  Line 

Ground  Level 
at  Opposite 

Coordinates 
End  of  Line 

Segment 
Length  (km) 

X  (1) 

Y  (1) 

X  (2) 

Y  (2) 

1 

440.800 

3684.430 

440.560 

3684.180 

0.347 

2 

440.560 

3684.180 

440.050 

3683.750 

0.667 

3 

439.760 

3683.980 

440.050 

3683.750 

0.370 

4 

439.760 

3683.980 

439.940 

3684.140 

0.241 

5 

439.760 

3683.980 

438.530 

3685.230 

1.754 

6 

438.530 

3685.230 

438.230 

3685.580 

0.461 

7 

438.530 

3685.230 

438.130 

3685.250 

0.400 

8 

438.130 

3685.250 

438.080 

3684.820 

0.433 

9 

438.130 

3685.250 

438.130 

3685.500 

0.250 

10 

438.130 

3685.500 

438.230 

3685.580 

0.128 

11 

438.230 

3685.580 

437.870 

3685.900 

0.482 

12 

438.230 

3685.580 

437.520 

3686.380 

1.070 

13 

437.870 

3685.900 

437.520 

3686.380 

0.594 

14 

437.520 

3686.380 

437.730 

3686.380 

0.210 

15 

437.730 

3686.380 

438.400 

3686.380 

0.670 

16 

438.400 

3686.380 

438.830 

3686.380 

0.430 

81 


Table  C-8.  Taxiway  Segment  Geometries  (Update) 


Line 

No. 

Ground  Level  Coordinates 
of  One  End  of  Line 

Ground  Level  Coordinates 
at  Opposite  End  of  Line 

Segment 
Length  (tcm) 

X  (1) 

Y  (1) 

X  (2) 

Y  (2) 

1 

437.570 

3684.600 

437.900 

3684.920 

0.460 

2 

437.900 

3684.920 

438.040 

3684.790 

0.191 

3 

437.900 

3684.920 

437.910 

3685.330 

0.410 

4 

437.910 

3685.330 

438.050 

3685.330 

0.140 

5 

437.910 

3685.330 

437.920 

3685.720 

0.390 

6 

437.920 

3685.720 

437.600 

3686.050 

0.460 

7 

437.600 

3686.050 

437.220 

3686.460 

0.559 

8 

437.220 

3686.460 

437.460 

3686.460 

0.240 

9 

437.460 

3686.460 

438.100 

3686.460 

0.640 

10 

438.100 

3686.460 

438.530 

3686.460 

0.430 

11 

439.100 

3685.870 

438.900 

3685.650 

0.297 

12 

438.900 

3685.650 

438.530 

3685.330 

0.489 

13 

438.530 

3685.330 

438.270 

3685.330 

0.260 

14 

438.530 

3685.330 

438.410 

3685.150 

0.216 

15 

438.410 

3685.150 

438.040 

3684.790 

0.516 

16 

438.410 

3685.150 

438.270 

3685.330 

0.228 

17 

438.270 

3685.330 

438.050 

3685.330 

0.220 

18 

438.050 

3685.600 

438.050 

3685.330 

0.270 

19 

437.920 

3685.720 

438.050 

3685.600 

0.177 

20 

438.410 

3685.150 

439.470 

3684.090 

1.499 

21 

439.470 

3684.090 

439.730 

3683.830 

0.354 

22 

439.730 

3683.850 

440.220 

3684.300 

0.665 

23 

440.220 

3684.300 

440.470 

3684.500 

0.320 

24 

439.470 

3684.090 

439.580 

3684.230 

0.178 

Table  C-9.  Annual  Aircraft  Arrivals  and  Departures  (SRI) 


Runway 

F5 

T37 

T38 

TRAN 

30L 

Arrivals 

0 

34515 

0 

0 

Departures 

0 

3451  5 

0 

0 

30C 

Arrivals 

2960 

4353 

24917 

7950 

Departures 

2960 

4353 

24917 

7950 

3  OR 

Arrivals 

5263 

0 

50818 

0 

Departures 

5263 

0 

50818 

0 

12K 

Arrivals 

0 

704 

0 

0 

Departures 

0 

704 

0 

0 

12C 

Arrivals 

60 

704 

508 

324 

Departures 

60 

704 

508 

324 

12L 

Arr  ivals 

107 

0 

1037 

0 

Departures 

107 

0 

1037 

0 

Table  C-10. 

Annual 

Aircraft 

Arrivals 

and  Departures 

(Updat 

Runway 

F5 

T37 

T38 

F4 

Cl  3  OH 

C141 

30L 

Arrivals 

0 

24608 

0 

0 

66 

18 

Depnr tures 

0 

25140 

0 

0 

60 

18 

30C 

Arrivals 

5084 

0 

21962 

1270 

0 

0 

Departures 

6356 

0 

26884 

1286 

0 

0 

30R 

Arr  ivals 

1180 

0 

4396 

0 

0 

0 

Departures 

0 

0 

0 

0 

0 

0 

12R 

Arrivals 

0 

786 

0 

0 

0 

2 

Departures 

0 

764 

0 

0 

2 

4 

12C 

Arr ivals 

260 

0 

654 

40 

0 

0 

Departures 

310 

0 

802 

34 

0 

0 

12L 

Arr i vals 

42 

0 

168 

0 

0 

0 

Departures 

0 

0 

0 

0 

0 

0 

— 

-  .  —  - - - 

_ _ _ _ __  . 

-  - - 

_ 

_ 

_ 
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Table  C-ll.  WAFB  Runway  Geometries  (SRI) 


Runway 

Coordinate  (UTM) 

X  Y 

Angle 

(Deg) 

Length 

(km) 

30L 

439.94 

3684.14 

315 

3.17 

30C 

440.56 

3684.18 

315 

3.12 

30R 

440.80 

3684.43 

315 

2.85 

12R 

437.73 

3686.38 

135 

3.17 

12C 

438.40 

3686.38 

135 

3.12 

12L 

438.83 

3686.38 

135 

2.85 

Table  012.  WAFB  Runway  Geometries  (Update) 


Coordinate  (UTM) 
Runway  X  Y 


Angle 

(Deg) 


Length 

(km) 


30L 
30C 
3  OR 
L2R 
L2C 
1 2L 


439.58 

3684.23 

316 

3.08 

440.22 

3684.30 

316 

3.03 

440.47 

3684.50 

316 

2.76 

437.46 

3686.46 

136 

3.08 

438.10 

3686.46 

136 

433.53 

3686.46 

136 

2.76 

Tabic  C-13. 

Aircraft  Movement 

over  Taxiway  Segments 

Aircraft 

Annual 

Inbound 

Parking 

Annual 

Outbound 

Runway 

Name 

Arrivals 

Taxi  Segments 

Area 

Departures 

Taxi  Segments 

30L 

T37 

24608 

B.7 

1 

25140 

6,19,18,17,16 

20,24 

C130H 

66 

8, 7,6,5 

2 

60 

4,17,16,20,24 

C140 

18 

8,7,6, 5 

2 

18 

4,17,16,20,24 

30C 

F5 

5084 

9, 8, 7. 6, 5 

2 

6356 

4,17,16,20,21 

22 

T38 

10981 

9,8, 7, 6, 5, 3 

2 

13442 

3,4,17,16,20, 

21,22 

T38 

10981 

9, 8, 7, 6, 5, 3,1 

3 

13442 

1,3,4,17,16, 

20,21,22 

F4 

1270 

9, 8, 7, 6, 5 

2 

1286 

4,16,16,20,21 

22 

30R 

F5 

1180 

10,9,8,7,6,5 

2 

0 

T38 

2198 

10,9,8,7,6,5,3 

2 

0 

T38 

2198 

10,9.8,7,6,5,3, 

1 

3 

0 

12R 

T37 

786 

24,20,16,17,18, 

19,6 

1 

764 

7,8 

C130H 

0 

2 

2 

5, 6, 7, 8 

C141 

2 

24,20,16,17,4 

2 

4 

5, 6, 7, 8 

12C 

F5 

260 

22,21,20,16,17, 

4 

2 

310 

5, 6, 7, 8, 9 

T38 

318 

22,21,20,16,17 

4,3 

2 

401 

3, 5, 6, 7, 8, 9 

T38 

318 

22,21,20,16,17, 

4,3,1 

3 

401 

1 ,3, 5, 6, 7, 8, 9 

F4 

40 

22,21,20,16,17, 

4 

2 

34 

5, 6, 7, 8, 9 

12L 

F5 

42 

23,22,21,20,16, 

17,4 

2 

0 

T38 

84 

23,22,21,20,16, 

17,4,3 

2 

0 

T38 

84 

23,22,21,20,16, 

17,4,3,1 

3 

0 
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Table  C-14.  Service  Vehicle  Emissions  in 
Kilograms  per  Operation 
(Arrival  or  Departure) 


CO 

THC 

NOX 

PM 

SOX 

Gasoline : 


F5 

0.674 

0.040 

0.002 

0.005 

0.001 

T37 

0.133 

0.008 

0.001 

0.001 

0.001 

T38 

0.674 

0.040 

0.002 

0.005 

0.001 

F4 

0.674 

0.040 

0.002 

0.005 

0.001 

C130H 

0.674 

0.040 

0.002 

0.005 

0.001 

C140 

0.674 

0.040 

0.002 

0.005 

0.001 

JP4 : 

F5 

0.072 

0.007 

0.014 

0.007 

0.002 

T37 

0.003 

0.001 

0.001 

0.0001 

0.001 

T38 

0.072 

0.007 

0.014 

0.007 

0.002 

F4 

0.072 

0.007 

0.014 

0.007 

0.002 

C130H 

0.072 

0.007 

0.014 

0.007 

0.002 

C141 

0.072 

0.007 

0.014 

0.007 

0.002 

Table  C-15.  Refueling  Information 
(JP4-Jet  Fuel) 


F5 

T37 

T38 

F4 

C130H 

C141 

Refueling  Value  2146 

(liters  per  fillup) 

587 

1105 

1893 

1893 

1893 

Fuel  Spillage  4 

(Liters  per  fillup) 

10 

4 

4 

4 

4 

Fuel  Venting  2 

(Liters  per  arrival) 

0 

2 

2 

2 

2 

Fuel  Venting  0 

(Liters  per  departure) 

2 

0 

2 

2 

2 

APPENDIX  D 


AIRCRAFT  OPERATIONS  DATA 


1.  OBJECTIVE 

The  objectives  of  this  appendix  are:  first,  to  develop  plans  and 

procedures  for  the  acquisition  and  reduction  of  aircraft  time-in-mode  and 
aircraft  activity  data;  second,  to  present  a  summary  of  an  aircraft  time-in- 
mode  study  performed  at  Williams  AFB;  and  third,  to  present  a  summary  cf 
aircraft  activity  at  Williams  for  the  13-month  period  of  the  EPA  air  quality 
monitoring  experiment.0-^ 


2.  BASE  OPERATIONS 

Williams  Air  Force  Base  is  a  relatively  remote,  air-training  command  base 
near  Phoenix,  Arizona.  The  main  mission  of  the  base  is  the  training  of  jet 
fighter  pilots.  There  are  three  squadrons  at  this  airbase,  each  with  its  own 
type  of  aircraft  (Fig.  D-l),  and  these  aircraft  account  for  most  of  the  air 
traffic:  approximately  2100  T-37  flights,  2300  T-38  flights,  and  550  F-5 

flights  per  month,  depending  on  the  number  of  students  in  training  at  a 
particular  time. 

The  layout  of  runways  and  parking  areas  for  squadron  aircraft  is  illus¬ 
trated  in  Fig.  D-2.  Only  three  runways  are  used  at  any  giver,  time.  Runways 
prefixed  by  a  12  are  used  only  when  the  northwesterly  wind  exceeds  10  knots. 
The  suffixes  L,  C,  and  R  denote  left,  center,  and  right,  respectively.  Runway 
30L  is  used  almost  exclusively  by  T-37  aircraft.  Runways  30C  and  30R  are 
used  jointly  by  T-38  and  F-5  aircraft.  Transient  aircraft  operations  are 
normally  performed  on  runway  30C.  A  schedule  of  aircraft  operations  by 
aircraft  type  and  runway  is  provided  in  Table  D-l. 


3.  DATA  ACQUISITION  AND  REDUCTION  OF  AIRCRAFT  TIME-IN-MODE 

Time-in-mode  is  the  time  spent  in  modes  of  ground  operations.  A  mode  is 
simply  a  subdivision  of  an  aircraft  ground  (or  air)  operation.  These  times 
can  be  characterized  statistically  as  an  average  for  each  group  of  time-in- 
mode  observations.  Naugle0-^  has  collected  aircraft  time-in-mode  data  and  has 
given  estimates  of  time-in-mode  for  several  aircraft  types  at  five  Air  Force 
bases.  This  present  effort  was  made  because  of  the  need  for  accuracy  in 
source  inventory  data  to  be  used  for  Williams. 

The  results  of  a  six-day  study  of  aircraft  time-in-mode  at  Williams  Air 
Force  Base  are  presented  in  Table  D-2.  Only  time-in-mode  observations  for 
squadron  aircraft  ground  operations  were  made.  Observations  of  idle  at 
start-up  times  were  made  from  the  parking  areas;  all  other  time-in-mode 
observations  were  made  from  the  control  tower.  Approximately  1200  separate 
aircraft  mode  observations  were  made  during  the  study. 

The  time-in-mode  observations  were  stratified  according  to  aircraft,  run¬ 
way,  taxiway  segment,  pad  area,  and  aircraft  mode.  A  diagnostic  investigation 
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F-5  Aircraft 


Fig.  D- 


.  Aircraft  Squadrons  at  Williams  Air  Force  Base 


b8 


Table  D-l .  Schedule  of  Aircraft  Operations 


Schedule 

Runway 

Ai  rcraf t 

Operat ion 

L 

30L 

T-3  7 

Takeoff 

2 

301 

T-3  7 

Land ing 

3 

30  L 

T-3  7 

Touch-and-Go 

4 

12R 

T-37 

Takeof  f 

5 

12R 

T-37 

Land ing 

6 

12R 

T-37 

Touch-and-Go 

7 

30R 

T-38 

Takeof  f 

8 

30R 

T-38 

Land ing 

9 

30R 

T-38 

Touch-and-Go 

10 

30C 

T-38 

Takeoff 

11 

30C 

T-38 

Landing 

12 

30C 

T-38 

Touch-and-Go 

13 

12L 

T-38 

Takeof  f 

14 

12L 

T-38 

Landing 

15 

12L 

T-38 

Touch-and-Go 

16 

12C 

T-38 

Takeof f 

17 

12C 

T-38 

Landing 

18 

12C 

T-38 

Touch-and-Go 

19 

30R 

F-5 

Takeoff 

20 

30R 

F-5 

Landing 

21 

30R 

F-5 

Touch-and-Go 

22 

30C 

F-5 

Takeoff 

23 

30C 

F-5 

Landing 

24 

30C 

F-5 

Touch-and-Go 

25 

12L 

F-5 

Takeof  f 

26 

12L 

F-5 

Landing 

27 

12L 

F-5 

Touch-and-Go 

28 

12C 

F-5 

Takeoff 

29 

12C 

F-5 

Landing 

30 

12C 

F-5 

Touch-and-Go 

31 

30C 

Transient 

Takeof  f 

32 

30C 

Trans ient 

Landing 

33 

12C 

Transient 

Takeoff 

34 

12C 

Trans ient 

Landing 

of  the  observations  for  each  aircraft  mode  was  performed  to  provide  insight 
into  the  distributional  characteristics.  The  median,  mean,  standard  devia¬ 
tion,  and  the  number  of  observations  were  computed.  Observations  of  t  ime- 
in-tnode  were  contaminated  by  queueing  situations,  aborted  operations,  and 
emergency  situations. 

An  examination  of  the  size  of  the  standard  deviation  relative  to  the 
mean  reveals  that  the  inbound  modes  immediately  preceding  takeoff  and  the 
outbound  modes  immediately  following  landings,  with  the  exception  of  taxiing 


Table  D-2.  Time-in-Mode  Data 


Sumbe  r 

Throttle  of  Standard 

operation  -Aircraft  Mode _ Setting  Observations  Median  Mean  Pevininn  Kxplanati-.n 

1-  17 Aircraft 


Outbound 

Idle  at  Startup 

ah 

24 

26  1.00 

266 . be 

49 .  c  V 

Vaxi 

S5 

n 

9Q i . on 

3  /  1.61 

.w  .96 

'•one 

rn!  laxx  to  hast  Chance 

Taxi 

55 

16 

152.00 

160  289 

37.42 

Cene 

ra!  Taxi  to  Segment  5 

l’axt 

55 

14 

I 89 . 00 

190.9 1 

36 . 4  3 

Taxi 

i .  •  S’.-int  b  froir.  Pad  1/ 

Taxi 

bb 

17 

no.  oo 

1 34 . 40 

In.  12 

Taxi 

t  Point  5  f  r or.  Pad  i 

Taxi 

5b 

41 

21 b .00 

2 16. 10 

2  1 . 9  b 

l  a  y  1 

fi  x  Pt.  5- Last  Chance 

I'nxi 

bb 

b 

— 

1  66 . 8(1 

2  2.96 

[  ,1X1 

f  re:  Pad  \i  »o  P  ml  21 

Last  Chance 

40 

10 

89.00 

106.01 

6  1.'’ 

Preflight  Check 

4  b 

3.. 

49.5" 

101.30 

1  /  ,  <.  4 

Taxi 

4  4 

71.00 

22.48 

i.  56 

’)  f  >.  i 

.'Ml'  Fl.'WflV 

Kngtiu*  Check 

100 

4  t 

10.00 

Vi.  14 

12.05 

tt  i  i  i 

a  r  v 

Runway  Ho l 1 

ion 

44 

21.00 

70.58 

1  .99 

M  i  1 1 

.  i  r  v 

In'  1  lehr 

C t  tmbout 

inn 

... 

— 

-- 

-- 

Approach 

HO 

--- 

— 

-- 

-- 

I  nbiuiml 

landing  Roll 

AO 

bb 

124 .00 

121  .  34 

i  /  .  38 

1’axl 

55 

b(l 

**4.  SO 

I’S.nn 

:..v  i 

c .  lv*‘  r  ?  1 1  ghr 

Host f light 

40 

b  7 

64.00 

72.  6' 

36.  -.2 

Taxi 

bb 

4  1 

119.50 

129.14 

19.6.4 

Taxi 

Inh'  ‘mil  (i!«ni* t a  1  ) 

Taxi 

ib 

9 

164.00 

164.67 

18  .  •.  7 

i'.m 

t  Pad  At  *  a  f'Oone  r  i  i 

1  axi 

bb 

n 

104.00 

JfH.02 

21  .  6  ’ 

Taxi 

t*.  1'ad  I1 

Taxi 

bb 

14 

15b. 50 

158.  36 

32. 37 

Taxi 

t.*  Pad  C 

Shut  down 

40 

S 

VI .  sn 

9  l  .00 

5.9  1 

T  -  38 

At rcraft 

Out  haund 

Idle  at  Startup 

4b 

28 

2  70 . 00 

282.06 

6  .  ‘lb 

Taxi 

60 

40 

120.00 

111. 4b 

42  ,  7  6 

Mineral  Taxi  to  Point  5 

Taxi 

60 

lb 

414.00 

421 .80 

7  1. 87 

!  a.x : 

to  ha>t  i  harue  (  5  ner  i 

Taxi 

60 

s 

— 

79  .00 

is.  r 

l  a  x  i 

to  Piiint  *  f  i  v:  -  5  a  J.  11 

Taxi 

60 

Is 

102.56 

10b. 79 

25.57 

T.  si 

tc  lVmr  ■'  ‘  t  ■-  pad  i 

Taxi 

60 

lb 

1 17.00 

144.27 

32. 2  7 

Taxi 

•  i*  :  -*f  '■  :  r.  -  T.i  1 

Taxi 

bO 

7 

-- 

185.29 

2  6 .  "  ! 

1  a.-M 

t  *  ;  vi:i!  5  f  r  ?r  !\i  :  K 

Taxi 

60 

22 

265.00 

260.  r 

28.  - 

1 .1X1 

f  r-1-  Pt .  b-L^st  Chance 

hast  Chance 

bn 

49 

?  1 . 00 

74.13 

12.7  1 

Taxi 

60 

10 

10 . 00 

34.03 

1  6  .  S  1 

Tax: 

t--*  Pref  leht 

P re fl  i gh t 

50 

29 

96.00 

108. 16 

7  3.-4  0 

l’axi 

60 

46 

2  4.00 

2  4 .  V 

6.6  3 

I  ix; 

■nLc  Sunwav 

Engine  (‘beck 

100 

46 

2  7 . 00 

32.5  3 

16. 54 

Mi  1  i  tar  -' 

Runway  Roll 

1 00+Af terburner 

18 

I9.no 

19  n4 

2 .  ’2 

Mil: 

irv  and  At  ri  rK'.n.’  v 

Inf  1 i ght 

<  i  unbent 

10n+Af terburner 

— 

-- 

-- 

— 

V  i  I  i 

irv  m  :  After  hir^r 

Approach 

SR 

— 

-- 

-- 

~ 

Inbound 

( .inding  Hr>  1 1 

bb 

3; 

89.00 

87.50 

16.76 

lax  i  <30R) 

60 

li 

90.00 

H9  .  Q  1 

•  0  ,  T9 

Taxi 

tc  !V  ;t  flight 

Taxi  (JOC) 

60 

2  l 

bl  .00 

51.  74 

J  7  .  30 

lax  l 

to  }’<•••  t  f  1  Ighl 

Host  f  1  i ghr 

bO 

40 

45.00 

4  « .  60 

M  ,6  1 

Tax  i 

bO 

46 

402.50 

391.03 

79 . 6  1 

I  ax  l 

Inbound  (4o:i-r.i!  1 

Taxi 

bt) 

19 

410.00 

386. 79 

'‘■’.0  1 

I  a  x  i 

to  Pad  t: 

Taxi 

bn 

7 

-- 

347.00 

1  1  ■  .  6  1 

iaxt 

to  laii  M 

Taxi 

bO 

4 

-- 

36  1  .  b() 

2  I.S  3 

l'axi 

to  Tad  I 

Tax  3 

bo 

-- 

4 1 1.00 

413.  37 

b  7 . 60 

T  a  x  l 

to  Pad  7 

Taxi 

bO 

446. 50 

44  3.67 

58.  Jh 

Taxi 

t  '  Pad  »* 

F— b 

Aircraft 

Out  bound 

Idle  at  Startup 

bO 

16 

2 4 ".60 

26b. 19 

io_qi 

i  a  x  i 

6  b 

tt) 

U-..00 

)  )6  .  8  7 

1  .1X1 

Out  h.  md  0  -em  i-'ii 

’axi 

bb 

2  l 

79  .  'll! 

8b  .  09 

2  1  . 

1'  a  x  l 

f  r  *r  Pal  t  •  •  1 

I  ix i 

6  b 

4 

7 ib.no 

.’r*l  .  78 

If, .  a  5 

Segment  6  t  .•  We  a:  -• 

^Papons  Arm-  ‘ 

' a s t  Chance 

bO 

1  1 

107.no 

161 .64 

16 5.42 

1  ex  i 

6  b 

n 
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operations,  are  the  most  highly  contaminated  d is t r ibu t uious .  Such  contami¬ 
nated  distributions  create  problems  in  choosing  a  spatial  distribution  or  most 
typical  value  to  associate  with  the  t imo- i n-mode  . 

t'reflight  and  postflight  modes  for  the  F-5  aircraft  exhibited  the  largest 
standard  deviations.  The  extent  of  contamination  oi  Lhe  observation-,  ol 
t  iir.e- 1  n'-mode  depends  on  the  number  or  aircraft  in  a  flight,  the  ground  service 
ere w,  pilot  proficiency,  and  whether  anning  or  dearming  is  to  be  perform,  j. 
This  range  illustrates  the  fact  that  there  may  be  more  than  one-  1  ndepenu.-;  : 
mechanism  functioning  in  an  aircraft  operational  mode. 

a.  DK TAILED  AIRCRAFT  ACTIVITY  AT  W1 1, LIAMS  AIK  FOKCK  BASE  F('R  THE 
PERIOD  1  JUNE  1976-30  JUNE  1977 

The  aircraft  activity  information  transcribed  from  flight  forms  A !  3:>5, 

ATC  90,  and  TWI2A  consists  primarily  of  the  dates  and  times  oi  actual  takeoffs 
and  landings  of  all  aircraft  making  use  of  the  Williams  Air  force  Base  (WAFH) 
facilities.  Form  ATC  355,  the  Runway  Supervisory  Unit  (RSI!)  hog ,  pertains  : . 
operations  involving  successful  takeoff  or  landing  of  local  aircraft;  tow.  ATC 
90.  the  Aircraft  Traffic  Log,  pertains  to  the  activities  ot  transient  air- 
era!  t;  and  Form  IV124,  the  Flight  Deviation  Report,  pertains  to  operations 
which  involve  deviations  from  scheduled  takeofl.  These  Air  Force  flight 
forms,  used  during  the  air  pollution  monitoring  program,  were  the  most 
reliable  source  on  which  to  base  hourly  emissions  activity  data  needed  for  the 
Air  Quality  Assessment  Model  (AQAM)  evaluation  effort. 

The  completed  transcription  process  provided  the  a ircraf t  activity  data 
base,  which  consists  of  the  daily  records  described  in  Table  D-3  where  each, 
day's  activity  begins  with  one  Head-  x  giving  fin  war,  month ,  day,  and 

fnl  I  traffic  count  on  each  airbase  runway.  Any  uuirher  of  Activity  Cards  ma; 
follow  the  Header  but  each  one  must  restate  the  date,  specify  only  one  of 
five  possible  operations,  denote  the  runway  in  use,  give  lhe  aircraft  type 
i evolved,  and  list  the  military  times  of  the  opetation. 


a.  Aircraft  Mix  at  Williams  Air  Force  Base 

A  review  of  the  activity  data  base  for  the  period  1  Tune  1976-30  June 
1977  shows  that  aircraft  at  WAFB  involves  scheduled  and  nonsctiedu i ed  flight 
operations  of  approximately  50  different  aircraft  types.  While  most  aircraft 
ate  small  twin-engine  trainer  types,  others  are  helicopter  types  and  largti 
cargo  types  equipped  with  either  turboprop  engines  or  jet  entities. 

Table  R-  4  provides  a  monthly  summary  ...  f  a  ire  is  ft  types  using  WAFB 
.'•ting  the  study  period.  This  in  format  ion  also  appears  1  at  er  as  part  >i 
t!ie  Monthly  Activity  Summary  tound  <n  Attachment  H. 


b.  Cla-s  i  f  i  cat  ion  of  Aircraft  Types  at  Williams  Air  Tone  t'.'-ic 


For  summary 
similar  aircraft 


report  purposes, 

uie  grouped 


sev 
■.  b  e 


<  ■  n 

i  r 


a  .  ’  raid  rj  ism  ’  w  >  1  ■  emerge  when 
• :  -.y  ■  ne  -  cl''  i  •-  i  ■  n  <  h  'i  i  at  t  *■  r  i  >  '  :  s  . 


f  .  * 


Table  D-3 .  Aircraft  Activity  Data  Format 


Header 

Card  -  One  Per  Day 

Col . 

Content 

1-6 

YYMMDD  = 

Year,  month,  day 

7 

C  => 

Total  aircraft  counts 

8-U 

30L 

l 

CM 

r—4 

30  C 

16-19 

30R 

► 

Full  traffic  count  on  these 

20-23 

12R 

Right  justify  integers 

24-27 

12C 

28-31 J 

12L 

Aircraft  Activity  Cards 

Col . 

Content 

1-6 

YYMMDD 

7 

N  => 

Normal  operations  -  times  are 
local  takeoff  and  landing  times. 

Slashes  separate  sorties. 

I  => 

Transient  inbound  1  Times  are  in  GMT 

/X 

H 

O 

Transient  outbound J  d*e.>  Zulu) 

< 

(1 

Deviation  -  times  are  local  engine 
start  and  takeoff  times. 

Slashes  separate  sorties. 

/\ 

II 

< 

Non-weather  aborts  -  times  are  local 
engine  start  times 

8-10 

Runway 

11-20 

Aircraft 

type  -  left  justified 

21-24  ' 

\ 

26-30 

1  Local  24- 

hour  times  separated  by  commas  or 

. 

1  slashes  as  appropriate,  (Zulu  time  for 

, 

1  transient 

aircraft. ) 

76-79. 

1 
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Table 

D-4 . 

Monthly 

Summary  of 

Ai  rcra  ft 

Mix  at  WAFB 

Month 

Aircraft  Type 

June  1976 

K5 

F4 

T37 

T38 

KF4P 

C 140 

A4 

A7 

T39 

A6P 

FI  00 

C 13  5 

A6 

HI 

AHl 

SH5 

OH  5  8 

H3 

UHl 

HH1 

Jul  y 

F5 

F4 

T37 

T38 

Cl  30 

C130A 

A4 

HI 

UHl 

H46 

F10 1 

VC  130 

A6 

T2 

T39 

T33 

C135 

KC135 

C4 

H3P 

TA4 

CH46 

TC4 

F14 

August 

F5 

F4 

T37 

T38 

C130 

F100P 

A4 

A3 

A37 

T33 

F105 

T2 

E2 

UHl 

T39 

Cl  35 

C9 

September 

F5 

F4 

T37 

T38 

Cl  30 

A21P 

A4 

T2 

T33 

T39 

FI  01 

FIDO 

C9 

HI 

UHl 

H53 

CH46 

CH53 

A7 

H46 

T28 

OH58 

Oc  t  ober 

F5 

F4 

T37 

T38 

T39 

Cl  30 

C9 

A7 

UHl 

TA4 

RF4 

C140 

H2 

A4 

A37 

T33 

T28 

VC  140 

HI 

T2 

A6 

A48 

November 

F5 

T37 

T38 

TRANHEL 

TRANTRN 

TRANRCP 

December 

F5 

F4 

T37 

T38 

300C 

BEECH 5 5 

A7 

A4 

1A1 

TA4 

C13  5 

BARON 

H3 

HI 

T39 

CH3 

OVIO 

F10 1 

A6 

T33 

AV8 

OUlO 

F100 

F14 

UHl 

OH  5  8 

January  1977 

F5 

T37 

T38 

TRANHEL 

TRANTRN 

TRANRCP 

A4 

T39 

February 

T37 

T38 

T39 

C130 

F5 

A4 

F100 

HI 

F4 

UHl 

T33 

A3 

OVIO 

A7 

T43 

C140 

1357 

T28 

T2 

H53 

AV8 

CF5 

March 

T37 

T38 

F5 

T39 

A4 

F4 

F106 

A7 

T33 

T2 

UHl 

OVIO 

FI  00 

C130 

C4 

TC4 

Apr  l  l 

T3  7 

T38 

F5 

A4 

AV8 

F4 

H58 

T39 

OH  5  8 

H3 

A3  7 

HI 

A6 

C130 

T2 

UH1 

H46 

Cl  3  5 

BEECH  55 

F10 1 

OV10 

T33 

C140 

Table*  D-4 .  (Cont'd) 


Month  Aircraft  Type 


May  1977 

T37 

T38 

F3 

Cl  30 

02  1 

T33 

OVIO 

F101 

A4 

HI 

UH1 

F4 

F100 

C9 

T41 

VC  140 
H46 

H53 

T39 

A7 

Cl  3  5 

June 

T37 

T38 

F5 

H46 

A4 

T39 

HI 

CESSNAl 

T33 

F4 

C130 

C140 

H3 

OV10 

Ab 

Classes  1,  2,  and  3  will  be  used  to  designate  aircraft  types  F3,  T37,  and  T38, 
respectively.  Classes  4,  5,  and  6  will  be  used  to  designate  three  classes  of 
transient  aircraft,  comprising  aircraft  types  that  are  similar  to  F4,  C130, 
and  C141,  respectively,  where  F4  represents  a  fighter  type  class,  C130  a 
turboprop  engine  class,  and  C141  a  multienginc  jet  class.  Class  7  will  b, 
reserved  for  transient  aircraft  of  the  helicopter  type. 

For  AQAM  evaluation  purposes,  only  the  first  six  classes  of  aircraft  will 
be  considered;  class  7,  the  helicopter  types,  will  be  omitted. 

Table  D-5  shows  a  listing  of  transient  aircraft  types  and  their  respec¬ 
tive  classification.  Note  that  all  aircraft-type  notations  are  as  given  after 
the  information  transcription  process. 


c.  Definitions  for  Aircraft  Activity  Operations  at  VAFB 

Aircraft  activity  data  compiled  for  WAFB  during  the  period  1  June  ]97fo-3u 
Juno  1977  provides  an  accounting  for  aircraft  operations  that  were  normal, 
transient,  and  aborted  or  otherwise  deviated  from  takeoff. 


(1)  Normal  Operations 

Normal  operations  are  defined  for  aircraft  types  F5 ,  T37,  and  T38 

only.  Local  military  time  was  recorded  at  the  moment  of  liftoff  (during 
takeoff  mode)  and  again  at  the  moment  of  touchdown  (during  the  landing  mode). 
For  landing  operations  involving  touch-and-go  maneuvers,  the  military  time 
recorded  was  the  moment  of  last  touchdown.  Whether  some  Air  Force  flights  did 
or  did  not  originate  at  WAFB,  their  takeolf  or  landing  time  was  reported  in 
this  normal  group. 


(2)  Transient  Operations 

Transient  operations  are  defined  for  all  aircraft.  Greenwich 
Mean  Time  (Zulu)  was  recorded  for  each  inbound  and  each  outbound  transient 
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Table  D-5.  Classification  of  Transient 
Aircraft  at  WAFB 


Class  4 

Class  5 

Class  6 

Class  7 

A4 

A3 

C4 

AH  7 

A6 

Beech  55 

C9 

CH3 

A7 

Beech  Baron 

C135 

CH46 

A37 

C130 

C140 

CH56 

AV8 

C130A 

C141 

HI 

B57 

Cessna  1 

CF5 

H2 

F4 

Cessna  300 

KC135 

H3 

F14 

E-2 

T2 

H46 

F100 

OV-10 

TC4 

H53 

F101 

T-28 

VC140 

H58 

F105 

T-41 

HH1 

F106 

T-43 

OH58 

RF4 

TA-1 

SH5 

Tee 

TRANRCP 

U10 

T39 

300C* 

U21 

TA4 

UH7 

TRANTRN 

TRANHEL 

A48* 

H3P* 

A6P* 

A21P* 

F100P* 

RF4P* 

♦This  aircraft  type  is  questionable:  possible  data 
recording  error. 


operation.  These  operations  at  WAFB  account  for  nonscheduled  aircraft 
arrivals  and  departures. 


(3)  Abort  Operations 

Abort  operations  are  defined  for  aircraft  types  F5,  T37,  and  T38 

only.  Local  time  was  recorded  at  the  moment  of  aircraft  engine-start.  This 
is  the  only  recorded  information  in  regard  to  nonweather  abort  operations. 


96 


(4)  Deviation  Operations 

Deviation  from  scheduled  takeoff  is  an  operation  defined  for 
aircraft  types  F5,  T37,  and  T38  only.  Local  time  was  recorded  if  the  time  at 
engine-start  deviated  more  chan  5  minuter,  from  scheduled  start  time  and  again 
if  the  time  at  takeoff  deviated  more  than  10  minutes  from  scheduled  takeoff. 
Upon  completion  of  the  takeoff  procedure,  the  deviation  operation  becomes 
identical  to  a  normal  operation;  only  exceeding  it  in  duration  by  several 
minutes  of  ground-time.  The  total  ground-time  being  the  full  time  interval 
between  engine-start  and  final  liftoff 

d.  Edits  to  the  Aircraft  Activity  Data  Base 

Some  of  the  entries  resulting  from  the  transcription  process  were  found 
to  be  erroneous  and  were  resolved  by  making  a  best-guess  estimate  of  actual 
time  values.  For  example,  some  normal  operation  time  entries  (sorties)  would 
extend  beyond  the  hour  of  midnight,  while  other  takeof 1/ land ing  times  for 
sorties  seemed  transposed. 

In  the  cases  where  time  entries  would  extend  beyond  the  midnight  hour, 
the  sortie  was  assumed  to  begin  and  then  end  on  consecutive  days.  This  was 
achieved  by  substituting  the  time  -1  (i.e.,  unspecified)  as  the  landing  time, 

then  adding,  for  the  next  day,  a  sortie  whose  takeoff  time  was  -1  and  landing 
time  was  used  from  the  previous  record. 

In  the  cases  where  sortie  times  seemed  to  be  transposed,  the  time  entries 
were  interchanged. 

(1)  Edits  to  Normal  Operations  Data 

Table  D-6  provides  the  list  of  edits  made  to  the  normal  operations 
portion  of  aircraft  activity  data  for  WAFB.  The  date,  runway,  and  aircraft 
type  are  given  in  columns  1,  2,  and  3  respectively.  Uncorrected  takeoff  and 
landing  times  are  shown  in  columns  4  and  5  beside  the  corrected  times  shown  in 
columns  6  and  7.  The  indicator,  -1,  shows  a  nonoperation  so  these  are 
not  included  when  the  operations  are  tallied. 

(2)  Edits  to  Deviation  Operations  Data 

Table  D-7  provides  the  list  of  edits  made  to  the  deviation  opera¬ 
tions  portion  of  aircraft  activity  data.  The  date,  runway,  and  aircraft  type 
are  given  in  columns  1,  2,  and  3,  respectively.  Uncorrected  aircraft  engine- 
start  and  takeoff  times  are  shown  in  columns  4  and  5  and  their  corrected  times 
appear  in  columns  6  and  7. 

e.  Full  Traffic  Count  on  Airbase  Runways 

The  layout  of  airbase  runways  is  shown  in  Fig.  D-2.  Runways  prefixed  by 
a  12  are  used  only  when  tin*  northwesterly  wind  exceeds  10  knots.  The  letters 


Table  D-6.  Edits  to  Normal  Operations  Data  at  WAFB 


Date 

Knwy 

Aircraft 

Type 

Uncurrec t  ed 

Co  r  r  t 

T 

c  ted 

L 

760630 

300 

F5 

1600 

15  54 

1600 

1654 

760702 

30C 

T38 

-1 

0043 

760702 

300 

T38 

-1 

0058 

760702 

30C 

T38 

0005 

0057 

760701 

30C 

T38 

2343 

0043 

2343 

-1 

760701 

30C 

T38 

2350 

0058 

2350 

-1 

760701 

30C 

T38 

0005 

0057 

760702 

30C 

T38 

-1 

0008 

760702 

30C 

T38 

-1 

0018 

760702 

30  C 

T38 

-1 

0030 

760702 

30C 

T38 

-1 

0024 

760701 

30C 

T38 

2308 

0008 

2308 

-1 

760701 

30  C 

T38 

2318 

0018 

2318 

-l 

760701 

30C 

T38 

2327 

0030 

2327 

-1 

760701 

30C 

T38 

2330 

0024 

2330 

-1 

760720 

30C 

T38 

2311 

0007 

2311 

-1 

760721 

30C 

T38 

-1 

0007 

760720 

30C 

T38 

2306 

0005 

2  306 

-1 

760720 

30C 

T38 

2308 

0007 

2308 

-1 

760721 

30C 

T38 

-1 

0005 

760721 

30C 

T38 

-1 

0007 

760712 

30C 

T38 

2322 

0019 

2322 

-1 

760713 

30C 

T38 

-1 

0019 

760712 

30C 

T38 

2311 

0006 

231  1 

- 1 

760713 

30C 

T38 

-1 

0006 

760712 

30C 

T38 

2308 

0004 

2308 

01 

760713 

30C 

T38 

-1 

0004 

760712 

30C 

T38 

2342 

0032 

2  342 

-1 

760713 

30C 

T38 

-1 

0032 

760712 

30C 

T38 

2314 

0001 

2314 

-1 

760713 

30C 

T38 

-1 

0001 

760712 

30C 

T38 

2335 

0036 

2  335 

-1 

760713 

300 

T38 

-1 

0036 

760819 

30C 

T38 

2312 

0006 

2312 

-1 

760819 

300 

T38 

2314 

000  5 

2314 

-1 

760820 

300 

T38 

-1 

0006 

760820 

300 

T38 

-1 

0005 

760818 

30  C 

T38 

2  122 

0004 

2  122 

-1 

760818 

300 

T38 

2  130 

001  1 

2  130 

-1 

760818 

300 

T38 

2330 

0013 

2  130 

-1 

760819 

300 

T38 

-1 

0004 

760819 

300 

T38 

-1 

001  3 

760819 

300 

T  38 

-1 

001  3 

Tab  If  D-6  .  fCnnt'd) 


Date 

;<nwy 

Aircraft 

Type 

One  jr 

T‘> 

r  t>c  t  »■  d 

Lh 

<,.>r  rn< 

T 

t  ed 

L 

760812 

30C 

T38 

231  7 

0074 

2  317 

- 1 

760812 

300 

T38 

2329 

0026 

2329 

-i 

760812 

300 

T38 

2334 

0026 

2  3  34 

-i 

760813 

30C 

T38 

-1 

0024 

760813 

300 

T38 

-1 

002  6 

760813 

300 

T38 

-1 

0026 

760811 

300 

T38 

1909 

16  21 

15  09 

1621 

760802 

300 

T  38 

2304 

000a 

2  <04 

-1 

760803 

300 

T38 

-1 

0004 

761108 

30K 

T38 

1311 

1252 

1211 

1252 

761108 

30  L 

T37 

1120 

024  7 

1  120 

1247 

761201 

30L 

T37 

1836 

1  728 

1836 

1928 

761203 

301, 

T3  7 

1640 

1451 

1  340 

1451 

761210 

300 

T38 

1509 

1410 

1309 

1410 

770112 

300 

T38 

1617 

1533 

1617 

17  33 

770120 

300 

F5 

1513 

1317 

11)3 

1317 

770222 

300 

T38 

2224 

1907 

2224 

-1 

770301 

300 

F5 

1303 

0901 

1303 

1401 

770301 

300 

F5 

1750 

0849 

1750 

1849 

770307 

30C 

T38 

1403 

12  3  7 

1203 

12  37 

770308 

30C 

T38 

1329 

0442 

1329 

1442 

770308 

300 

T38 

1249 

0620 

1249 

1620 

770314 

300 

R-5 

1620 

1408 

1620 

1708 

770331 

300 

T38 

1805 

1351 

180  5 

1851 

770408 

300 

T38 

1233 

1125 

1033 

1125 

770408 

30C 

T38 

1233 

1125 

1033 

1125 

770415 

30L 

T37 

1939 

1058 

09  19 

1058 

770504 

30C 

T38 

1744 

1003 

0844 

100  1 

770523 

300 

F5 

1  1  32 

1013 

1  132 

1213 

770607 

300 

T  3  8 

1  2  59 

1062 

1  2  59 

14  02 

770614 

300 

F5 

082  2 

0000 

082  2 

-1 

770620 

30C 

T38 

3  34  3 

004  7 

3  34  3 

-J 

776020 

300 

T38 

?  34  3 

0026 

2  34  3 

-1 

776020 

300 

T38 

0003 

0049 

770620 

300 

T38 

0003 

0036 

770621 

30  C 

T38 

-1 

004  / 

770621 

300 

T38 

-1 

0026 

770621 

30C 

T38 

0003 

(.1049 

770621 

300 

T38 

0003 

0036 

770629 

30C 

T38 

2228 

0000 

2228 

-1 

aT  “  Takeoff  Time 
^L  =  Landing  Time 
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Table  D-7.  Edits  to  Deviation  Operations  Data  at  WAFB 


Date 

Rnwy 

Aircraft 

Type 

Oncorrec  ted 

ESa  Tb 

Correc  ted 
h'S  T 

760621 

30C 

T38 

0654 

0619 

0604 

0619 

760818 

30C 

T38 

1658 

1620 

1658 

1720 

760804 

30C 

T38 

1723 

1540 

1723 

1740 

760901 

30C 

T38 

1439 

1352 

1339 

1352 

760902 

30C 

T38 

1545 

1401 

1545 

1601 

760929 

30C 

T38 

1656 

1620 

1606 

1620 

761027 

30C 

T38 

1630 

1605 

1605 

1630 

761027 

30C 

T38 

1334 

1251 

1251 

1334 

761103 

30R 

T38 

1555 

1514 

1455 

1514 

761109 

30R 

T38 

1510 

1447 

1410 

1447 

761122 

30R 

T38 

1545 

1403 

1345 

1403 

761123 

30L 

T37 

1558 

1402 

1558 

1602 

761215 

30C 

T38 

1156 

1116 

1056 

1116 

770201 

30C 

T38 

1147 

1101 

1107 

1121 

770519 

30L 

T37 

1717 

1628 

161  7 

1628 

770523 

30C 

T38 

0750 

0728 

0728 

0750 

770601 

30C 

T38 

— 

1332 

1312 

1322 

770614 

30L 

T37 

— 

1556 

1546 

1556 

770620 

30C 

T38 

2130 

1249 

2130 

2149 

770624 

30L 

T37 

1327 

0357 

1327 

1357 

770628 

30C 

T38 

1630 

1045 

1630 

1645 

aES  =  Takeoff  Time 


bL  ■  Landing  Time 


L,  C,  and  R  denote  left,  center,  and  right,  respect ive ly ,  Runway  30L  is  used 
almost  exclusively  by  T37  aircraft;  runways  30C  and  30R  are  used  jointly  by 
T38  and  F5  aircraft. 

The  full  traffic  count  on  a  runway  is  determined  as  the  sum  of  landing 
plus  takeoff  operations,  plus  twice  the  number  of  touch-and-go  operations  on 
that  runway.  It  is  unfortunate  that  these  daily  runway  traffic  counts, 
gathered  during  the  WAFB  experimental  period,  were  inaccurate  because  of  the 
several  methods  used  to  record  aircraft  operations.  Specifically,  three 
methods  exist  by  which  runway  traffic  counts  were  recorded,  thus  the  data 
base,  in  effect,  is  partitioned  into  three  subsets,  in  regard  to  daily 
traffic  counts:  (l)  the  actual  count  or  the  daily  averagp  based  on  estimated 
monthly  activity,  depending  on  the  runway,  (2)  the  actual  count  on  each 
runway,  and  (3)  the  actual  count  or  the  daily  average  based  on  actual  monthly 
activity,  depending  on  the  runway. 


(I)  Fa  L  l  Traffic  Count  for  June-  and  July  1976 

The  methods  used  to  record  June  and  July  1976  full  traffic  counts 
are  runway  dependent.  For  runway  30L  (and  12R),  which  are  used  almost 
exclusively  by  T37  aircraft,  the  actual  daily  activity  count,  is  recorded. 
Traffic  counts  for  runways  30C  (and  12C)  and  30R  (and  1 2  L ) ,  used  jointly  by 
T38  and  F5  aircraft,  are  average  daily  traffic  counts  derived  from  estimated 
monthly  activity. 


(2)  Full  Traffic  Count  for  August  1976-Januarv  1977 

For  the  6-month  period,  August  1976-January  1977,  the  traffic 
count  for  each  runway  is  given  as  the  actual  daily  activity  count. 


(3)  Full  Traffic  Count  for  February-June  1977 

The  methods  used  to  record  February-June  1977  full  traffic  counts 
are  runway  dependent.  For  runways  30L  and  12K,  used  almost  exclusively  by  T37 
aircraft,  the  counts  were  supplied  by  the  Flight  Records  Branch,  WAFB,  and 
represent,  for  most  days,  the  actual  daily  activity  count.  However,  there  are 
several  days  for  which  the  traffic  count  is  an  estimated  value. 

The  traffic  counts  for  runway  30C,  used  jointly  by  T38  and  FI 
aircraft,  were  recorded  on  AFCS  Form  5a  (Control  Tower  and  CCA  logs).  The  use 
of  the  Tower  and  CCA  logs  provided  an  estimate  of  daily  traffic  counts  since 
the  Tower,  in  addition  to  controlling  the  center  runway  (30C),  always  assumes 
control  of  other  runways  whenever  shutdowns  of  respective  KSUs  occur.  During 
shutdown  periods,  no  distinction  is  made  between  runway  use  by  aircraft. 

The  traffic  counts  for  runway  30R,  used  jointly  by  T3  8  arid  F5 
aircraft,  were  supplied  by  the  Air  Traffic  Control  Superintendent  at  WAFB. 
Daily  traffic  counts  were  derived  by  equally  distributing  the  supplied  monthly 
totals  of  RSU  station  name  GASSER  among  the  weekdays  of  the  month,  i.e.,  both 
weekends  and  holidays  were  excluded. 


f.  Touch-and-go  Training  Flights 

Touch-and-go  training  flights  are  carried  out  by  student  pilots  using 
T37,  T38,  and  F5  aircraft.  In  performing  these  maneuvers,  an  aircraft 
approaches  and  lands  on  a  runway,  travels  down  the  runway  for  several  seconds, 
accelerates,  and  lifts  off.  The  runway  in  use  depends  on  the  aircraft 
involved.  Runway  30L  is  used  exclusively  by  T37  aircraft  and  runways  30C  and 
30R  are  shared  by  T38  and  F5  aircraft. 

Unlike  the  hour  of  landing  and  hour  of  takeoff  described  tor  normal 
operations,  we  deline  an  hourly,  discrete  time  history  of  touch-and-go  opera¬ 
tions  specific  to  runway  and  aircraft  type  as  a  function  of  the  total  number 
of  operations  (landing  and  takeoff)  on  that,  runway  by  the  relation: 
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TG .  . 


1  . 

2  I 


Full  Traffic  Count 


where : 

HAij 


HT 


FTC 


N  24 


Y  V  (HA. .  +  HT.  .  )' 

j=!  A  lJ  1J  j 


the  numoer  of  landing  operations  on  the  runway  by 
aircraft  type  j  during  hour  i, 

the  number  of  takeoff  operations  on  the  runway  by 
aircraft  type  j  during  hour  i,  and 

the  daily  total  of  landings,  takeoffs,  plus  twice 
the  touch-and-go  operations  on  the  runway. 


In  consideration  of  the  Full  Traffic  Count  as  derived  by  the  several 
methods  stated  above,  it  is  realized  that  the  imbedded  touch-and-go  operations 
specific  to  runway  30C  or  runway  30H  do  not  distinguish  between  the  aircraft 
types  F5  or  T38.  Thus,  all  or  none  of  these  operations  could  be  performed  by 
either  aircraft  type.  Therefore,  we  note  that  both  F5  and  T36  aircraft  are 
equipped  with  two  J85  aircraft  engines  and  that  all  significant  operational 
parameters,  such  as  their  altitude  at  final  phase  of  approach;  their  speed  at 
point  of  touchdown,  etc.,  are  much  the  same.  It  can  be  concluded  that  the 
above  relationship  is  sufficient  for  this  application  of  the  AQAM,  though  it 
may  not  suffice  for  more  detailed  studies,  such  as  air  quality  impact  analysis 
of  individual  aircraft  types. 


g.  Summary  Report  of  Aircraft  Activity  at  WAFB  for  the  Period 
1  June  1976-30  June  1977 

In  order  to  demonstrate  the  type  of  activity  data  available  for  input  to 
the  Air  Quality  Assessment  Model,  a  computerized  program  was  developed  that 
summarizes  the  Aircraft  Activity  Information  collected  at  WAFB,  Arizona.  The 
input  data  to  the  program  is  formatted  as  shown  in  Table  D-3;  resulting  sum¬ 
maries,  formatted  with  different  headings  but  identical  bodies,  are  produced 
for  the  following: 

(i)  Hourly  Aircraft  Activity  (A  one  day  example  summary 
is  presented  in  Attachment  A), 

(ii)  Daily  Aircraft  Activity  (not  presented  here), 

(iii)  Monthly  Aircraft  Activity  (Attachment  B) , 

(iv)  Grand  Totals  of  Aircraft  Activity  (Attachment  C). 

Summaries  (i)  and  (ii)  are  headed  by  each  date  having  a  day-label  such 
as,  weekday  (for  Monday  through  Friday),  weekend  day  (for  Saturday  or  Sunday), 
weekday-holiday,  or  weekend-holiday.  it  is  instructive  to  note  that  a  holiday 
pattern  is  realized  when  very  low  aircraft  activity  periods  are  examined.  The 
following  table  illustrates  the  observed  pattern  tor  the  dates  that  reported 
no  aircraft  activity  at  all. 
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Date  s 

Hoi iday 

Day  of  Week 

1 

July 

19  76 

4 

duly 

19  76 

Independence  bay 

Sunday 

5 

du  1  y 

1976 

6 

Sept 

.  1976 

Labor  Day 

Monday 

1 1 

Gc  t  . 

1976 

Columbus  Day 

Monday 

25 

Oc  t . 

1976 

Veterans  Day 

Monday 

24 

Dec  . 

1976 

25 

Dec. 

1976 

Christmas  Day 

Saturday 

26 

Dec  . 

1976 

31 

Dec  . 

1976 

1 

Jan . 

1977 

New  Year's  Day 

Saturday 

2 

dan  . 

1977 

29 

May 

1977 

30 

May 

1977 

Memorial  Day 

Monday 

Thanksgiving  Day,  Thursday,  25  November  1  976,  although  a  holiday,  was 
found  to  have  one  outbound  transient  fightet  operation. 

In  the  body  of  each  summary  we  see  that  the  column  "Full  Traffic  Count" 
is  an  aggregate  of  all  aircraft  types,  but  that  each  of  the  other  columns  have 
aircraft  type  indicators  with  the  corresponding  meaning  as  follows: 

F5  Class  1 

T37  Class  2 

T38  Class  3 

F  Class  4  (F4  -  Fighter  Class) 

T  Class  5  (C130  -  Turboprop  engine  clsss) 

M  Class  6  (C141  -  Multiengine  jet  class) 

H  Class  7  (Helicopter  class) 

The  repeated  indicators  are  then  grouped  by  the  mode  of  the  operation 
M and i ng- inbound ,  takeoff-outbound,  ground  time,  touch-and-go),  and  then, 
where  appropriate,  by  the  previously  defined  activity  operations  Normal, 
Deviation,  Abort,  or  Transient. 

The  frequency  of  occurrence  of  each  of  the  activity  operations,  in  mode 
of  the  aircraft,  type,  is  given  for  each  airbase  runway.  Hie  ground  time  of 
deviation  operations,  however,  is  given  in  minutes  on  the  runway. 

Six  runways  are  identified  from  which  all  airbase  aircraft  operations  are 
performed.  The  computer  program  is  coded  so  as  to  assure  that  the  following 
assumptions,  regarding  runway  use,  are  preserved: 

(1)  All  T37  operations  (landing,  takeoff,  and  touch-and- 
go)  are  performed  on  the  inside  runway  (301.  or  ]2R). 

(2)  All  C 1 30  and  all  C141  operations  (transient)  are 
performed  on  301,  (n.b.,  only  301,  can  handle  the 
heav  v  a i rc  r a f t  )  . 

(3)  All  T3H  and  K3  operations  are  performed  on  either  the 
center  runway  ( 30C  or  12C)  or  the  outside  runway  (30k 
or  121.)  with  the  shared  pattern  as  follows: 


(i)  All  takeoffs  (T38  and  F5)  occur  from  the  center 
runway  (30C  or  12C). 

(ii)  88%  of  T38  landings  occur  on  the  outside  runway 
(30R  or  12L),  the  remainder  occur  on  the  center 
runway  ( 30C  or  12C). 

(iii)  80%  of  F5  landings  occur  on  the  outside  runway 
(30R  or  12l),  the  remainder  occur  on  the  center 
runway  (30C  or  12C). 

Just  as  was  stated  above  about  the  holiday  pattern  of  aircraft  activity, 
we  note  here  a  weekend  pattern  of  aircraft  activity.  In  general,  there  are 
cases  of  weekend  airbase  aircraft  activity,  bet  these  events  are  not  recorded 
as  ar*  the  weekday  events.  Some  specially  requested  weekend  data  were 
supplied  covering  the  weekends  of  December  1976  and  January  1977  and  these 
data,  though  limited,  served  to  define  the  weekend  pattern  of  events  at 
WAFB . 
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ATTACHMENT  A:  HOURLY  AIRCRAFT  ACTIVITY  SUMMARY;  A  ONE-DAY  EXAMPLE,  6  AUGUST  1976 
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ATTACHMENT  C:  GRAND  TOTALS  OF  AIRCRAFT  ACTIVITY,  1  JUNE  1976-30  JUNE  1977 
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APPENDIX  E 


ESSAY:  VERIFICATION  OF  SHORT-TERM  AIK  QUALITY  MODELS 

USING  THE  GAUSSIAN  DISPERSION  APPROACH 

by  Karl  Zeller,  EPA* 

(January  1977) 


l.  INTRODUCTION 

A  need  exists  to  define  an  acceptable  approach  for  the  verification  of 
short-term  Gaussian  air  quality  simulation  models  for  single  and  multiple 
source  complexes.  The  demand  for  answers  to  air  quality  questions  is  so 
strong  that  some  models  have  been  formulated  and  applied  extensively  prior  to 
proper  evaluation.  It  is  hoped  that  proposed  models  would  be  subjected  to  a 
verification  procedure  to  enable  the  model  developer  either  to  defend  a 
particular  use  of  his  model  or  to  indicate  possible?  situations  in  which  his 
model  produces  invalid  or  questionable  results. 

Each  complex  Gaussian  air  quality  model  is  comprised  of  submodels  to 
describe  source  emissions,  pollutant  dispersion  and  transport,  plume  rise, 
source-receptor  geometry,  and  meteorology.  In  the  modeling  of  complex  geo¬ 
metric  situations,  it  is  currently  necessary  to  use  empirical  dispersion 
parameter  values  (standard  deviations  of  plume  spread  as  a  function  of  down¬ 
wind  distance  or  time  and  atmospheric  stability)  derived  from  continuous  point 
sources.  Sometimes  these  values  are  applied  with  subjective  modifications. 
For  example,  models  used  for  airports  or  highways  use  the  same  dispersion 

parameters  as  models  used  for  elevated  area  or  point  sources. 

One  logical  approach  would  involve  the  separate  verification  of  each  of 
the  submodels  based  on  specific  experiments;  however,  expenses  involved 
usually  dictate  that  accuracy  of  the  overall  model  be  examined. 

There  are  many  opinions  as  to  how  a  Gaussian  model  that  uses  subjec¬ 
tive  mod  i  f  ic  at  ions  to  handle  the  various  configurations  encountered  sb'-ui  ■:  re¬ 
validated,  verified,  or  calibrated.  In  this  discussion,  a  specific  approach 
for  Gaussian  models  is  suggested.  The  data  requirements  for  model  verifica¬ 
tion  include  meteorology  and  emission  data  (the  necessary  information  tor 

model  concentration  calculations),  and  measured  air  quality  concentrations 
(including  background  information)  with  which  to  compare  the  calculations. 

Previous  verification  programs  comparing  short-term  (1-3  hours)  cal¬ 
culations  (Koch  1971;  Rote,  et  al.,  1973)  with  observations  have  demonstrated 
the  difficulty  in  coming  to  specific  conclusions  using  any  one  comparison 

technique  or  statistic.  The  scatter  diagram  approach  prevalent  in  the 
literature  has  not  been  particularly  useful  for  verification  because  of 
wide  scatter  and  low  correlation  between  measured  and  calculated  concentra¬ 
tions.  The  reason  direct  comparison  of  calculations  to  observations  does  not 
provide  good  results  is  the  statistical  nature  of  the  complex  Gaussian  model. 
For  instance,  the  Gaussian  approach  assumes  steady-state  conditions  over  a 
specific  period  of  time,  usually  one  hour.  In  reality  a  continuous  point 
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source  plume  is  subjected  to  many  variations  of  wind  velocity  during  any  given 
hour ,  When  modeling  intermittent  or  moving  sources  such  as  automobiles  or 
aircraft,  simplifying  assumptions  are  necessarily  made;  for  instance,  taxiing 
aircraft  are  usually  modeled  as  line  sources  considered  to  be  continuous  over 
periods  relatively  long  compared  to  individual  aircraft  taxi  time.  In  real¬ 
ity,  there  is  not  a  continuous  Gaussian  plume  stretched  out  downwind  of  such 
sources  but  intermittent  ones  which  are  locally  distorted  by  variable  wind 
speed  and  direction  (turbulence).  Data  collection  for  such  comparisons  is 
presently  accomplished  by  a  network  of  monitors  at  fixed  horizontal  and 
vertical  locations.  The  total  number  of  monitors  or  stations  in  most  studies 
is  usually  less  than  five  due  to  monetary  restraints.  Considering  these 
limitations,  it  is  not  surprising  that  short-term  calculations  and  observa¬ 
tions  do  not  correlate  well.  The  comparison  of  cumulative  frequency  distribu¬ 
tions  on  the  other  hand  has  enabled  some  modelers  to  make  general  statements 
about  the  overall  performance  of  their  specific  models. 

2.  RECOMMENDED  VERIFICATION  APPROACH 

It  is  assumed  that  a  large  data  base  is  taken  over  a  period  long  enough 
(at  least  one  year)  to  provide  information  under  the  varying  conditions  of 
meteorology  and  emission  modes.  Assuming  such  a  data  base  is  available,  the 
following  procedures  for  model  verification  are  proposed. 

a.  Discussion  of  Data  Base 

A  discussion  of  the  data  base  should  include  efforts  made  to  provide 
background  information  and  quality  assurance  during  the  actual  collection  of 
data.  A  discussion  of  the  adequacy  of  the  numbers  and  locations  of  actual  air 
quality  monitoring  sites  should  be  included,  together  with  specific  geometry 
relating  sources  and  receptors  to  the  general  meteorological  conditions  during 
the  data  collection  period. 

b.  Comparisons  for  Analysis 

A  thorough  verification  requires  that  data  and  concentrat ion  calcula- 
t  ions  (model  predictions)  be  compared  under  a  number  of  different  circum¬ 
stances.  Because  of  the  conglomerate  nature  of  the  complex  Gaussian  model, 
certain  meteorological  conditions  or  receptor  locations  may  give  better  or 
worse  results  compared  to  others  in  the  same  model  application.  Each  monitor¬ 
ing  location  should  therefore  be  evaluated  separately  in  addition  to  overall 
evaluation  of  the  model.  Evaluating  each  monitoring  location  separately  will 
give  the  modeler  an  indication  of  the  performance  of  the  dispersion  submodel 
tinder  different  situations.  This  is  important  because  at  the  present  time  the 
dispersion  submodels  of  most  Gaussian  complex  models  have  been  altered  and  are 
subjective  in  nature  and  not  based  on  applicable  dispersion  experiments. 
There  is  reason  to  believe  that  a  great  deal  of  the  wide  scatter  in  observed 
and  predicted  values  already  discussed  is  due  to  the  dispersion  submodel 
as  well  as  faulty  emission  estimates.  The  pollution  data  should  he  stratified 
in  at  least  the  following  categories: 


High  and  low  emission  density  periods. 

Periods  of  the  day  (can  include  more  than  one  hour) 
dependent  upon  emission  operations  and  meteorology. 

Seasons  (this  category  is  optional  -  depends  on  the 
nature  of  the  data  and  problems). 

Meteorological  categories: 

wind  direction  -  increments  depend  on  situation. 

wind  speed  -  increments  depend  on  situation. 

stability  category  -  Pasquill  categories  o:  in  some 

cases  stable,  neutral,  and 
unstable  will  he  sufficient. 

mixing  height  -  two  or  three  categories  based 

on  a  chosen  mixing  height 
dependent  upon  the  scale 
of  the  model  applicat ion.  For 
instance,  in  the  case  of 
airport  models:  mixing  height 
below  100  m;  mixing  height 
above  100  m;  and  unlimited 
mix  ing . 

Note  that  the  above  categories  are  not  meant  to  be  mutually  exclusive. 
Also,  only  data  that  are  above  the  noise  level  of  the  particular  pollutant¬ 
monitoring  instruments  used  should  be  evaluated.  Although  there  is  a  specific 
interest  in  pollutant  levels  approaching  National  Ambient  Air  Quality 
Standards,  all  levels  should  be  evaluated. 

c.  Cumulative  Frequency  Distributions  and  Error  limits 

The  data  should  first  be  sorted  according  to  the  above  scheme;  for 
instance,  all  cases  in  a  specific  wind  direction  category,  all  cases  during 
the  morning  hours,  etc.  (approximately  26  separate  categories).  The  data 
should  then  he  displayed  on  cumulative  frequency  distribution  diagrams, 
similar  to  the  one  in  Fig.  K-l  ,  if  at  least  30  observations  within  a  givin 
category  are  available.  If  there  are  less  than  30  observations,  a  scatter 
diagram  of  observed  versus  calculated  concentrations  should  he  prepared. 
Longer  term  averages,  for  instance,  monthly  or  vearlv,  will  inherently  have 
fewer  data  points  for  comparisons  and  therefore  will  require  the  scatter 
diagram  approach  for  verification.  The  use  of  the  scatter  diagram  will 
give  a  qualitative  feel  for  the  performance  of  the  model  even  though  in  some 
cases  resultant  coefficients  for  calibration  purposes  may  bo  lacking  in 
s  ign  i  f  ic  anc  e  , 

It  is  desirable  to  have  a  number  measure  of  model  performance  for  spe¬ 
cific  cases.  A  cumulative  frequency  distribution  of  the  calculation  error 
should  he  constructed  to  allow  the  modeler  to  report  the  percentage  of 
trials  for  wh ich  his  estimate  was  within  a  specific  range  of  the  observation. 
A  similar  approach  was  used  by  Koch  and  Thayer  (1971);  frequency  distribut  ion 
plots  of  absolute  value  of  overpredictions  and  underpredictions  similar  to 
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those  in  Fig.  E-2  should  be  prepared.  The  approach  sup  re- ted  in  this  discus¬ 
sion  is  to  subtract  each  observation  from  each  corresponding  model  prediction. 
This  will  provide  information  on  the  distribution  oi  ovt rpred i ct ions  and 
unde rpred ic t ions  in  addition  to  the  overall  tendencv  of  the  model  to  over-  or 
under-predict  in  that  specific  category.  Subtracting  the  cumulative  percent¬ 
age  of  overpred ic t ions  from  the  cumulative  percentage  of  underprediction 
(see  Fig.  E-2  )  at  a  specific  concentration  difference  wi 1 i  give  the  percent  of 
comparisons  that  fall  within  that  specific  concentration  error  limit. 
The  error  limits  and  corresponding  percentages  can  then  be  plotted  or,  a 
frequency  distribution  diagram  similar  to  Fig.  !•-').  Tin-  modeler  will  then  be 
able  to  report  that  in  the  given  situation  his  predictions  were  within 
x  (p  g/m^)  y  percent  of  the  time  for  that  particular  data  category  using  the 
verification  data  set.  It  must  be  recognized  that  the  percentage  of  error  is 
not  displayed  in  the  above  technique;  i.e.,  the  iminbet  difference  between  a 
predicted  value  of  2  and  measured  value  of  4  is  the  same  as  between  that  of 
100  and  102  but  the  percent  error  is  quite  different. 


d.  Percent  Error  Distribution 

In  order  to  evaluate  the  percentage  error  and  make  the  verification 
results  usable  for  all  ranges  of  concentrations,  the  following  procedure 
should  be  followed: 

Using  the  same  categories  previously  discussed,  prepare  frequency 
diagrams  of  the  percent  error,  E ; ,  as  shown  in  Fig.  E-4: 

X  -  (X  -  X  ) 

K ;  =  — — = — ^— = — —  *  100, 

(X  -  X  ) 
o  b 

wh  e  re : 

E;  =  percent  error  per  case, 

=  model  calculation  without  background, 

X  =  observed  air  quality,  and 

Xb  =  estimate  of  background  air  quality. 

The  display  of  Fig.  K— 4  will  also  enable  the  model  evaluator  to  discuss  model 
bias,  F. ,  randomness,  o  y*-,  and  overall  variance,  o  j- ,  tor  each  category: 


who  re : 


1  n 

E  -  —  V  E  .  (bias), 

n  i 


i  -  -  2 

,,  •  --  y  (K  -  !•)“'  represents  the  randomness,  and 

u  E  n  i 
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Th  i  ,s  i  n  f  o  rmat  ion  .mild  also  he  displayed  on  a  percent 
frequency  diagram  similar  to  Fig.  E-5.  In  this  diapiam 
and  position  from  50  percent  frequency  of  the  line  would 
normality,  and  bias,  respectively. 


error  --  cutr.ul  at  ive 
the  slope,  curvature, 
indicate  randomness, 


In  Fig.  E-5  a  horizontal  line  through  zero  percent  error  represents  a 
perfect  model.  A  very  good  model  with  no  bias,  allowing  for  ihe  random  nature 
of  the  model  input  parameters  would  possibly  look  likP  the  distribution 
represented  bv  the  dashed  line.  An  example  application  of  this  diagram  would 
be  as  follows:  suppose  the  modeler  was  trying  to  substantiate  a  specific 

model  calculation  of  8.  He  could  say  that  in  the  category  represented  by 
Fig.  E-5,  solid  tine,  95  percent  of  the  time  bis  model  has  an  error  of  20 
percent  or  less.  Therefore,  only  in  one  out  of  twenty  cases  when  8  is  pre¬ 
dicted  (by  the  model)  will  the  actual  value  Xo  be  10  (i.e.  —  2  =  C  ii—  Xo )  /  X0  J  or 
greater,  or  there  is  a  9  5’<  confidence  in  the  value  being  10  or  less. 


e.  Emission  Sub-model  Adjustment 

Because  each  of  the  inputs  to  the  model  have  uncertainties  expressed 

as  a  standard  deviation  (  o),  the  model  output  cannot  be  expected  t.o  match 

perfectly  with  the  measured  concentrations.  However,  if  the  actual  variations 
of  the  input  parameters  are  independent  of  each  other,  and  the  input  para¬ 

meters  are  the  true  mean  values  of  their  distributions,  then  the  percentage 
error  assuming  Xb  =  0,  E  =  (x  c  model  -  x0  observation),  should  be  normally 
distributed  with  mean  zero  and  a  total  standard  deviation,  aT  .  If  the  mean 
E  is  significantly  different  from  zero,  as  determined  by  a  statistical  test, 
then  the  question,  What  is  causing  the  error?  should  he  asked.  Assuming  that 
adequate  quality  assurance  is  maintained  for  the  measured  pollution  concentra¬ 
tion  and  meteorological  data,  the  largest  uncertainties  will  probably  be 
either  in  the  emission  inventory  or  within  the  model  itself. 

Although  in  theory,  the  measured  product  ot  the  concentration,  x0 , 

observation  and  the  wind  velocity,  u,  is  a  flux  (mass/ area-t  ime )  and  can  be 
used  with  a  model  to  compute  the  emissions  Q  (gm/sec),  this  inverse  of  the 
standard  modeling  technique  is  not  practical.  A  possible  technique  to  test 
the  emission  inventory  would  be  to  break  t  lie  data  sets  into  two  s  ign  i  f  icant ly 
different  meteorological  conditions  such  as  stable  and  unstable.  If  the 
emission  inventory  is  in  error  and  the  model  is  consistent,  E  should  show 
similar  bias  from  zero  for  both  periods. 

If  the  K  values  are  significantly  different  lor  both  conditions  the 
emission  inventory  and  the  model  mav  both  be  in  error. 

Such  an  analysis  can  be  used  to  justify  modifying  the  emission  inventory 
Q  to  set  E  to  zero  for  both  data  sets  as  a  whole.  [f  this  adjustment  of 
Q  is  valid,  then  a  th^rd  'lata  set  ,  independent  of  the  two  used  as  above,  an 
hi'  used  to  determine  K.  If  E  is  not  significantly  different  from  zero  for 
this  new  data  set,  then  the  adjusted  Q  values  would  appear  to  be  the  valid 
ones  to  use  for  further  study  of  {he  model. 
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Fig.  E-5.  Percent  Error  -  Cumulative  Frequency  Diagram 


I  .  Model  Verification  I)  i  hi.  u  ■,  s  i  on 


A  subjective  discussion  should  accompany  each  com  par  i  son ,  pointing  out 
various  aspects  specific  to  the  modi  1  and  data  base  being  used.  Foe  discus¬ 
sion  wi l 1  depend  upon  the  interpretation  of  the  investigator.  Next,  the  data 
should  be  mult  istrat  i  tied  using  the  above  single  c  jtogut  i  ;  t  instance,  one 
stratum  may  be  all  cases  during  winter  morning  periods  with  wind  flow  from 
a  northerly  direction,  atmospheric  conditions  stable,  and  low  wind  speeds.  In 
this  case  the  categories  would  be  mutually  exclusive.  Ibis  data  comparison 
should  also  be  presented  in  both  frequency  distribution  and  scatter  diagrams 
of  observed  versus  predicted  if  there  are  more  than  30  examples;  and  scatter 
diagrams,  only,  if  there  are  less  than  30  examples.  Eat h  distribution  and/or 
scatter  diagram  should  be  evaluated  and  subjectively  discussed. 

It  is  recogn'zed  that  there  is  a  great  expense  involved  in  extensive 
model  calculations  to  fulfill  the  above  requirements.  Perhaps,  in  practice, 
models  can  be  evaluated  generally  in  the  broad  single  category  groups  using  a 
subset  of  the  total  data  set  to  evaluate  cumulative  frequency  distributions. 

A  few  of  the  mul  t  i  s t r at i f ied  categories  should  then  be  selected  for  specific 
model  performance,  preferably  those  with  the  most  data.  When  the  model 
is  then  used  in  a  specific  application  and  that  application  is  in  question  and 
not  previously  verified,  the  data  base  would  be  reevaluated  for  the  specific 
situation  in  question. 

3.  CONCLUSIONS 

Upon  completion  of  the  above  analysis,  assuming  a  good  air  qualitv  data 
base,  the  modeler  would  have  a  good  idea  of  how  valid  his  model  really  is.  In 
addition,  an  independent  model  user  or  Environmental  Impact  Statement  (EIS) 
reviewer,  would  be  able  to  scan  the  verification  report  and  obtain  a  quick 
feel  of  how  the  model  performs  compared  to  a  real  data  set. 

The  recommended  verification  procedure  discussed  above  does  not  require 
s  pc  c  i  tic  mo  <1  el  correlation  coefficients  or  any  other  number  on  which  to 
dele  nr.  i  ne  a  pass/ fail  grade  for  the  model. 

Trie  tinal  decision,  therefore,  as  to  whether  the  model  is  valid,  will  not 
depend  upon  any  specific  performance  limitation;  it  will  depend  upon  the 
user's  ability  to  support  the  use  of  the  model  in  any  instance,  based  on 
the  verification  analysis. 
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APPEND  I /.  I- 


SI  Ai  I.STK.AE  ANALYSIS  oh  AIK  PDl.i.ll'I  ION  Ij/VIA 
ANN  TIIK  VALIDATION  Of  THE  AQAM 


1.  INTRODUCTION 

This  section  first  discusses  several  past  airport  air  pollution  programs 
and  previous  attempts  to  validate  air  pollution  models.  Statistical  methods 
for  testing  several  hypotheses  related  to  both  the  analysis  of  air  pollution 
measurements  and  the  problem  of  model  validation  are  proposed.  Included  in 
the  discussion  are  descriptive  statistics  and  distributional  cons ideral sons , 
time  series  analysis,  and  multiple  linear  regression  analysis,  all  of  v,nich 
are  proposed  in  the  diagnostic  invest  igat  ion  and  analysis  jf  observe.’  and 
predicted  air  pollution  concentrations  for  use  in  defining  Lhe  limits  ot 
accuracy  of  the  AQAM .  Kmphasis  has  been  placed  on  the  details  of  general 
interest  and  justification  of  these  techniques  rather  than  on  the  presentation 
of  the  rigorous  mathematical  formalism  behind  each  one.  The  essential  math¬ 
ematical  details  and  computational  methods  employed  in  these  techniques  are 
both  lengthy  and  quite  complicated,  and  should  be  described  subsequent  to 
their  imp  1 emen t a t i on  in  another  report  by  Argonne  National  Laboratory.  for 
the  present.,  however,  the  reader  should  consult  references  to  the  literature 
that  are  provided  for  specific  mathematical  details  of  the  techniques  de¬ 
scribed  here. 


2.  PAST  AIRPORT  AIR  POLLUTION  AND  MODEL  VALIDATION  PROGRAMS 

Past  airport  air  quality  monitoring  studies  including  Thayer, ^ ^ 
Platt  Rot  e ,  and  MacWalters  et  a  1 .  have  failed  to  observe  concen¬ 

trations  of  air  pollutants  significantly  above  local  background.  In  order  to 
understand  the  complexity  of  validating  or  verifying  or  delining  the  limits  ui 
accuracy  of  airport  air  pollution  models,  one  must  appreciate  the  difficulties 
involved  in  monitoring  ambient  air  quality  levels  at  airports  and  in  attempt¬ 
ing  to  relate  them  to  aircraft  activity.  An  investigation  of  recent  airport 
air  pollution  monitoring  programs  was  therefore  performed. 

In  the  final  report  to  the  EPA-NERC,*  "Air  Quality  Measured  at  Atlanta 
International  Airport  before  anil  during  Experimental  Aircraft  Taxiing  Opera¬ 
tions,'’  MacWalters  et  a  1  . ^ have  conc  luded  that  : 

it  is  not  possible  to  detect  statistically  significant 
reduction  in  air  quality  levels  which  could  be  attributed 
to  the  emissions  control  procedures...  and  the  analysis  of 
air  quality  measurements  failed  to  show  a  statistically 
significant  improvement  as  a  result  of  lhe  application 
of  aircraft  emissions  control  procedures...  In  view  of 
the  limited  application  of  emissions  controls  and  the  result¬ 
ant  small  change  in  emissions  during  any  given  hour,  a 
substantially  longer  sampling  period  would  be  required  to 
detect  a  st.it  ist  ii  .illy  significant  >  t  t  ec  t  . 
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Plume  rise  effect  and  masking  of  aircraft  aerolite  trie  signals  t>y  back¬ 
ground  were  advanced  as  possible  explanat ions  for  the  results  of  the  Atlanta 
study. 

One  of  the  earliest  attempts  to  study  the  sensitivity  and  accuracy 
of  an  air  quality  computer  model  for  airports  was  performed  by  Platt  et  al.*--^ 
on  a  model  called  "A1REC."  Meteorological  parameters  were  subjected  to  a 
sensitivity  analysis  using  the  average  and  maximum  value  (the  maximum  obtained 
as  the  average  of  the  three  highest  values)  as  the  two  levels  tor  each  tactor 
being  considered.  The  resulting  changes  in  the  calculated  concentrations  were 
then  evaluated  as  a  percentage  of  the  data  values: 

The  results  of  the  sensitivity  analysis  for  the  meteorologi¬ 
cal  parameters  have  a  direct  bearing  on  the  selection  of 
short-term  periods  used  to  calculate  maxima  concentrations  in 
the  impact  analysis.  Hence,  time  periods  that  exhibit  low 
wind  speeds,  stable  atmospheric  conditions,  and  high  wind 
persistence  in  the  prevailing  direction  should  be  selected. 

The  accuracy  of  the  model  was  investigated  indirectly  in  the  analysis 
of  various  control  techniques.  Source  parameters  including  strength  and 
location  of  emissions  sources  with  respect  to  a  receptor  were  studied.  These 
i nves t iga t ions  reinforced  the  following  limitations  to  the  validity  of  the 
mode  1 : 

(1)  the  model  is  limited  to  time  periods  that  are  much 
longer  than  the  characteristic  times  of  individual 
aircraft  activity;  and 

(2)  the  model  is  limited  to  receptors  that  are  not  in  prox¬ 
imity  to  the  emission  sources  (closer  than  100  meters). 

Such  limitations  of  spatial  and  temporal  resolution  of  any  air  quality 
computer  model  need  to  be  established  in  defining  the  accuracy  of  the  model. 

Under  EPA  contract,  GEOMET  has  completed  a  program  of  air  quality 
measurement,  model  development,  and  model  validation  at  Washington  National 
A  i rpor  t  . 

Statistical  analysis  techniques  for  model  validation  of 
modi  f  ied  NREC  models^-^  are  developed  in  the  Geomet  Task  1 
Model  Verification  --  Aircraft  Emissions  Impact  on  Air  Qual  it'j. 
tion  of  the  model  involves: 

(1)  Comparison  of  versions  of  the  NREC  model  —  performed  by 
comparing  frequency  distributions  for  observations  and 
predictions  for  specified  pollutant-location-averaging 
time  combinations.  Other  methods  for  analyzing  the  error 
of  the  predictions,  such  as  the  Wile  oxen  test  or  the 
t-test  for  paired  comparisons  are  suggested,  based  on  an 
examination  of  the  data. 

(2)  In  comparing  predicted  and  observed  annual  means,  the 
bias  of  predictions  are  proposed  to  be  tested  -isnv  the 
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b  i  ii  i  >m  i  .‘i  I  test  to  determine  i  i  t  li  e  I  e  in  a  significant 
deviation  from  the  expected  value. 

O)  Comparison  of  the  number  <>|  times  t  tie  predicted  and 
observed  values  exceed  air  quality  standards.  This 
required  fitting  lognormal  distributions  to  the  avail¬ 
able  values  of  predicted  and  observed  data,  and  com¬ 
paring  these  to  the  standard  in  order  to  calculate  the 

number  of  times  during  the  year  that  the  standard  is 
exceeded  for  each  distribution. 

(4)  Determination  of  the  conditions  under  which  the  NRfX 
model  operated  best  and  worst  bv  stratification  of 
the  data  into  subsets  of  specific  me t eorol og : ca 1  and/or 
emissions  conditions  and  using  the  Ko  1  umogorov-Sir.i  ruov 
test  (referenced  below!  to  compare  frequency  d i str ibut  ions 
of  observed  and  predicted  concentrations. 

In  the  NREC  and  the  modified  NKhl.  (GEOMI-.T) 

mode  1  ... 

The  steady-state  Gaussian  plume  model,  which  is  the 

kernel  of  the  concentration  calculation  and  assumes 
steady  conditions  during  the  period  ol  calculation 
(in  this  case  one  hour),  is  c  h  a  x  ac  t  e  r  i  s  t  i  c  a  1  1  y 
not  expected  to  give  good  results  when  compared  to 

observations  on  a  paired-comparisons,  hour-bv-hour 

basis.  Such  a  model,  when  performing  well,  will 
reproduce  distributions  and  means  to  a  reasonable 
degree,  and  this  will  be  highly  useful  tor  most, 
purposes,  including  those  of  making  comparisons  with 
standards,  and  studying  the  impact  of  various  types  of 
contributing  sources  ... 

Several  measures  (descriptive  statistics)  and  statistical  tests,  in¬ 
cluding  the  t-test  for  paired  means  ( Grayb  i  1  1  ;  Itrown  1  ee  ,  simple  linear 
regression  (Kao),^~®  the  Spearman  rank  correlat  ion  coefficient  and  matched 
pair  sign  rank  test  (Kendall)  ,^"9  and  the  Ko  lomogorov-.Smi  rnov  goodness  -i 
fit  lest  (  Ko  lomogor  ov  )  ^  “  10~*  2  wi,  re  discussed.  These  techniques  wore  proposed 
for  the  statistical  treatment  ol  the  data,  but  are  not  related  to  the  testing 
of  specific  hypotheses  concerning,  the  model  valulat  ion.  The  discussion  o! 
e .ach  technique  includes  the  statement  of  the  null  hypothesis  and  abbreviated 
discussion  of  the  underlying  assumptions  and  ine thodo 1 ogy ,  along  with  the 
format  of  the  printout  for  each  test.  While  it  is  rei  'gi'.ized  that  smh 
techniques  are  the  subject  of  numerous  textbooks  on  statistics,  no  retereuces 
were  provided  for  further  examination  of  discussions  ot  these  methods. 

T1hj  appropriateness  ot  different  measures  and  tests  is  not  adequately 
developed  and  neither  are  the  inferences  made  subsequent  to  the  application  ol 
the  tests.  The  fact  that  paired  comparisons  were  "not  generally  bad,"  or  that 
the  t-tesls  were  "good"  tolls  the  reader  very  little.  The  p-v.iluos,  <„•  i  lov>  I 
o  I  statistical  i  on  f  i  ■!<  nc e  at  wli  1 1  h  the  reversal  of  the  null  hypothesis  lakes 
place,  were  not  even  men!  join'd.  The  assumpt  ions  ol  normal  it v  and  i ndependt  in  e 
underlying  tin1  applies! ion  of  correlation  and  regression  analysis  in  particu¬ 
lar  appeared  no  place  to  he  supported  by  preliminary  investigations. 


At  Washington  National.  Airport  the  contributions  of  aircraft  and  airport- 
related  sources  make  up  only  a  very  small  fraction  of  the  total  air  pollution 
burden  in  the  airport  vicinity.  Since  the  aircraft  and  airport  contributions, 
which  are  the  most  important  elements  in  the  model  study,  represent  such  small 
fractions  of  the  measured  concentrations  and  since  no  correlation  between  them 
was  demonstrated,  one  cannot  be  sure  that  the  model  as  it  relates  to  these 
factors  has  been  validated. 

Only  a  few  air  quality  computer  models  have  been  developed  specifically 
for  simulating  the  impact  of  the  airport  complex  on  ambient  air  quality 
(Thayer;F-l  Platt  et  al.;^-^  and  Rot e^“3 i 1 3 i ) .  a  recent  study,  A  Survey  of 
Computer  Models  for  Predicting  Air  Pollution  from  Airports ,  prepared  for  the 
FAA  by  Haber, F~15  concludes  that  no  properly  controlled  verification  of 
airport  models  has  been  performed  to  date: 

Airport  air  pollution  models,  recognizing  the  uncertainties 
inherent  in  the  transport  and  diffusion  process,  use  the 
Gaussian  type  description.  These  have,  however,  not  been 
properly  validated...  Adequate  validation  has  been  hampered 
by  the  difficulties  in  obtaining  reliable  data  and  by  the 
complex  nature  of  the  problem. 

Among  the  defects  in  airport  air  pollution  computer  programs,  Haber*'-*'’ 
notes  the  large  volumes  of  complex  input,  the  need  for  large  computers  lor 
the  comprehensive  air  quality  computer  programs,  the  lack  of  appropriate 
attention  to  the  relative  significance  of  different  pollutants,  and  the 
program's  failure  to  recognize  large  uncertainties  in  both  the  input  data  and 
the  models  implemented  in  the  programs.  He  also  concludes  that: 

Since  meteorological  variables  can  be  observed  but  not 
controlled  emissions  data  should  be  substantially  under  the 
control  of  the  investigator. 

By  conducting  a  controlled  experiment  at  a  relatively  remote,  high 
volume  military  airport  where  accurate  statistics  on  aircraft  type  mix  and 
activity  are  available,  Haber^"*^  notes  that: 

The  basic  Gaussian  transport  and  diffusion  equations  can  be 
validated  and  any  shortcomings  ident i f ied .. .quest  ions  regard¬ 
ing  the  accuracy  of  sets  of  values  for  diffusion  coefficients 
and  their  appropriateness  for  line  sources  (as  has  been 
suggested  by  ANL)  should  be  validated  empirically. 

The  problems  of  evaluation  and  validation  of  air  pollution  models  are 
discussed  by  Brier^-*^  in  Statistical  Questions  Relating  to  the  Validation 
jf  Air  Quality  Simulation  Models ,  a  report  to  EPA-RTP .*  (Haber*"* ^  has  de¬ 
scribed  model  validation  as  "reproduc ibly  generating  adequate  estimates.") 
Brier's  discussion  enlarges  as  well  upon  the  appropriateness  of  statistical 
and  analysis  techniques  such  as  the  computation  of  the  mean  square  error  (MSE) 
and  regression  analysis  (including  robust  linear  regression).  He  fails  to 
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merit  inn,  However,  that 
a  location  parameter, 
results  of  regression 
tie  fact  that  a  great 
n i que  . 


tin  MSK  is  a  dispersion  parameter  win  rras  the  median  is 
In  discussing  the  effect  of  outlier  points  on  the 
and  correlation  techniques,  in  graphically  illustrates 
deal  of  thought  must  go  info  tin  selection  of  a  tech- 


Briers'  general  guidelines  and  suggestions  for  a  validation  procedure 
include  a  sensitivity  analysis  to  indicate  any  internal  ineonsisf  <  no  id?  in  tie- 
model  and  to  understand  parameters  that  dominate  the  process,  followed  by  MSh 
and  regression  analyses  supported  by  appropriate  diagnostic  investigati  .ns  o! 
input-output  errors  and  of  the  assumptions  concerning  normality. 

As  Koopmans  (personal  comnun ic  at  ion )  (see  also  Ref.  K-36)  points  out  : 

The  use  of  robust  methods  to  lit  regression  linos  here  woiiic, 
conspire  to  make  the  model  look  much  better  Ilian  it  reallv 
is.  A  scatter  diagram  may  well  be  useful  for  displacing 
this  data,  but  fitting  a  regression  line,  witli  all  of  the 
corresponding  problems,  appears  a  very  poor  choice  ot  method, 
since  much  of  the  important  information  lies  in  the  outliers. 

Possibly,  a  combination  of  regression  lines  with  independent 
explanation  of  the  outliers  would  he  useful. 

3.  DESCRIPTIVE  STATISTICS  AND  1>  1  STK 1 BUT IONAL  CONS  RATIONS 

A  great  deal  of  attention  has  been  tocused  on  the  issue  of  determining 
the  distributional  characteristics  ot  ambient  air  pollution  concentrations. 
Most  researchers  will  agree  that  ait  pollution  data  appear  to  have  a  skewed 
underlying  distribution.  There  are  several  distributions  that  could  te  used 
for  expressing  this  type  of  skewed  data.  The  actual  or  empirical  distribution 
ot  air  pollution  concentrations  do  not  have  negative  values  and  ar--  usually 
asvmmettival  with  the  longer  tails  extending  above  the  median  concent i at  i  ■  - 1 :  . 
The  lognormal  distribution  lias  been  proposed  as  accurateiv  describin'..  ,i  i  r 
pollutant  concentrations  by  Zimmer  and  l.arscn,  I  a  r  s  en  ,  F  ~  1  .  1  .mb 

Mage .  During  the  EPA  Symposium  on  Statistical  Aspects  ot  Air  Oua  1  1  t v 

Data  (October  1974),*  Benarie,  Knox,  Pollack,  Turner,  and  1. viin  all  demon.,  t  rat  <  d 
the  usefulness  of  the  frequency  distribution  in  the  analysis  oj  air  pollutant 
concentrations.  Standard  statistical  calculations  mav  he  applied  il  we  first 
transform  the  concentrations  to  their  logarithms  base  in.  Hie  geometric  mean 
is  i  a  1  c  ti  I  a  t  <  •  d  as  the  a  u  t  i  1  o  g .  i  r  i  t  1 1  m  id  the  mean  logarithm.  the  standard 
geometric  deviation  is  the  an!  i  logarithm  of  tlu  si. unlaid  deviation  of  the 
logarithms.  It  tin-  data  follow  a  lognormal  distribution,  then  tin  eumulat  i v . - 
percentage's  fall  approximately  on  a  straight  line  on  I  or. -pi  obab  i  1  i  t  v  pa  pet  . 
The  b(l  percentile  point  is  given  by  the  geometric  mean  and  the  slope  of  the 
1 i ne  is  given  by  the  geometi ic  standard  deviation. 

Larsen*1-'  *  has  concluded  that  "concentrations  are  appr  ox  im.it  <■  1  y  lognor¬ 
mal  I  y  distributed  for  all  pollutants  in  all  cities  lor  all  averaging  times," 
based  on  actual  air  quality  data.  Larsen  contends  that  the  lognormal  disti  i~ 
but  ion  provides  a  simple  means  foi  predict  ion  purposes  especially  useful  in 
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air  pollution  control  applications.  Mage*"~2<)  by  treating  diffusion  as  a 
stochastic  process  has  obtained  the  censored  three-parameter  lognormal  model, 
of  which  Larsen's  standard  two-parameter  model  is  a  special  case.  The  new 
model  has  been  compared  with  actual  air  pollutant  concentration  data  and 
appears  superior  to  the  standard  lognormal  model,  especial  1 y  in  its  approxima¬ 
tion  for  the  upper  percentile  points  of  the  cumulative  distribution  t unci  ion. 
It  is  probable  that  this  improvement  results  not  just  from  addition  of  another 
parameter  but  from  the  fact  that  the  new  model  is  more  appropriate  ior  the 
underlying  physical  processes  that  generate  atmospheric  pollutant  concen¬ 
trations. 

The  method  for  estimating  the  parameters  of  the  censored  three-parameter 
lognormal  model  is  rather  complicated  in  comparison  to  that  for  the  two- 
parameter  model  and  involves  numerical  optimization  for  determining  a  best¬ 
fitting,  cumulative  distribution  function,  using  a  weighted  least  squares 
approach.  The  model  is  "censored",  meaning  that  it  will  predict  negative 
concentrations  that  are  simply  taken  to  he  zero.  The  censored  three-parameter 
model  proposed  by  Mage  appears  to  agree  or  f  it  the  actual  data  in  the  upper 
tail  of  the  distribution  more  satisfactorily  than  does  Larsen's  standard 
two-parameter  model.  Wh i 1 p  Larsen's  model^-21  may  be  adequate  for  making 
assumptions  about  control  strategies,  the  place  for  more  complicated  distribu¬ 
tions  like  the  one  proposed  by  Mage,^-^  £s  £n  fitting  large  bodies  of  data  so 
as  to  produce  a  better  estimate  of  the  behavior  of  the  upper  tail  ol  the 
d i s  t  r ibut ion . 

An  important  aspect  of  any  sampling  program  is  to  determine  the  fraction 
of  time  in  wh  ich  a  certain  concentration  is  exceeded.  The  empirical  cumula¬ 
tive  frequency  distribution  plot  can  be  used  to  get  the  percent  of  the 
sample  population  exceeding  any  specified  concentration.  However,  one  must 
keep  in  mind  the  accuracy  and  resolution  obtained  in  estimates  associated  with 
the  tail  of  the  distribution.  A  very  few  observations,  sometimes  relerred  to 
as  extreme  values,  are  usually  available  on  which  to  base  the  extreme  percen¬ 
tile  points.  Since  the  standards  promulgated  by  the  Knv  i  r  mmenla l  Protection 
Agency  are  stated  in  terms  of  exceedance,  the  frequency  ol  occurrence  of  high 
concentrations  are  of  major  importance.  For  this  reason,  the  straight  line  or 
more  complex  fitted  cumulative  frequency  equations  neon  to  be  critically 
examined  and  interpreted,  and  supported  by  diagnost  ic  invest igat tons  ot  the 
basic  assumptions  concerning  normality  ami  i  nde pendem  e .  Aitehison  and 
Brown^-22  provide  a  general  discussion  ol  tie  lognormal  theory  and  distribu¬ 
tion  methods,  and  should  be  consulted  for  general  discussions  of  lognormal 
d i s t  r ibu  t  ions  . 

Another  technique  known  as  the  boxplot  has  been  suggested  by  J.  W, 
rukey*"^!  for  examining  distributional  characteristics.  The  boxplot  technique 
consists  of  building  a  st em-and-1 eaf  data  structure,  similar  to  a  histogram, 
in  which  stems  correspond  to  higher  order  digits  and  leaves  to  lower  order 
digits  in  the  observations.  The  structure  is  used  to  compute  plotting 
positions  for  outlier  points  and  extreme  values,  and  a  box  indicating  the 
interquartile  range  and  the  median.  The  boxplot  produces  a  horizontal 
schematic  display  that  indicates  the  distributional  character  of  the  observa¬ 
tions,  L'n  r  tin-  assumption  o  1  no  rma  1  i  t  v  ,  a  certain  pe  r  i  n  t  ago  of  the 
observations  fall  mi' s  i  te  ,  >  t  f  tie  i  n  t  <■  r  u  i  .u  r  t  .  1  <■  range.  <  )u  t  1  i  i  s  and  ext  r  erne 
/al  i"'.  are  ri- 1  e  rr.  i  re  d  u  ■  a  result  "  t  t  1  ■  »  i  I  i  ng  ■  i  it  s  i  u  ■  <  •  1  t  h  •  i  u  l  •  r  q  ua  i  t  i  1  e 


range  ami  the  number  of  siu.li  data  values  exceeding  the  proportion  of  observa¬ 
tions  normally  found  in  the  tails  of  the  normal  d  '  s 1 1  i  bu t i on  .  The  application 
ol  the  boxplot  technique  to  the  original  data  and  to  ft."  log  10  transte  .1  01 
the  data  will  provide  insight  into  the  a :•  sump ;  ion  that  the  data  are  intrin¬ 
sically  norma l ,  that  is  to  say,  lognormal.  This  technique  otters  new  perspec¬ 
tives  on  the  distributional  character  of  the  data,  i:  at; lies  outliers,  and 
provides  estimates  of  several  descriptive  statistics. 

Sa  1 1  zina  n^-"^  lias  suggested  a  nomographic  representation  tor  estimating 
confidence  intervals  for  mean  values  and  determination  of  confidence  intervals 
for  proportions  of  a  population  exceeding  a  spec  if  led  emu  ent  rat  ion  i  1  th<j 
data  follow  a  normal  or  lognormal  distribution.  Johns,  <11  ag-i  Leone*"---5  provide 
a  discussion  of  confidence- i ntorval  estimation  of  the  por,el  at  lon-t  urn  j  1  a  t  i  ve- 
distribution  function.  Such  methods  may  prove  useful  wfi-  ri  comparing  two  or 
more  cumulative  frequency  distributions  for  significant  dilterenc.es,  and  also 
should  provide  insight  ns  to  the  accuracy  of  specific  percentile  points  that 
are  based  on  observed  or  predicted  data  being  compared  to  air  quality 
s  t  andards . 


Most  natural  phenomena  are  observed  as  statistical  distributions  and 
can  be  represented  by  a  frequency  curve,  the  parameters  of  which  can  be  esti¬ 
mated  from  the  observations.  As  R.  McCammon^-^^  has  suggested: 


Natural  events  are  the  product  of  compound  effects  resulting 
from  multiple  causes  and  the  observed  distributions  are  not 
always  representative  of  homogeneous  populat ions, but  rather  a 
mixture  of  two  or  more  subpopulations  each  of  wh ten  possesses 
unique  parameters.  An  observation  therefore  can  represent  a 
sample  from  one  of  several  popul  at  i  ons .  ...The  statistical 
problem  that  arises  is  the  dissection  ot  a  mixed  frequency 
distribution  into  its  individual  components... 


Background  concentrations  for  various  pollutants  superimposed  on  concen¬ 
trations  for  which  modeled  emissions  activity  are  known  to  occur  can  often 
cause  <  he  empirical  frequency  distribution  to  appear  to  he  bimodal.  This  i 
normally  reflected  in  the  cumulative  distribution  function's  depart  no 
from  a  straight  line  when  plotted  on  log  paper.  The  inability  to  relate 
pollution  levels  to  emissions  activity  can  usually  be  explained  by  the  masking 
effect  the  background  concentrations  have  on  the  pollution  signal.  A 

graph  ica 1 -ana  1 vt ic  technique  for  dissecting  normal  or  lognormal  distributions 
by  graphical  estimation  of  the  parameters  for  the  two  distributions  has 
been  proposed  bv  MoCammon  .  Dissecting  the  distributions  is  accomplished 

through  the  use  of  graphical  first  estimates  ol  the  lognormal  components 
followed  by  a  regression  analysis  that  provides  final  est  i mates .  A  cn i -square 
test  cart  he  used  to  lest  the  hypotheses  that  the  data  !■  present  a  mixture  of 
two  lognormal  distributions.  Such  an  approach  could  le-lp  to  identify  the 
distribut  ion  ot  hack  ground  concent  rat  ions  and  at  the  same  1  ime  would  prodm  <•  a 
better  straigtit  line  fit  for  t  tie  higher  percentiles  it  tie-  cumulat  ive  fre¬ 
quency  distribution  that  is  of  pa  r  t  i  <  u  1  ar  interest  i  r  1  d<  - 1  1  n  1  ng  the  limits  of 

.icnir.'ti  v  of  air  pollution  predi<  t  ions.  Having,  ideiil  if  led  a  background  g  1  vi; 

using  such  a  technique,  descriptive-  statistics  fOt  background  subtracted 

aerometric  signals  can  be  .imputed  easily. 


Robust  techniques  for  the  estimation  of  the  mean,  median,  and  confi¬ 
dence-interval  estimates  have  been  explored  and  developed  by  Dixon  and 
Tukey,^-27  Tukey  and  Me  Laugh  1 i n , ^ -2&  and  Gross,  ^"29  The  computation  of 
descriptive  statistics  are  often  made  under  the  assumption  that  tiie  data 
follow  a  normal  or  intrinsically  normal  distribution.  In  practice  this 
assumption  is  often  false.  A  procedure  called  the  wave- i nterva 1  ,  proposed  by 
Gross*1-^  lor  the  construction  of  a  confidence  interval  that  is  not  dependent 
on  the  normality  assumpt ion, has  been  found  to  perform  very  well  and  is  recom¬ 
mended  for  general  use  in  place  ol  the  t-interval  for  sample  sizes  of  eight  or 
more.  It  is  assumed  that  the  distribution  of  the  data  is  unimodal  and  sym¬ 
metric  and  is  most  likely  long-tailed.  This  technique  is  recommended  for  the 
estimation  of  descriptive  statistics  such  as  the  mean,  median,  and  confidence 
limits  for  small  sample  sizes  as  a  means  of  compensating  for  extreme  points 
and  outliers.  Small  sample  sizes  will  undoubtedly  arise  in  the  stratification 
of  available  data  on  factors  such  as  level  of  emissions  activity,  time  of  day, 
etc  . 


It  is  important  to  emphasize  the  role  of  graphics  related  to  the  previous 
discussions.  Plots  of  actual  and  fitted  cumulative  frequency  distributions, 
boxplots,  and  other  descriptive  statistics  are  essential  aids  in  the  analysis 
and  interpretation  of  air  pollution  data.  Graphics,  as  we  shall  see  in  the 
next  section  are  also  extremely  important  in  time  series  analysis  and  in 
plotting  the  data  itself. 


4.  TIME  SERIES  METHODS 

This  section  describes  time  series  analysis  as  it  relates  to  the  statis¬ 
tical  analysis  of  air  pollution  and  meteorological  data.  The  time  dependent 
nature  of  such  data  and  the  liklihood  of  serial  correlation  found  in  such 
measurements  make  obvious  the  need  for  analyzing  such  data  in  terms  of 
its  frequency  composition.  This  section  provides  background  on  time  series 
concepts,  attempts  to  justify  the  use  of  time  series  analysis,  and  discusses 
specific  analysis  techniques,  providing  the  reader  with  references  tot  further 
reading . 

Cleve land^--^  has  proposed  the  use  of  the  inverse  autocorrelation  of  a 
time  series  in  the  analysis  of  air  pollution  data.  Two  different  methods  tor 
estimating  the  inverse  autocorrelation  arise  from  two  different  methods  of 
estimating  the  spectral  density  —  autoregressive  and  periodogram  smoothing. 
The  estimates  of  the  inverse  au t oc or r e 1  a t i on s  can  be  used  to  assist  in 
identifying  a  parsimonious,  moving  average,  autoregressive  model  tor  the 
series.  The  feasibility  of  the  autoregressive-moving  average  models  lias  been 
demonstrated  by  Lee.^'^l  The  use  of  time  series  in  the  analyst  of  air  pollu¬ 
tion  monitoring  data  has  also  been  suggested  by  Harlow  and  Si ngpurwa 1 1  a , ^ 

Sal  tzman,^-^  Marcus,  ^-34  an(j  xiao  and  Hamming.  ^-35 

The  value  of  a  time  series  at  any  given  time  is  called  the  amplitude 
component  at  that  time.  A  useful  measure  of  the  "activity"  of  the  time 
series  is  its  power.  The  amplitude  is  a  positive  number  measured  in  the 
amplitude  units  of  the  problem  at  hand.  The  phase  is  a  dimensionless  para¬ 
meter  that  measures  the  displacement  of  the  sinusoid  relative  to  the  given 
t ime  origin. 
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An  extrusive  treatment  of  the  spectral  a  ria  1  v  s  :  s  oi  t  ime  series  h  a  ■ 
been  given  by  Koopmans  .  ^  Among  his  views  on  times  s-ries  are: 

In  the  same  sense  that  various  colors  ot  light  are  composed 
of  a  blend  of  monoch romat ic  component  ■  and  nusical  tones  are 
formed  by  a  supe  r  impos  i  t  ion  of  pure  harmonic--,  time  series 
can  be  constructed  by  composing  a  number  of  harmonic  func¬ 
tions  with  varying  frequencies,  amplitudes  and  phases.  livery 
time  series  has  an  explicit  spectral  representation  in  terms 
of  frequence,  amplitude  and  phase. 

Time  series  data  analysis  and  heorv  also  is  discussed  extensively  dv 
Br  i  1  1  inger  and  the  Fourier  analysis  oi  time  series  bv  B1  ootnf  it-  ld^ 

provides  considerable  detail  on  spectrum  analysis  and  discussion  of  the 
general  framework  ol  time  series  analysis  methods.  Additionally,  a  bibli¬ 
ography  of  time  series  analysis  papers  has  been  given  by  Wold.^~39 

The  principal  goal  ol  spectrum  analysis  is  to  decompose  the  power  of 
the  given  time  series  int'>  its  harmonic  components,  Lha.  is,  to  estimate  the 
power  spectrum  from  the  available  data.  The  estimated  spectrum  can  then  be 
used  to  gain  information  about  the  mechanism  that  generated  the  data.  The 
statistical  theory  of  spectral  analysis  has  been  based  on  certain  hypotheses, 
that  the  underlying  process  is  stationary  and  Gaussian  and  that  the  process 
m< an  is  zero  and  the  spectrum  is  continuous.  Some  ot  the  preprocessing 
operations  performed  on  the  data  before  a  spectrum  analysis  is  performed  ar> 
intended  to  bring  the  data  into  reasonable  conformity  with  these  hypotheses. 
A  non-zero  mean  and  possible  trend  in  the  mean  can  be  either  an  indication  of 
non-s tat  ionar i ty  or  of  a  drift  in  the  instrumentation  systems .  The  estimation 
and  removal  of  linear  or  quadratic  trends  that  best  fit  the  data  in  tne  least 
squares  sense  can  be  subtracted  trom  the  original  data  with  the  resulting 
residual  conforming  more  closely  to  the  stationary  hypothesis. 

Several  import  ant  reasons  for  t h e  use  of  time  set  its  analysis  are 
■  nit  1  i  noil  he  1  nw  : 


As  a  diagnostic  and  exploratory  tool  in  checking,  the 
assumpl ions  ol  independence  and  normality  needed  lot 
ot  he  i  statist  ical  tests,  particularly  those  related  to 
significant  testing,,  and  frequency  distribution  construe- 
t  i  on  t  och  n  i  q  tie  s  ; 

As  a  means  of  analyzing  the  serial  correlations  in  u  .e 
data; 

As  a  method  of  examining  frequency  dependent  cor i e 1  at  i ons 


and  i  nva  I  i  - 


(called  coherence  between  two  limes  series);  correlation 
methods  in  the  time  domain  —  for  example,  examining 
i  nni'l  it  ion  by  using  straight  line  methods  can  !  •  verv 
sen-  it  ive  to  the  assumption  of  independence  and  invali¬ 
dated  by  (tie  presence  o|  serial  correlation; 

As  an  '■  >p  I  or.it  "tv  tool  in  searching  for  per  iodic  i  t  v 
and  gross  vuii.it  toil  ot  power  with  frequency  to  provide 
insight  into  the  mechanism  behind  the  data,  and 


(5)  That  it  is  frequently  possible  to  describe  complex 
physical  mechanisms  better  in  the  frequency  domain  than  in 
the  time  domain. 

By  examining  periodicities  in  the  data,  we  have  another  way  of  looking 
at  some  of  the  fundamental  mechanisms  that  produce  diurnal  patterns  and 
pollution  episodes.  The  correspondence  between  aircraft  emissions  and  pollu¬ 
tion  concent r at  ions  should  be  readily  observable  through  an  examination  <>i  the 
spectral  parameters. 

Two  models  of  spectral  analysis  are  now  described  --  the  univariate 
spectral  analysis  and  the  bivariate  spectral  analysis  of  time  series.  Multi¬ 
variate  spectral  analysis  has  not  been  considered  for  the  present. 

Univariate  spectral  analysis  involve:  the  computation  of  the  power 
spectrum  or,  as  it  is  also  called,  the  estimated  spectral  density.  The  power 
spectrum  is  the  Fourier  transform  of  the  autocorrelation  function.  The 
Fourier  transform  of  a  sequence  of  autocorrelations  contains  precisely  the 
same  information  found  in  the  original  autocorrelations  themselves. 

Many  statistical  analysis  methods  that  operate  on  the  data  in  the  time 
domain  make  no  use  of  the  time  dependency  in  the  data.  For  example,  in  the 
estimation  of  the  frequency  or  cumulative  frequency  distributions,  the  data 
are  sorted  into  non-decreasing  order,  and  information  in  the  time  dependency 
of  the  data  is  no  longer  available.  An  important  property  of  time  series 
analysis  is  that  much  of  the  time  dependence  in  the  data  is  preserve.!.  The 
accumulated  spectrum  is  essentially  the  first  integral  of  the  spectrum  and  can 
be  used  to  indicate  a  range  of  frequencies  that  contain  a  substantial  propor¬ 
tion  of  the  total  power  of  the  time  series. 

Univariate  time  series  analysis  is  also  an  important  tool  in  the  design 
of  digital  filters.  By  examining  the  spectrum, the  effect  of  filtering  can  be 
seen  in  terms  of  the  frequency  composition  of  the  time  series.  In  this 
fashion,  digital  filters  can  be  developed  that  remove  frequency  components 
that  might  be  regarded  as  noise  prior  to  the  analysis  of  trends.  For  example, 
tf  the  original  data  are  sampled  every  five  minutes,  but  only  hourly  averages 
are  required,  then  the  frequency  components  in  the  data  due  to  noise  in  the 
five-minute  data  can  be  removed  or  filtered  out,  resn  1 1  inp  in  a  smoothed 
five-minute  series  that  can  be  used  to  calculate  the  one-hour  data.  In  a 
similar  fashion,  the  effect  of  the  intake-manifold,  averaging  bottles  can  be 
investigated  by  examining  the  spectrum  of  already  smoothed  pollution  concen¬ 
trations.  Similarities  between  the  spectra  of  two  series,  such  as  peaks  at 
similar  frequencies,  may  raise  the  possibility  that  the  series  are  related. 
To  investigate  such  possibilities,  we  compute  estimates  of  the  cross  spectrum 
of  the  two  series.  This  is  an  extension  of  the  definition  of  the  spectrum  and 
is  usually  estimated  by  smoothing  the  cross-per iodogram ,  and  is  performed  in 
the  bivariate  spectral  analysis  of  the  two  series. 

Bivariate  spectral  analysis  for  the  spectral  decomposition  01  two 
series  into  various  spectral  parameters  is  equivalent  to  descriptive  statis¬ 
tics.  These  include  the  estimated  spectral  density  of  each  series,  the  phase 
function  that  indicates  the  relative  shift  of  harmonic  components  of  the  two 
time  series  at  the  same  frequency,  the  transfer  function  that  indicates  the 
linear  relationship  between  th>-  two  series  in  the  same  wav  regression 


cue  f *  i  c  i  <-nt  s  relate  the  linear  i  e  1  at  i  onsli  l  p  between  two  data  sets,  and  the 
coherence  fnnct  ion  that  provides  an  est  imate  oi  the  f  r  t-ci  uenc  y-d  epe  nd  *  -n  t 
correlation  between  the  tw-'  series. 


Because  of  the  vast  amounts  of  numerical  information  contained  in  the 
estimates  of  spectral  parameters,  graphs  should  be  obtained  tor  the  spectral 
parameters  of  interest.  Vor  the  univariate  spectral  analysis  these  i  tic  luce 
the  estimated  spectrum  and  the  accumulated  spectrum.  In  the  :>  i  va  r  i  a  t  spec¬ 
tral  ana  lysis  are  included  the  spectral  densities  tor  ea- h  serifs,  as  well  a- 
the  phase,  transfer  and  coherence  functions.  Plot-  of  the1  origin.:!  and 
corrected  data  traces  should  also  be  obtained. 

A  general  outline  of  the  application  of  time  series  is  given  below: 

(1)  Collect  all  data  tor  times  series  analysis,  making 
efforts  to  ensure  that  the  data  are  feasible  and  if 
necessarv  adjust  the  data  for  comparability. 

(2)  Obtain  a  graph  of  the  t  ime  series. 

(3)  Perform  trend  analysis. 

(4)  Accommodate  and  adjust  for  seasonal  variations  if 
r equ i red . 

(3)  Adjust,  the  data  for  trend. 

(b)  Compute  and  graph  all  spectral  parameters  of  interest. 

Tt  is  proposed  that  spectral  analysis  be  applies  to  the  air  pollution 
and  meteorological  data  traces  for  each  station.  Missing  data  can  be  <.  or.:  «-••.- 
satfcd  for  by  simple  linear  interpolation,  which  does  nut  affect  the  estimated 
spectral  parameters.  The  purposes  of  these  time  series  analysis  techniques 
have  been  discussed  in  the  section  of  this  report  dealing  with  statements  of 
hypotheses . 


3.  "I'l/ril’l.h  I.  INK  AH  KKGKKKS  I  ON  ANALYSIS 

There  are  i  wo  main  uses  nl  (In-  multiple  linear  repression  technique. 
Tin’  first  provides  an  analysis  t  ••rhn  i  que  for  the  evaluation  ol  the  relat  ive 
effects  and  .signif  icance  ol  I  at  tors  at  led  ing  the  dependent  vaiiable.  In  this 
case,  the  dependent  variable  is  the  response  under  study  and  the  independent 
variables  are  the  actual  values  ol  parameters  related  to  the  responses.  The 
second  usage  piovides  an  important  framework  for  planning,  ami  analysis  of 
sensitivity  experiments.  This  technique  is  especially  useful  in  situations  in 
which  each  lespopn  or  dependent  variable  may  be  expensive  or  difficult  to 
obi  a i n  . 

Consider  t  lie  regression  <  >r)  ■  i.r  t  ion  of  the  form: 

Y  ts(1  +  B  |  X  |  t  B  j  X  ■]  +  ....  +  I: 
wb  e  re  : 

Y  -  the  dependent  variable,  which  is  the  response  quant  itv 

under  '  oils  i  do  r  a  t  ion. 
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=  independent  variables, 

6 j  =  regression  coefficients  (sometimes  railed  partial 
coefficients),  and 

E  =  the  random  error  term, 

A  discussion  of  the  assumptions  underlying  the  application  of  regression 
methods  can  be  found  in  most  statistics  texts,  such  as  Fitting  rq  ua  ti  ones  to 
! at' 1  by  Daniel  and  Wood,F-40  which  discuss  the  least  squares  method.  Those 
assumptions  concern  a  knowledge  of  the  correct  form  of  the  regression  equa- 
t  ion,  use  of  typical  data  in  the  analysis,  and  assurance  that  the  dependent 
variables  are  statistically  uncorrelated  (or,  more  generally,  are  said  to  be 
independent).  Less  important  assumptions  are:  that  the  dependent  variables 

have  the  same  (though  unknown)  variance,  that  the  values  of  the  independent 
variables  are  known  without  error,  and  that  the  random  error  is  normally 
distributed. 

The  use  of  regression  methods  and  exper iment a  1  designs  in  the  context 
of  performing  sensitivity  analysis  is  to  estimate  the  effects  of  the  independ¬ 
ent  variables  of  interest  (and  possibly  their  interactions)  in  a  reasonable 
number  of  experimental  runs  or  values  of  the  dependent  variable.  Fractional 
factorial  designs,  discussed  extensively  by  Cochran  and  Cox^-^  can  be  used 
for  this  purpose.  These  designs  provide  estimates  of  the  main  effects  and,  if 
desired,  some  low-order  interactions  with  a  minimal  number  of  computer  runs  or 
values  of  the  independent  variable. 

The  model  using  two  levels  of  the  factors  is  called  a  first-order  model 
and,  hence,  contains  no  higher  order  terms.  If  for  some  factor  a  quadratic 
component  is  to  be  estimated,  the  factorial  designs  must  be  renlanned  to  take 
that  factor  at  three  or  more  levels.  In  practice,  however,  the  quadratic 
effects  are  often  neglected  in  the  earlier  stages  of  experimentation  in  order 
to  make  decisions  about  future  experiments  on  the  basis  of  the  first-order 
mode  1  .  Orthogonal  polynomials  are  suggested  when  quadratic  effects  are  to  be 
e s t ima  t ed  . 


Tlie  fitted  equation  could  fail  to  be  statistically  signiticant  because 
1 1)  the  assumed  form  of  the  equation  is  not  the  true  or  (2)  the  variation  of 
observed  points,  though  random,  is  so  large  that  the  fitted  equation  could 
have  arisen  by  random  sampling  from  a  population  with  all  partial  regression 
coefficients  equal  to  zero. 

The  hypothesis  that  all  true  partial  regression  coe f f ic vents , 1  s  ,  equal 
:ero,  can  be  tested  by  an  F-test  of  the  variance  accounted  for  by  regression 
relative  to  the  error  variance.  The  F-value  calculated  from  the  data  can  be 
compared  with  the  tabulated  F-value,  ancj  the  hypothesis  that  all  true 
partial  regression  coefficients  equal  zero  can  be  rejected  if  the  calculated 
F-value  is  larger. This  same  hypothesis,  if  rejected,  also  allows  testing 
to  determine  whether  the  regression  accounts  for  a  significant  amount  of  the 
variation  observed  in  Y,  the  response  quantity. 

It  is  recommended  that  the  modified  Gram-Schmidt  analysis  of  variance 
he  used  to  perform  l nc  least  squares  solution  in  the  fact oi  izat i on  process  in 
solving  the  regression  equation.  In  the  modi  lied  Gram-Sc hm i d t  procedure,  all 
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columns  of  tho  cl  us  i  >  n  m  <1 1  :  ix  .*  r  <  or  I  linjciri.j  I  i  z<-d  nc  c  >■  sc  ;  1  v  will,  respect  to 
the  previous  or  tno/on.i  1  cot  umn  i>  the  overall  process.  'Jne  advantage  xn\  n*-<j 
is  in  dealing  with  i  1  l -to  ml  i  l  i  oned  hi:-  n  n-s  .  For  such  cases,  lire  mod  if  it i 
Gram-Schmidt  method  works  adeciuat  1  v  .  1 ~u  4  ’  '  ■' 

b.  RKC0MMENDAT1 ONS 

It  is  re comm ended  that  the  one -hour  aircraft  act ivit  v  tape  be  reduced 
from  the  discrete  flight  time  'lata  and  that  explanations  of  assumptions 
underlying  the  estimation  of  such  quantities  as  the  number  of  touch-and-go 
operations  per  hour  be  provided. 

Since  the  actual  aircraft  modes  discussed  here  will  very  likely  oe 
different  from  the  modes  to  he  used  in  the  AQAM ,  some  explanation  of  the 

rationale  behind  the  choice  of  AQAM  modes  should  be  developed. 

Since  a  vast  amount  of  numerical  information  will  be  involved  in  the 

observed  and  predicted  data  sets  and  will  also  be  generated  during  the  A'^AM 
model  validation,  it  is  recommended  that  computer  graphics  capabilities  be 

utilized  extensively  in  all  phases  of  the  investigation.  Graphic  routines  for 
general  two-dimensional  plotting  (scatter  diagrams,  lime  histories,  etc.), 
wind  and  pollution  rose,  effective  source  mapping,  and  two-d imens i ona  1  con¬ 
touring  will  be  required.  The  importance  of  graphics  capabilities  cannot 
be  overstressed.  Any  development  costs  will  be  substantially  offset  by  the 
contributions  these  techniques  will  make  in  the  analysis  process. 

ft  is  recommended  that  ail  statistical  analysis  of  data  from  the  Williams 
AFB  air  pollution  monitoring  program,  including  the  reduction  of  the  raw  data 
tapes  and  the  computation  of  cumulative  frequency  distributions,  be  coordi¬ 
nated  with  the  AQAM  model  validation  program  performed  by  Argonne  National 
Labor  at  or  v  . 

Time  series  analysis  is  recommended  tor  use  in  the  investigation  or 

assumptions  and  hypotheses,  for  the  analysis  of  serial  correlations  in  tic- 
data,  and  for  the  determination  of  frequency  dependent  correlations  (cohe¬ 
rence)  between  two  time  series. 

Multiple  regression  analysis  is  recommended  as  a  means  of  analyzing 
experimental  designs  in  the  determination  of  the  relative  significance  and 
importance  of  emissions  (including  at  least  total  on-base  sources,  aircraft, 
and  background)  and  meteorology  (including  data  stratified  on  wind  direction, 
wind  speed,  atmospheric  stability,  and  diurnal  patterns)  on  the  accuracy  ol 
AQAM  predictions.  The  following  considerations  are  recommended  in  the  overall 
detinit  ion  of  arcui.icy  of  the  AQAM: 

(I)  l(U'iu  i  t  i  i  at  i  on  ol  po  I  I  utant  (s  >  ,  time  period':  and  condi¬ 
tions  under  which  the  model  performs  best  and  wo i s t . 

(?)  Comparison  of  predicted  and  observed  exceedances  of  the 
air  quality  standards. 


F-l  . 

K- 2  . 

F-3  . 

F-4  . 

F-5  . 

F-6 . 

,■  -  7 . 

F-8  . 

F-9  . 

F-  10  . 

F-l  L  . 

F-L2  . 


!)  General  comparison  of  predictions  and  observations 
during  worst  case  air  pollution  episodes.  These  epi¬ 
sodes  should  include  time  periods  of  substantial  air¬ 
craft  activity  and  periods  of  persistent  a  Linos  pile  i  i  c 
stability  by  an  analysis  of  data  from  the  William;;  AKB 
air  pollution  monitoring  program. 
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Fig.  G-2e.  Cumulative  Frequency  Distributions  of  THC  at  Station 


OBSERVED  (THC-CH4) 


PROBABILITY 

Cumulative  Frequency  Distributions  of  NMHC  at  Station 


OBSERVED  (THC-CH4) 
AQAM  1  THC 
AQAM  2  THC 


PROBABILITY 
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Fig.  G-Aa.  Cumulative  Frequency  Distributions  of  NO  at  Station 
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Fig.  G-6b .  Frequency  Distribution  for  the  Upper  vlO%  of  Observed  NHMC  Concentrations  at  Station 

Williams  AFB  Hourly  n^ta 


Williams  AFB  Hourly  Data 


O'O  £ 

(J 


r 


SXN3A3  30  'ON 


185 


Williams  AFB  Hourly  Data 


Frequency  Distribution  for  the  Upper  ^10%  of  Observed  NMHC  Concentrations  at  Station 

Williams  AFB  Hourly  Data 


II/Observation 


Williams  AFB  Hourly  Data 


EVENTS 


Lop,  (AQAM  II /Observation) 

Frequency  Distribution  for  the  Upper  v j 0 Gf  Observed  NO  Concentrations  at  Station 


APPENDIX  H 


ANALYSIS  OF  HIGHER  REPETITION  RATE  DATA 


1  .  INTRODUCTION 


During  normal  operations  of  the  meteorological  and  aeromett  ic  monitor  i ng 
facility  at  Williams  AFR,  each  station  and  thus  each  experimentally  measure. I 
quantity  was  sampled  once  per  minute.  This  repet  it  ion  rati  was,  at  lirst  , 
considered  questionably  low,  especially  in  view  of  the  short  t  itno  scales 
associated  with  aircraft  operations  or  atmospheric  turbulence.  Thus,  the 
question  was  investigated  of  how  hourly  averages  ol  once  per  minute  compare 
with  hourly  averages  computed  on  the  basis  ol  higher  repet  it  ion  rate  data. 

The  data  sample  for  this  comparison  was  obtained  du.  1  ng  Mav  2  5  —  27,  1*17/, 

bv  operating  the  Monitor  Labs  9400  Data  Acquisition  System  at  a  higher  sam¬ 
pling  rate  at  the  expense  of  number  of  channels  sampled.  The  samples  chosen 
for  analysis  consist  of  17  hr  of  wind  speed  (WS),  wind  direction  (WD),  NO  and 
N0X  data  taken  at  6-sec  intervals,  and  3  hr  of  vert  ical  wind  speed  data  taken 
at  sampling  rates  as  high  as  once  per  2  sec. 

The  basic  analysis  approach  was  to: 

(1)  calculate  hourly  averages  using  rapid  scan  data; 

(2)  subdivide  this  rapid  scan  data  into  subsets  containing 
data  at  1-min  intervals; 

(3)  calculate  hourly  averages  using  the  once-a-minute 
data  in  these  subsets;  and 

(4)  compare  the  two  calculations. 


These  data  suggest  that  the  following  uncertainties  related  to  sampling 
.ate  mav  be  associated  with  the  once-pe r-mi nu t e  Williams  AIR  measurements. 


Pollutants: 


Meteorological  quantities: 


NOx 
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2X 

NO 

< 
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3%) 

WS 
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3  % 

WD 

8“'. 

(of  Og) 

°e 
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12% 

(usually  < 

8%1 

ov 

8% 

°w 

17% 

In  the  last  section  of  this  appendix,  the  "one-fifth"  power  law  is 
investigated.  This  law,  describing  the  time-dependent  increase  in  o„  with 
averaging  time  raised  to  the  power  0.2,  is  an  assumption  inherent  in  many  air 
quality  dispersion  algorithms,  including  AQAM.  These  data  suggest  that 
betvH-'r  averaging  times  of  3-60  min  this  jxiwer  is  0.27  +  0.12  in  agti.-im  ni 
wit:,  'he  one-fifth  hypothesis. 
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2.  lNHKI'KNUKNCK  UK  COMI'UTKI)  HOURLY  QUANT  IT  IKS  FROM  SAMPLING  FKKQIIKNGY 

From  the  rapid  scan  data  taken  on  days  145,  146,  and  147  of  the  experi¬ 

ment,  17  hours  of  windspeed  and  wind  direction  measurements  are  available. 

This  data,  taken  every  six  seconds,  can  he  used  to  compute  various  hourly 
quantities  such  as  hourly  average  windspeed,  hourly  o0  ,  etc.  If  all  the 
6-sec  data  values  for  1  hr  are  used  in  these  calculations,  the  sampling  fre¬ 
quency  if  the  data  is  once  per  6  sec,  whereas  the  sampling  time  ( a  )  is  1  hr  , 
since  hourly  quantities  are  computed. 

If,  however,  one  only  wanted  data  with  a  sampling  frequency  of  S=1 
minute-^,  one  could  simply  pick  out  every  tenth  data  point  from  the  S=(6 
sec)-!  data.  However,  because  there  are  10  of  these  data  points  every 
minute,  there  are  10  different  times  on  which  to  start  the  first  minute  of  the 
6=1  min-!  data.  So,  for  the  subsequent  analyses,  the  S=(b  sec)-!  data  is 
broken  up  into  10  independent  subsets  of  data  for  eacii  ol  the  17  hr.  Eacn 
subset  of  data  therefore  has  a  sampling  frequency  of  once  per  minute. 

L,e  t 

4)  i  =  the  data  points  that  make  up  the  complete  rapid  scan 
data  set,  with  sampling  frequency  S=(6  sec)-!.  i  is 
an  index. 

and 

=  the  data  points  that  make  up  the  j*-!1  subset  of  the 

S=1  min  data,  k  is  an  index  while  j  indicates  the  subset. 

Then  the  process  whereby  the  10  subsets  of  S=(]  min)-!  jata  are  ex_ 
traded  from  the  S=(6  sec)-!  data  can  be  represented  by: 

,fjl  -  4>  j 

*if  *j+10 

*jj  =  *1+20 

*  j  k  =  *  j  +  10(k-l ) . 


which  can  be  repeated  for  each  j,  j=l , 2 , . . . . 10 . 
the  data  has  a  sampling  frequency  of  1  min-!. 


The  asterisk  indicates  that 


The  object  of  this  part  of  the  analysis  is  to  compare  the  values  ol 
various  hourly  quantities  obtained  by  using  the  S=(6  sec)-!  data  (all  4 j  ) 
to  the  values  obtained  by  using  the  S= 1  min-!  data  (the  $*^). 
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Let  St(Xi)  be  some  statistical  process  that  operates  on  tin'  set  ol  data 
points  for  the  sampling  period  t.  For  example,  ST(Xj)  eon  Id  represent  an 
hourly  average  of  the  X£,  in  which  case: 


S  (X.) 
T  1 


s 


1  hr  (X 


1  N 

i)  ~  N  £  *i 


i=l 


where : 


N  =  the  number  of  X^  in  1  hr. 

Then  define  the  following  operations  on  the  meteorological  data,  and  let 
the  variables  indicated  represent  the  desired  hourly  quantities: 

®  =  Si  [lr  ( $  £ )  =  value  obtained  from  the  statistical 

process  ST  using  the  S=(6  sec)-^  data 
during  a  1-hr  sampling  interval  (t  =  1  lir). 

4^  =  Si  tlr  =  value  obtained  from  the  statistical  process 

ST  using  the  S=1  min-^  data  in  the  j^h 
subset  during  a  1-hr  sampling  interval. 

To  show  the  independence  of  Sj  ^ r ( X £ )  on  sampling  frequency  S,  one  must 
compare  the  value  $  (the  hourly  quantity  computed  from  S=(6  sec  )"*  data)  to 
the  4j  values  (the  hourly  quantity  computed  from  each  of  the  S= 1  min-^  data 
subsets).  But  because  there  are  ten  4*  values  (one  for  each  subset  of  v=  1 
min~l  data),  it  is  convenient  to  define  some  method  for  comparing  these  values 
without  having  to  look  at  each  4*.  However,  since  each  subset  j  k  h.  s 
the  same  chance  of  being  selected  as  a  unique  set  of  minute-averaged  data 
(assuming  a  random  start  time  for  data  collection),  all  4*  values  should  he 
considered  when  comparisons  with  4  are  made.  Thus  it  is  reasonable  to  dot  me 
a  "standard  deviation"  that  measures  the  distribution  of  the 
the  4  value  : 


4  j  va  lues  a roun d 


10 


To  I 

j  =  l 


(4*  -  4  ) 2 


1/2 


~U) 

wh e  re  : 

o  ( X  )  =  the 

A  1 ow  va 1 u< 
the  S= 1  min~l  d; 
therefore,  that 
o f  t  he  data  X  ;  . 

one  may  further  define  the  mean  value  <i|  the 
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( H-2  ) 


As  a  measure  of  how  well  the  value  #*  represents  the  individual  quantities  $j, 
use  the  standard  deviation: 


o(4*)  = 


10 


TZ>  I 

j-i 


(*. 


**y 


1/2 


(H-3) 


In  cases  where  S  defines  an  arithmetic  average,  $  will  be  exactly  equal 
to  ♦* .  For  other  processes,  however  4*  can  be  used  as  a  representative  value 
for  the  individual  fcj's  and  may  differ  from  the  value  t . 

As  an  example  of  this  approach,  suppose  that  ST(X^)  represents  the  hourly 
average  of  X^,  and  that  the  are  the  rapid  scan  measurements  of  wind  speed. 

Then 
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Oi  suppose  that  ST(Xj)  represents  the  hourly  standard  deviation  of  the 
w  i  ml  direction,  ciy  .  Then: 
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If  a  computed  quantity  ♦  (such  as  WS,  WD,  Og)  is  independent  of  the 
sampling  frequency  S  (at  least  to  the  extent  of  S  being  either  (6  sec)-^  or 
1  min'*),  then  one  would  expect  that  o(t)  will  have  a  low  value,  indicating 
that  the  $j's  are  very  close  in  value  to  ♦. 

Let  Vj:  represent  each  individual  wind-speed  measurement  and  let  0!  repre¬ 
sent  each  individual  wind-direction  measurement.  Then,  the  equations  actually 
used  to  compute  each  of  the  hourly  quantities  are  as  follows: 

average  wind  speeds: 
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vector  mean  wind  speed: 
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vect  '*•  mean  wind  direction: 
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The  standard  deviation  is  defined  by: 


i  I  lxi  -  «2 
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i /: 


,  and 


(H-b) 


( H—  7  ) 


In  computing  the  standard  deviations  of  the  above  quantities  one  should  note 
that  _for  wind  directions  the  smallest  angle  between  and  X  is  used;  that  is, 
X{  -  X  is  always  less  than  or  equal  to  it  radians. 


The  results  of  these  calculations  for  the  17  hr  of  rapid  scan  meteoro¬ 
logical  data  are  shown  in  Tables  H-l  through  H-6.*  The  important  thing  to 
notice  in  these  tables  is  the  low  value  of  o(®)  in  most  cases  —  implving  good 
agreement  between  the  $  and  values.  In  the  case  of  the  hourly  average 
wind  direction  Wl),  the  largest  value  of  o(WD)  one  encounters  is  2.4°.  For 
hourly  VMWD,  the  maximum  o(VMWD)  is  2.2°.  Both  values  are  well  within  the 
rangt  of  instrument  error.  For  WS  the  largest  percentage  deviation  (o($ g)/f  x 
1001.)  is  37,  and  for  VMWS  is  4%,  but  most  values  of  ait)  are  less  than  or  equal 
to  2  ■?.  of  tlie  value  of  #  itself.  In  the  case  of  eg ,  the  relative  values  of 
o(cg)  are  larger,  ranging  between  1  and  \2X  but  with  over  1/2  of  the  values  at 
or  h-low  bZ .  Therefore,  the  value  one  computes  for  a  quantity  during  a  given 
hour  will  not  depend  greatly  on  whether  the  data  used  for  the  computations  was 
taken  at  six-second  intervals  or  one-min  intervals  during  that  hour. 


In  addition,  from  the  values  of  the  standard  deviation  of  it  is 

apparent  that  the  values  computed  from  the  S=1  min-^  data  subsets  ( t>  * )  do  not 
ditfer  widely  from  their  mean.  This  lack  of  spread  implies  that  ♦*  is  repre¬ 
sentative  of  the  individual  t:  and  that  the  starting  time  of  the  first  minute 
i.s  not  important  --  any  set  of  minute-averaged  data  recorded  during  an  hour  i: 
just  as  valid  for  computat iona 1  purposes  as  any  other  set  that  starts  at  a 
d i  f  fe rent  t imo . 


So  it  is  apparent  that  any  m i nute -aver aged  data,  regardless  of  the 
starting  point  during  an  hour,  is  just  as  acceptable  for  use  in  computing 


★Tables  and  figures  appear  consecutively  at  the  end  of  this  appendix. 
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hourly  quantities  as  data  taken  with  a  higher  scan  rate,  at  least  down  to 
S  =  (6  sec)~l.  Hence,  if  computations  are  to  be  done  only  for  1-hr  sampling 
P<  riods,  nothing  substantial  is  gained  by  using  a  scan  interval  of  6  sec 
instead  of  l  min. 

In  the  same  manner,  hourly  average  pollutant  concentrations  can  be 
checked  for  dependence  on  the  scan  interval  S.  Along  with  the  meteorological 
data,  NO  and  NOy  were  recorded  at  6-see  intervals  during  the  17  hr  on  days 
145,  146  and  147. 

Tables  H-7  and  H-8  contain  the  results  of  the  same  calculations  using 
this  NO  and  NQy  data.  Again,  notice  the  good  agreement  indicated  by  the  low 
values  of  a  ( 4> )  for  most  hours.  For  N0X ,  the  value  of  o(N0x)  is  never  more 
than  2%  of  the  value  for  N0X  itself.  The  values  for  o(NO)  are  relatively 
larger,  ranging  from  1  to  12%  of  the  value  of  NO.  For  most  hours,  however, 
the  ratio  is  under  3%. 

Therefore,  as  in  the  case  for  the  hourly  meteorological  quantities,  data 
'  aken  with  a  scan  interval  of  one  minute  as  opposed  to  six  seconds  seems  to 
"..ake  very  little  difference  in  the  values  computed  for  the  hourly  average 

-  vic-.-ntrat  ion . 


3.  THE  EFFECT  OF  SCAN  RATE  ON  THE  AVERAGE  VERTICAL  WIND  SPEED  AND  THE 
STANDARD  DEVIATION 

Let  W  represent  the  vertical  velocity  of  the  air  as  measured  by  a  verti¬ 
cal  propeller.  The  data  used  to  compute  W  (the  average  vertical  wind  speed) 
and  ow  (the  standard  deviation  of  the  vertical  wind  speed)  comes  from  five 
hours  of  measurements  taken  at  various  sampling  frequencies.  The  notation 
used  here  is  the  same  as  used  in  the  previous  section,  with  the  exception  that 
now  there  is  more  than  one  rapid  sampling  rate.  However,  the  object  still  is 
to  ccxnpare  the  results  obtained  using  data  with  a  short-scan  interval  (the 
ompuced  quantity  denoted  by  ® )  to  the  results  obtained  using  S=1  min-1  data 
(the  quantities  denoted  by  $*). 

As  before,  the  S=1  min-*  data  is  obtained  by  selecting  points  from  tiu- 
rapid- scan  data  at  1-min  intervals.  Because  the  shortest  scan  interval  varies 
for  different  hours  there  will  be  a  different  number  of  subsets  of  minute- 
interval  data  in  each  case.  In  general,  for  this  minute-interval  data 

1=1 ,  2 . 60s,  where  S  is  the  sampling  frequency.  The  use  of  the  standard 

deviation  o(«)  to  compare  the  values  to  $  is  very  helpful  here  because 

■-? ( t )  accounts  for  the  different  numbers  of  values. 

Tables  H-9  and  H-10  show  the  values  calculated  for  and  W.  A  sampling 
period  of  1  hr  is  used  here,  and,  as  in  the  last  section,  the  relation  between 
Q  and  the  is  emphasized.  In  this  case,  however,  notice  that  the  agree¬ 
ment  between  the  $  values  and  the  values  is  very  poor,  as  indicated  bv 
the  large  values  of  a ( ♦ ) .  For  the  values  of  o(o^)  are  on  the  order  ot 

10%  of  the  value  itself.  For  W,  the  relative  values  of  o(W)  are  much 
higher,  ranging  between  13  and  73%  of  the  value  for  W.  Thus ,  it  is  obvious 
that  the  scan  rate  for  vertical  wind  speed  data  does  influence  the  i;ilc  elated 
,)t,rly  quantities.  This  influence  could  be  due  to  the  large  dev  i  at  ions 
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f  lie  'Min  t  i- r  til  in  the  vrl  ii  .il  w  i  ml  spued  --  not  i  c  e  tli.it  is  always  larger  than 
W.  A  mu— mi  nut  e-  sc an  interval  may  not  provide  enough  cl  < » tail  of  the  W  trace 
lor  usi>  in  calculating  an  hourly  quantity  that  is  representative  (i.e.,  one 
that  is  independent  of  the  starting  time  of  the  first  minute).  So  for  verti¬ 
cal  wind  speeds,  a  scan  rate  higher  than  once  per  minute  is  desirable. 

One  other  interesting  point  to  notice  is  the  diurnal  variations  of  W. 

Of  course,  no  definite  conclusions  can  be  drawn  from  five  hourly  averages,  but 
it  is  reassuring  to  see  that  the  highest  values  of  W  occur  in  the  afternoon, 
when  the  temperature  lapse  rate  is  maximum^  whereas  at  night  and  in  the  morn¬ 
ing,  when  conditions  are  more  stable,  the  W  values  are  lower. 


A.  THE  STANDARD  DEVIATION  OF  THE  WIND  DIRECTION  AS  A  FUNCTION  OF 
SAMPLING  TIME 


It  has 
c  q ,  depends 


been  suggested  that  the  standard  deviation  of  the  wind  direction, 
on  the  sampling  time  in  the  following  way: 


°0(T) 

w 


(H-8) 


whore  : 


a0< To) 
o0(  r) 


a  reference  sampling  time, 

standatd  deviation  of  the  wind  direction  calculated 
with  sampling  time  t0,  and 

standard  deviation  of  the  wind  direction  calculated 
with  sampling  time  x. 


The  17  hr  of  rapid-scan  meteorological  data  provide  a  data  base  with  which  to 
tesi  this  equation,  since,  even  for  short  sampling  times,  there  will  still  be 
a  suitable  number  of  data  points  available  for  the  calculations  of  o0  . 

Equation  (H-8)  will  be  tested  for  three  sampling  times:  x  =  3  min, 

10  min,  and  1  hr.  In  all  comparisons,  x  is  taken  as  the  shorter  sampling 
time.  However,  to  compare  0g ( x )  and  0g(xo),  the  average  of  the  0g(x)'s 
during  the  longer  sampling  period  x0  is  used  as  the  value  to  be  compared 
with  o0( x0  )  . 

Assuming  that  the  form  of  Eq .  (H-8)  is  correct  (i.e.,  that  it  is  a  power 
law  relation),  the  exponent  can  be  optimized  with  respect  to  our  data  in  order 
to  check  the  agreement.  Generalize  Eq .  (H-8)  by  using  a  variable  p  as  the 
exponent : 


VT) 


a 


e 


(T  ) 

o 


(H-9) 


Taking  the  Log  of  both  sides,  rewrite  Eq .  (H-9)  as: 
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log 


( H  —  1 0 ) 


\  e  o  / 


For  each  set  of  calculated  Og's  [oq(t)  and  Og ( tq ) J  a  value  for  p  can  be 
determined.  lire  mean  of  these  values  for  p,  denoted  as  p,  can  be  used  as  the 
value  for  the  exponent  in  Eq .  (H-9).  The  standard  deviation  is  then  an  indi¬ 
cation  of  how  representative  p  is  for  different  conditions. 

In  Tables  H— 1 1  to  H-13  a  value  for  p  is  calculated  for  each  of  the  17  hr 
of  rapid-scan  data  using  different  combinations  of  sampling  times.  Note  that 
the  mean  values  for  p  determined  from  data  are  higher  than  the  exponent  value 
used  in  Eq .  (H-8).  For  t=3  min  and  to=60  min,  {5=0.265  with  a  standard  devia¬ 
tion  of  0.116.  For  t=10  min  and  to=60  min,  {5=0.227  with  a  standard  deviation 
of  0.125.  For  t=3  min  and  to=10  min,  p=0.322  with  a  standard  deviation  of 
0.134.  This  last  comparison,  however,  is  between  the  averages  of  the  og 
(3  min)  values  and  Og(10  min)  values  during  each  hour. 

With  such  relatively  large  standard  deviations,  it  is  apparent  that  Eq . 
(H-9)  is  only  a  first  approximation  of  the  relationship  between  Og  and 
s  amp  ling  t ime . 


Instead  of  this  power  law  relation,  an  exact  identity  can  be  derived 
that  relates  a  og  calculated  for  a  given  sampling  period  to  the  og's  calcu¬ 
lated  for  a  set  of  shorter  sampling  periods  that  lie  within  the  longer  period. 
Specifically,  suppose  one  has  a  set  of  wind  direction  measurements  taken  with 
a  scan  interval  of  s  seconds.  If  the  longer  sampling  interval  is  denoted  by 
and  the  shorter  sampling  interval  by  Tg,  then  there  are  exactly  J=Tl/ts 
short  sampling  periods  contained  in  the  long  interval. 


If  j  is  an  index  that  represents  the  position  of  the  short  sampling 
period  during  the  longer  one  and  i  is  an  index  that  represents  the  position 
>>f  a  wind  direction  measurement  during  the  shorter  sampling  period,  then  the 
individual  wind  direction  measurements  can  be  represented  by  6jj,  where  i=l, 

2 . I  and  j  =  l ,  2 . J.  I  is  the  number  of  measurements  taken  during 

the  short  sampling  interval  Tg.  The  mean  wind  direction  during  the  longer 
sampling  period  is 


»UL> 


-III 


IJ 


I  J 


and  the  mean  wind  direction  during  each  short  sampling  period  is 


0  .(T 

J 


1  I 

I  i 


e . . . 

ij 


The  standard  deviation  of  the  wind  direction  lor  the  long  and  short  sampling 
periods  is  defined  as 


200 


w  = 


?  ?  (?u  )  -  e  . . )‘ 

1.1  i  j  L  lj 

,7  ?(?.(i  )-e..)‘l 

H  i  j  s  ij 


l  /  2 


1/2 


Bv  squaring  and  expanding  the  expression  for  0^(7),  adding  and  subtract¬ 
ing  the  second  term,  on  the  right  hand  side  of  Eq .  (H-ll),  and  substituting 
from  the  other  equations,  the  following  relationship  can  be  obtained: 


-(V 


=  1 1 


J  .1  6  j 


l 

+  — 
J 


I 
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(6(,l> 


V’s”' 


( H  —  1 1  ) 


This  is  an  exact  equation  that  relates  the  eg  '  s  calculated,  over  several  short 
sampling  periods  to  a  Og  that  would  be  calculated  if  these  short  periods  were 
combined  into  one  long  sampling  period. 


Using  one  hour  as  the  1  mg  sampling  period  i ^  and  i  and  10  minutes  as  two 
different  short  sampling  periods  tg,  values  tor  0j(tg)  and  Ogj(ig)  can  be  com¬ 
puted  from  the  rapid-scan,  wind  direction  data  using  the  same  methods  as  in 
Sec.  2.  Equation  (H-ll)  can  then  be  used  to  calculate  a  value  for  ce(l  hrj. 
Tables  H-14  and  H-l 5  compare  the  actual  hourly  Og ' s  and  the  Og's  predicted  by 
Eq.  (H-ll)  for  short  sampling  times  of  3  and  10  min. 

Since  Eq .  (h-ll)  is  exact,  the  close  agreement  between  the  calculated  and 
observed  oQ 1 s  are  not  unexpected.  With  one  exception,  the  two  values  lor  og 
are  within  1.5  degrees  of  each  other  for  the  case  tg=3  min  and  within  O.b 
degrees,  for  the  case  tg=10  min.  The  exception  in  both  cases  is  the  hour 
starting  at  0930  on  day  147  where  relatively  large  differences  between  the  Cg 
values  occur. 


All  these  differences  in  the  Og  values  are  probably  due  to  the  fact  that 
d  is  calculated  using  Eq .  (H-5)  of  Sec.  2  instead  of  the  linear  mean,  for  0 
used  in  tli  is  section  to  derive  Eq .  (H-ll).  The  larger  deviation  occurs  during 
an  hour  with  low  wind  speed  and  a  meandering  wind  direction. 


5.  MEAN  SQUARE  EDDY  VELOCITIES  AND  FILTERING  OF  THE  RAPID  SCAN  DATA 

In  order  to  define  what  one  moans  by  mean  square  a  long-wind  and  crosswind 
eddy  velocities  or  fluctuations,  one  must  first  define  t he  mean  velocities 
from  which  the  instantaneous  quantities  deviate.  In  the  case  of  the  a  long- 
wind  component,  the  instantaneous  velocity,  uj,  is  just  the  sum  ii  +  uj,  where 
li  is  the  mean  or  low  frequency  wind  component  and  17  is  the  eddy  or  high 
frequency  wind  component  . 

At  the  time  of  the  i^  sampling,  let  the  instantaneous  wind  speed  and 
direction  he  given  as  Vj_and  0£,  respectively.  (See  Fig.  H-l)*  If  the  "mean" 
quantities  are  given  as  V  and  Uj,  then  the  along-wind  fluctuation  is  just 

11 J  =  Vj  -  V  jCosO  [  ( H-  1  2a  ) 


♦Figures  and  tables  appear  consecutively  at  the  end  of  this  appendix. 
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while  the  cross-wind  fluctuation  is  just 


»  I 

v,  =  -vtsi n6  j 
win*  re : 


(11-1  2b) 


The  mean  square  eddy  velocities  are  then  given  by  the  expressions 


u 


,2 


N 


N 

I 


i=*  1 


(u!) 

i 


2 


(H-lJa) 


v’2  =  £  (vl)  2  ( H-l  3b ) 

i-1 

idi  ere  : 

N  =  the  number  of  events  available  during 
the  sampling  period  of  interest. 


To  obtain  the  "mean"  quantities  V£  and  0£,  a  number  of  averaging  (or 
tittering)  techniques  can  be  considered.  One  method,  of  course,  is  a 
straightforward  linear  average  of  the  Vj  and  6;  during  the  sampling  period. 

The  advantages  of  this  scheme  are: 

1)  it  sets  an  upper  bound  on  mean  square  eddv  velocities; 

2)  it  may  be  most  indicative  of  what  is  actually  happening 
to  the  pollutant;  and 

3)  it  is  a  standard  statistical  procedure  and  is  easy  to 
implement  and  interpret. 

Some  of  the  disadvantages,  however,  include: 

1)  The  mean  value  depends  strongly  on  the  duration  of  the 
sampling  period,  which  most  likely  is  of  arbitrarily 
selected  length. 

2)  Poor  connection  with  u  v  w  results,  especially  under 
light,  variable  winds. 

To  overcome  these  problems,  it  is  necessary  to  consider  a  moving  mean  for 
U j  and  vj.  In  this  wav,  large,  low  frequency  fluctuations  during  the  sampling 
period  do  not  introduce  large  deviations  in  the  aloug-wind  and  cross-wind 
velocities  (iij  and  v-)  because  those  low  frequency  components  become  part  ot 
the  moving  mean. 


One  way  of  obtaining  such  a  moving  mean  is  to  borrow  the  techniques  ot 
electronic  signal  filtering  and  actually  construct  a  low-pass  filter.  When 
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the  individual  9^  and  v^  are  run  through  the  filter,  the  output  is  a  time 
series  consisting  of  the  frequency  components  below  the  cut-off  frequency  of 
the  filter.  The  higher  frequency  components  will  have  been  removed,  and 
therefore  the  output  signal  is  suitable  as  a  "moving  mean"  in  that  when  it  is 
subtracted  from  the  complete  signal,  the  remaining  values  are  the  true  high 
frequency  turbulence  components  of  the  wind  speed  or  wind  direction. 

Consider  the  three  stage  low-pass  Butterworth  filter  of  Fig.  H-2,  where 
the  output  is  taken  across  resistor  R5 .  From  analysis  of  the  circuit, 
given  an  input  voltage  V£n(t),  the  output  current  1  through  resistor  R5  is 
governed  by  the  third  order  differential  equation: 


I  +  b  I  +  c  I  =  d  I 


,<t> 


( H— 1 4 ) 


where : 

the  dot  denotes  a  derivative  with  respect  to  time,  and  the  coefficients  have 
the  following  values: 

a  =  R1R5C2C4L3 

b  =  R3C2L3  +  R5C4L3  (H-15 ) 

c  =  K3R3C2  +  R1R5C4.  +  L3 
d  =  R3  +  R3 

To  make  this  filter  have  a  cut-off  angular  frequency  of  wc  (rad/s)  the 
values  of  the  components  should  be  set  to:* 

*1  -  1 


c2 


l3  = 


t;4  ~ 


"5  \> 

Actually,  u>c  is  the  half-power  point  on  the  frequency  response  curve  of  the 
filter,  which  means  that  at  oic ,  the  output  voltage  is  71%  of  its  value  in  the 
pass  band.  Because  of  the  construction  of  the  Butterworth  filter,  the  output 
voltages  even  in  the  pass  band  are  50%  of  the  input  voltage.  This  condition 
can  be  corrected  by  simply  doubling  the  output  of  the  filter. 

Using  the  above  values  for  the  filter  components  results  in  coefficients 
tor  the  differential  equation  of  the  form: 


*See  A  Handbook  on  Klectrical  Filters,  White  El ec t r omagne t ie s ,  Inc., 
Rockville,  Md . ,  1963. 
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=  _4_ 

C  '  “c 

d  =  2 

Thus,  the  differential  Eq .  (H-14)  becomes: 


2  •  •  •  4 

to  3  u  2 


I  +  21 


v.  (t) 
m 


l  H—  17) 


(H-18) 


VQut  then  obtained  from  the  relation 

Vout(t)  =  I(t)-R5-2  (H-19 ) 

since  the  output  voltage  is  developed  across  resistor  R5  and  must  be  doubled 
to  compensate  for  the  unavoidable  attenuation.  But  because  R5  =  1 , 

Vout(5)  =  2 1 ( t )  (H-20) 

where : 

I ( t )  *  the  solution  of  the  differential  Eq .  (H-18),  and 

vout  =  the  low  pass  filtered  signal  using  cut-off  frequency  uc . 

From  steady-state  analysis  of  the  filter,  the  steady-state  frequency 
response  is  determined  to  be 


(H-2 1  ) 


wh  ere  : 

A^n  =  amplitude  of  input  voltage  signal, 

Aquc  55  amplitude  of  output  voltage, 

S  «*  |  |zt  ]  |  =  (x2  +  y2)!/2, 

Zt  =*  Complex  impedance  *  x  +  iy , 
x  *  -bw2  +  d ,  and 

y  =  -aw^  +  cu. 

where : 

u  =  the  angular  frequency  of  the  input  signal,  and 
a,b,c,d  are  defined  in  Eqs.  ( H— 1 5 )  and  (H-17). 
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The  phase  angle  of  the  filter  circuit  is  given  by: 


b  =  tan-1  .  (>i-22) 

Tlie  steady- slate  response  of  tills  3-stage,  low-pa.-s  Hutterworth  filter  is 
shown  in  Fig.  H-3  for  a  cut-off  period  Tc  =  10  nun.  I  Tc  =  2»/<*>c.)  Notice 
that  the  filter  effectively  attenuates  the  frequency  components  greater  than 
ui(.  ,  even  though  the  rolloff  is  not  extremely  sharp. 

Now,  instead  of  using  strictly  voltage  levels  as  the  input  signal  to  the 
filter,  use  the  wind  speed  and  wind  direction  signals  as  if  they  were  volt¬ 
ages.  The  output  of  the  filter  will  then  be  the  lew-pass  filtered  wind  speed 
and  wind  direction  that  is  suitable  for  use  as  a  "moving  mean,"  since  it  only 
includes  the  low  frequency  components  of  the  input  signal.  The  advantages  of 
using  a  filter  tc  get  a  moving  mean  are: 

1)  the  method  ties  in  more  closely  with  classical 
turbulence  theory;  and 

2)  the  mean  follows  the  low  frequency  components 
of  the  signal  and  therefore  does  not  introduce 
arbitrarily  large  deviations  into  the  eddy 
veloc  it  ies  . 

However,  some  important  di sadvantages  are: 

1)  The  insight  into  total  pollutant  dispersal  power 
is  not  as  clear  as  with  a  fixed  mean;  and 

2)  The  filter  introduces  a  delay  time  into  the  output 
s  ignal . 

The  last  point  is  very  important  when  the  eddy  velocities  are  to  be 
computed,  since  a  matchup  must  be  made  between  the  total  signal  and  the 
filtered  signal,  if  their  difference  is  to  represent  only  true  deviations  of 
the  wind  speed  arid  wind  direction.  If  the  delay  time,  t  j  ,  introduced  bv  the 
filler  is  known,  then  the  input  signal  Vjn(t)  can  be  matched  with  the  output 
signal  Vout(t  +  t  )  to  get  the  desired  moving  average  of  Vjn(t)  at  time  t. 

Recall  that  Eq  .  (H-22)  defines  the  phase  angle  introduced  by  the  filter 
circuitry: 

e  =  tan-*  (-*) 

It  t(j  is  the  delay  time,  then  at  angular  frequency  w,  the  phase  angle  (which 
is  a  function  of  w)  is: 

6  ( co  )  =  wt^  .  Ui-23) 

Substituting  the  values  for  a,  b,  c,  and  d  from  Eq .  ( H — 1 7 )  into  the  values  for 
x  and  y  of  Eq .  (H-21),  one  obtains: 


and 


Expanding  this  and  neglecting  higher  order  terms  yields 


and  substituting  this  into  Eq  .  (H-23)  yields 


(H-24) 


tor  the  cases  of  u<<uic.  Also  note  that  td  can  be  rewritten  using  the  co¬ 
efficients  of  the  differential  equation: 


01-25) 


This  equation  gives  a  simple  approximation  for  tj,  but  it  is  strictly 
valid  only  for  the  frequencies  w<<uic.  For  frequencies  greater  than  u>c  ,  this 
is  not  a  problem  because  these  frequencies  will  have  been  filtered  out  of  the 
output  signal.  For  frequencies  around  wc  where  attomtat  ion  is  not  complete, 
t  j  is  underestimated  by  Eq .  (H-24). 

The  phase  angle  8  at  w  =  wc  is 

8  *  2.356  r ad i ans , 

from  Eq .  (11-2  3),  which  gives  the  exact  value  of  tj ,  one  obtains: 
d  u  wc 


206 


so 


f 


The  approximate  Eq .  (H-24)  would  predict 


which  amounts  to  an  18%  underestimate  in  tj  at  w(.  . 


For  example,  notice  in  Fig.  H-4  how  the  filtered  signal  of  the  wind 
speeds  (adjusted  tor  tj)  contains  only  the  low  frequency  components  of  the 
wind  speed  trace  and  serves  as  a  moving  average  for  the  u i ml  speed.  Also 
notice  that  the  approximation  for  tj  effectively  compensates  for  the  doin', 
time  introduced  by  the  filter  for  almost  every  frequency  component  of  the 
f i 1 1 e  red  signal  . 


Filtering  the  wind  direction  is  more  difficult  because  of  the  cut  at 
360°  .  So  instead  of  directly  filtering  the  0j  values,  sinOj  and  cosOj  are 
filtered  separately  and  then  recombined  to  obtain  0j  via  the  expression: 


( s in0  . ) 
_ if 

( cos 0  .  ) 
l  f 


(H-26) 


wh ere : 


f  denotes  that  the  filtered  values  of  sin0j 
and  cosO j  are  used. 

This  is  not  strictly  a  correct  procedure,  but  it  is  a  usetul  and  simple  tech¬ 
nique’  that  gives  reasonable  results.  In  Fig.  H-5,  notice  that  this  method 
returns  a  signal  for  0j  that  corresponds  well  to  what  would  be  expected  from 
a  moving  average  of  the  individual  0j.  Results  of  filtering  the  wind  speed 
and  wind  direction  data  from  a  different  hour  are  shown  in  Figs.  11-6  and  I  i  -  7 . 

Values  for  0j  and  Vi,  then,  can  be  obtained  by  filtering  the  individual 
0 |  and  Vj,  using  this  low-pass  Butterworth  filter.  The  actual  technique  used 
involves  numerically  solving  the  Differential  Eq  .  (11-18)  with  the  wind  speed 
signal  (Vj )  or  the  sine  and  cosine  of  the  wind  direction  signal  (sinG1  and 
cosOj)  as  the  input  forcing  function  Vjn(tK  The  filtered  output  VlM]t  is  then 
calculated  from  Eq  .  (H-2U).  The  value  for  Vj  is  obtained  Jirectlv  in  this 
way,  but  Oj  must  be  calculated  from  Eq .  (H-2b). 

To  solve  the  third-order  Differential  Eq .  (H-18),  it  must  he  rewritten  as 
three  first-order  differential  equations  that  are  then  solved  numerically  bv 
the  Runge-Kut l a  method.  This  technique  is  desirable  for  obtaining  the  filter 
output  because  the'  input  function  Vjn(l)  need  only  he  known  at  certain  t  ime 
intervals.  Hence,  a  sequence  ol  data  points  can  be  used  as  input  to  the 
numerical  method  without  knowledge-  of  the  continuous  limit  ion  these  points 
represent  .  Moreover,  this  method  produces  a  sequence  ol  points  a a  filtered 
output  that,  when  adjusted  for  delay  time,  can  he  used  as  a  rent  inuous  moving 
meat!  lor  display  purposes  (Figs.  II  4  through  11-7)  or  lor  ealculat  iutv.  . 

Equal  ions  (lt-12a  and  b)  c  an  now  be  used  to  ealculat  <•  the  uj  and  v,  . 

The  mean  square  eddy  veloeities  ( u 1  ^  and  v  1  ^  )  are  then  calculated  I  i  m;t  lip;, 

( 11-  l  la  a  nd  b  )  .  For  a  cut-off  period  of  10  min  and  a  sampling  t  ime  uf  mu 
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hour,  the  u 1  ^  and  v '  ^  computed  from  the  17  hr  of  rapid-scan  meteorological 
data  appear  in  Tables  H-16  and  H-17.  As  in  previous  sections,  the  computa- 
t  ions  are  done  for  both  s  =  (6  sec)-*  data  and  s  =  1  min-*  data. 

From  the  large  values  of  o($),  it  appears  that  the  tj  values  are  not 
necessarily  the  same  as  the  value  for  $ .  This  factor  indicates  that  averaging 
t  ime  has  a  noticeable  effect  on  the  calculation  of  u 1 ^  and  v 1 ^ ;  although  the 
filter  cut-off  period  Tc  might  also  have  an  effect,  since,  an  averaging  time 
of  one  minute  means  that  after  removal  of  the  moving  average,  only  fluctua¬ 
tions  with  periods  between  1  and  10  min  are  left.  This  band  is  not  very- 
wide,  and  it  may  not  contain  enough  information  to  compute  the  R.M.S.  odd\ 
v  e  1  oc  it  i  e  s  . 

The  selection  of  the  cut-off  period  for  the  low-pass  filter  Tc  =  10  min 
was  almost  arbitrary,  being  based  on  the  desire  to  center  Tc  between  the  sam¬ 
pling  time  (t  =  1  hr)  and  the  sampling  intervals  (6  sec  and  1  min).  In  this 
way,  the  filter  eliminates  fluctuations  with  periods  between  10  min  and  6  sec. 
These  periods  can  be  detected  with  an  averaging  time  of  6  sec  or  1  min  and 
will  contribute  to  the  values  for  the  mean  square  eddy  velocities.  Also,  Tt. 
is  less  than  1  hr  so  that  the  moving  mean  follows  the  fluctuations  with 
periods  between  1  hr  and  10  min  and  does  not  become  much  like  a  fixed  mean 
'which  would  happen  i f  Tc  >>  1  hr). 

For  comparison,  TabLes  H-18  and  H-19  show  the  mean  square  along  wind 
( u 1 2)  and  cross  wind  ( v 1  ^ )  eddy  velocities  computed,  using  a  fixed  linear  mean 
for  0j  and  vj.  The  hourly  average  wind  direction  (WD)  is  used  for  0j  and  the 
hourly  average  wind  speed  (WS)  is  used  for  v;. 

Notice  that  the  values  for  u ' 2  and  v'^  computed  with  a  fixed  mean  are 
higher  than  those  computed  using  a  moving  mean.  However,  the  values  of  o(4>) 
are  just  as  large,  indicating  poorer  agreement  between  the  ®  and  **  values 
in  the  case  of  u ' ^  and  v 1 ?  than  for  other  hourly  quantities. 

6.  APPROXIMATE  EQUATIONS  RELATING  VARIOUS  COMPUTED  PARAMETERS 

Manipulating  the  equations  used  to  define  the  hourly  quantities  WS,  WD, 
VMWS,  VMWD,  u 1  and  V 1  -  leads  to  some  new  equations  that  interrelate  some  of 
these  quantities.  Before  starting,  however,  it  is  necessary  to  derive  a 
method  for  determining  the  correlation  between  a  set  of  wind  speed  measure¬ 
ments  Vj  and  a  set  of  wind  direction  measurements  0j.  It  can  be  shown  that 
the  agreement  between  WD  and  VMWD  is  an  indication  of  how  uncorrelated  Vj 
and  0  j  are. 

From  Eq  .  ( H  —  7 )  in  Section  2: 

N 

£  Vj  sin«i 
-l  i«l 

VMWD  =  tan  — - - 

N 

Y  vi  cos8j 
i»  1 
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If  and  are  uncorrelated,  then 


"  N 

N 

- 

VMWD  -  tan”1 

I 

i-1 

Vi  I 

i-1 

s  ine  j 

N 

N 

I 

Vi  1 

cose  : 

Li-1 

i-1  J 

-  N 

- 

I 

s  in0  i 

-1 

=  t  an 

i*  1 

-  WD 

N 

I 

cose  ; 

_  i*  1 

ns  defined  by  Eq  .  (11-5)  of  Section  2. 

Thus,  if  VMWD  =  WD ,  then  V^  and  8j  are  uncorrelated.  This  result  is 
important,  because  the  separation  is  used  many  times  in  the  subsequent  deriva¬ 
tions.  From  Table  11-1,  notice  the  closeness  between  VMWD  and  WD;  if  one 
compares  the  differences  in  VMWD  and  WD  to  og  ,  the  largest  difference  as  a 
percentage  of  oe  is  found  to  be  13%,  but  most  values  fall  below  10%. 

A  relationship  between  o0 ,  WS,  and  VMWS  can  be  derived  as  follows: 

Given  that 


1/2 


,  and 


"s  -i  X  V 


VMWS 


l\  sln9i)2  Sv,  cos9i)' 


1/2 


Since  Vj  and  0j  are  uncorrelated 


VMWS  = 


-Tv' 

N  L  1 


(n  Z  sl,,9i^2  +  (l 


1/2 


/vMWsV 
\  ws  ) 


VMWS-)2  -  (Tine)2  >  (cosV)2. 

WS 


( H-2 / ) 
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If 


Chen : 


AS.  =0.  -  6  .  (Note  that  AS.  =  0.) 

l  l  l 


sin0£  =  sinS  cos(A0£)  +  cosS  sin(A0£>, 


and  so : 


s  in9 


—  V  sinScosCAS.)  +  cos0sin(A0.) 
N  A  l  l 


Since  sin0  is  an  odd  function,  A0 £  =  0  implies  that  X  sinA0£  =  0  also. 


Hence : 


sin0  =  sinfl  •  cos(A8) 
Similarly: 


cos0  =  cos0  •  cos(A8). 

Substitution  of  Eqs.  (H-28)  and  ( H— 2 9 )  into  Eq .  (H-27)  yields 


(  VMWsV 

\  WS  / 


(sin0)2  •  (cos( A0 ) )2  +  (cos0)2  •  (cos(A0))2 


=  (cos(A0))2. 

If  A0 £  is  small,  then  use  of  the  small  angle  approximation  gives 


(  VMWSj2 

V  MS  / 


(CO8(A0 ) )2 


'*?(» 


(A0.)2  (AB.)4 


iB  .  r  \ 


>  1  -  i  T  (A0.)2. 

N  L>  \ 

i 

This  equation  can  be  rewritten  by  using  the  definition  of  oq  : 


/  VMWsV 
V  MS  ) 


1  -  °02. 


wh  i.ch  is  valid  for  small  values  of  A0  £  and,  subsequently,  Og  . 


(H-28) 


(H-29) 


(H-30) 
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The  validity  of  this  equation  is  shown  by  Figs.  H-8,  H-9,  and  H-10  where  ctg 
calculated  from  Eq .  ( H— 3 1  )  is  plotted  against  the  oq  calculated  directly  from 

the  wind  direction  data.  In  Fig.  11-8,  the  hourly  quantities  used  are  those 
computed  from  the  six-second-scan  interval  data  (the  quantities  denoted  by  $) . 
In  Fig.  H-9,  the  means  of  the  quantities  computed  from  the  one-minute-scan 
interval  data  are  used  (the  quantities  denoted  by  ) .  In  Fig.  H-10,  the 
individual  quantities  computed  trom  the  one-minute-scan  interval  data  are  used 
(the  quantities  denoted  by  <fj). 


Notice  that  for  oq  values  below  40°,  Eq .  ( H — 3 1 )  predicts  a  value  for  ag 
that  corresponds  very  well  to  what  is  actually  observed.  For  o0  values 
smaller  than  20°,  Eq .  ( H— 3 1 )  predicts  a  value  that  is  within  10%  of  the 
observed  og  value.  For  ag  values  of  less  than  40°,  the  predictions  are  good 
to  within  20%.  In  almost  all  cases,  however,  Eq .  (H-31)  slightly  underpre¬ 
dicts  the  c g  value.  Figure  H-10  shows  that  even  the  individual  t>j  values 
give  good  results,  and  more  important,  that  the  scatter  of  points  is  along 
the  one-to-one  correspondence  line.  This  shows  that  Eq .  (H-31)  is  applicable 
to  any  set  of  hourly  averaged  quantities  computed  from  mi  nut e- int erva 1  data. 

Manipulating  Eqs.  (H-4)  and  ( H—  7 )  of  Sec.  2  and  Eqs.  (H-13a)  and  (H-13b) 
of  See.  5  in  other  ways  leads  to  some  different  relationships  between  these 
quantities.  Two  of  the  more  useful  ones  are  listed  here  without  derivation: 

I 

l/2,  and  ( H-32 ) 


V 


V  2 


u’^  =  oy2  (  1  -  oq2)  ( H-33) 

where  : 

V  =  average  wind  speed  (WS),  and 

ov  =  standard  deviation  of  the  wind  speed. 

These  equations  use  the  small  angle  approximation  to  simplify  expressions  with 
Oq  and  assume  that  u ' ^  and  v 1 ^  are  computed  using  a  fixed,  linear  mean,  as  was 
done  in  Tables  H-15  and  H-16. 

In  Fig.  H  —  1 1 ,  Og  calculated  from  Eq .  (H-32)  is  plotted  against  the 
calculated  directly  from  the  wind  direction  data.  In  Fig.  H  —  1 2 ,  Eq.  (11-33)  is 
tested  in  the  same  manner.  For  both  cases,  the  ♦  values  are  used  in  the  cal¬ 
culations.  Referring  to  Fig.  H— 11,  notice  that  Eq  .  (H-32)  is  accurate  to 
within  10%  up  to  Jq  values  of  20°  and  is  accurate  to  within  30%  up  to  40°. 

From  Fig.  11-12,  however,  it  is  apparent  that  Eq  ■  (11-33)  is  valid  over  only  a 
ve i v  narrow  range.  The  calculated  value  for  u ' ^  is  within  10%  of  the  actual 
value  up  to  0.2  in’’  /  s  •  and  is  accurate  to  within  40%  up  to  2.0  m^/  with  one 
except  ion  --  a  point  calculated  during  an  hour  with  very  low  wind  speed  falls 
well  out  of  the  40%  hand. 
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Fig.  H-l.  Definition  of  Coordinates 


•1 


Hg'  H~3>  ^equeney  Response  of  3-Stage  Low-Pass  Lutterworth  Filter 
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ATTENUATION  (dB) 


TIME  l 

Fig.  H-4.  Filtered  Wind  Speed  f 
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Table  H-l.  Average  Wind  Direction  ($  =  WD) 


D.y 

Hour 

WD 

(deg) 

o(VD) 

(deg) 

1! 

145 

1200 

213.5 

1.6 

213.5 

1.6 

1300 

208.9 

1.2 

208.9 

1.2 

1400 

258.3 

1.7 

258.3 

1.7 

1500 

268.6 

1.6 

268.6 

1.6 

1600 

276.8 

1.1 

276.8 

1.1 

1700 

297.2 

0.9 

297.2 

0.9 

1800 

331.9 

1.1 

331.9 

1.1 

1900 

348.6 

0.9 

348.6 

0.9 

2000 

336.7 

0.3 

336.7 

0.3 

2100 

337.6 

0.2 

337.6 

0.2 

2200 

77.4 

1.7 

77.4 

1.7 

146 

1650 

250.7 

1.0 

250.  7 

1.0 

1750 

268.7 

0.4 

268.7 

0.4 

1850 

302.2 

0.3 

302.2 

0.3 

147 

0830 

141.7 

1.1 

141.7 

1.1 

0930 

168.9 

2.4 

168.9 

2.4 

1030 

223.8 

2.0 

223.8 

2.0 

Tab le  H-2. 

Vector  Mean  Wind 

Direction  ($ 

=  VMWD) 

Day  Hour 

VMWD 

o (VMWD) 

VMWD* 

o  (VMWD*’ 

(deg) 

(deg) 

(deg) 

(deg) 

145  1200 

212.4 

1.5 

212.4 

1.5 

1300 

214.0 

1.4 

214.0 

1.4 

1400 

262.4 

1.7 

262.4 

1.7 

1500 

267.5 

1.5 

267.5 

1.5 

1600 

275.3 

1.1 

275.3 

1.1 

1700 

293.8 

0.8 

293.8 

0.8 

1800 

330.6 

1.2 

330.6 

1.2 

1900 

349.1 

1.1 

349.1 

1.1 

2000 

336.3 

0.3 

336.3 

0.3 

2100 

337.0 

0.2 

337.0 

0.2 

2200 

71.1 

2.2 

71.1 

2.2 

146  1650 

250.6 

0.5 

250.6 

0.9 

1750 

266.8 

0.6 

266.8 

0.6 

1850 

301.0 

0.4 

301.0 

0.4 

147  0810 

140.9 

1.1 

140.9 

1 . 1 

0910 

167.3 

1.9 

167.3 

1.9 

1030 

229.5 

1.6 

229.  5 

1.6 

Table  H-3.  Standard  Deviation  of  the  Wind 
Direction  ($  =  on) 


Day 

Hour 

°9 

(dag) 

o(o#) 

(dag) 

Oft* 

(dag) 

(dal) 

US 

1200 

54.3 

3.6 

54.1 

3.6 

1300 

50.3 

1.8 

50.3 

1.8 

1400 

51.1 

1.6 

51.0 

1.6 

1500 

25.3 

2.0 

25.2 

2.0 

1600 

29.5 

2.6 

29.4 

2.6 

1700 

30.9 

1.2 

30.9 

1.2 

1800 

18.3 

0.8 

18.3 

0.8 

1900 

11.4 

1.4 

11.2 

1.4 

2000 

5.1 

0.4 

5.0 

0.4 

2100 

6.4 

0.2 

6.4 

0.2 

2200 

54.0 

0.6 

54.0 

0.6 

146 

16S0 

19.5 

1.2 

19.4 

1.2 

1750 

23.3 

0.7 

23.3 

0.  7 

1850 

15.0 

0.4 

14.9 

0.4 

147 

0830 

22.2 

1.4 

22.1 

1.4 

0930 

69.7 

2.3 

69.6 

2.3 

1030 

42.1 

1.1 

42.0 

1.1 

Table  H-A. 

Average  Wind  Speed  ($  ■  WS) 

Day 

Hour 

vs 

(■/a) 

O(VS) 

(»/») 

WS* 

(■/a) 

*7  (VS*) 

(■/a) 

145 

1200 

2.67 

0.09 

2.67 

0.09 

1300 

2.55 

0.03 

2.55 

0.03 

1400 

3.52 

0.06 

3.52 

0.06 

1500 

3.99 

0.06 

3.99 

0.06 

1600 

3.57 

0.04 

3.57 

0.04 

1700 

3.07 

0.03 

3.07 

0.03 

1800 

3.54 

0.04 

3.54 

0.04 

1900 

2.61 

0.03 

2.61 

0.03 

2000 

2.59 

0.01 

2.59 

0.01 

2100 

2.45 

0.01 

2.45 

0.01 

2200 

1.41 

0.01 

1.41 

0.01 

146 

1650 

4.46 

0.0* 

4.46 

0. 09 

1  7  VI 

4.07 

o  m 

4.0  7 

out 

1850 

3.37 

0.07 

3.  37 

0.03 

147 

0830 

1.97 

0.05 

1.97 

0.05 

09  30 

1.82 

0.03 

1.82 

0.03 

1070 

2.77 

0.04 

7.  77 

0  04 

Table  H-'). 

Vector  Mean 

Wind  Speed 

O  - 

VMWS) 

Day  Hour 

VMWS 

(»/■) 

c(VMWS) 

(«/•) 

VMWS* 

(m/s) 

a (VMWS*) 
(m/s) 

145  1200 

2.09 

0.07 

2.09 

0.07 

1300 

2.03 

0.04 

2.03 

0.04 

1400 

2.95 

0.05 

2.95 

0.05 

1500 

3.69 

0.07 

3.70 

0.07 

1600 

3.23 

0.05 

3.23 

0.05 

1700 

2.67 

0.05 

2.67 

0.05 

1800 

3.36 

0.05 

3.36 

0.05 

1900 

2.56 

0.03 

2 . 5t 

0.03 

2000 

2.58 

0.01 

2.58 

0.01 

2100 

2.43 

0.01 

2.43 

0.01 

2200 

0.80 

0.01 

0.80 

0.01 

146  1650 

4.22 

0.10 

4.22 

0.10 

1750 

3.76 

0.04 

3.76 

0.04 

1850 

3.27 

0.03 

3.27 

0.03 

147  0830 

1.85 

0.04 

1.85 

0.04 

0930 

1.14 

0.04 

1.14 

0.04 

1030 

2.17 

0.06 

2.17 

0.06 

♦ 

Table  H-6.  Standard  Deviation  of  the  fable  H-7.  Hourly  Average  NO 


Wind 

f  peed 

($  = 

0  ) 

V 

(* 

=  NO) 

Day 

Hour 

"v 

U(v!  ) 

V 

T* 

v 

o(o  *) 
v 

Day 

Hour 

NO 

c(N0) 

NO* 

'  (NO*  3 

(m/s) 

(to/ s) 

(m/ s) 

(ra/s) 

Cppb) 

(ppb) 

(ppb) 

<ppb) 

145 

i?no 

1 .4*4 

0.08 

1.44 

0.08 

— 

...  _.  .. 

- . . 

-•  - 

145 

1200 

3  .60 

0.04 

3.66 

o.  V« 

1  300 

1.40 

0.06 

1.40 

0.06 

1  300 

3.53 

0.08 

3.  53 

0 . 08 

1400 

1.79 

0.11 

1.79 

0.11 

1400 

3.68 

0.05 

3.68 

0.0r- 

I  500 

1.20 

0.05 

1.20 

0.05 

1500 

3.62 

0.08 

3.62 

0.08 

1600 

1.44 

0.06 

1.43 

0.06 

1600 

4.13 

0.07 

4.13 

0.07 

1700 

1 .  31 

0.07 

1.30 

0.07 

1700 

3.40 

0.  06 

3.40 

0.06 

1800 

0.81 

0.02 

0.80 

0.02 

1800 

3.5? 

0.07 

3.5? 

0.07 

1900 

0.47 

0.03 

0.40 

0.03 

1900 

3.11 

(0  05 

3.1  1 

0.05 

.’000 

0.  37 

0.02 

0.37 

0.02 

2000 

3.35 

0.03 

3.3  5 

0.  03 

:ioo 

0.  3? 

0.01 

0.33 

0.01 

2100 

3.5/ 

0.03 

3.  57 

0.0  3 

2?oo 

0.31 

0.01 

0.31 

0.01 

2  200 

3 . 68 

0.04 

3.68 

0 . 04 

146 

1*50 

1 .  16 

0.08 

1  .  36 

0.08 

146 

16  50 

1  .48 

0.0  5 

1  .  4R 

0.05 

17  50 

f) .  84 

0.04 

0.84 

0.04 

17  50 

0.88 

O.  06 

0.88 

0.  06 

1850 

0.  5/ 

0.03 

0.51 

0.03 

18  50 

0.  78 

« » .  09 

0.  78 

.  f»‘i 

r*7 

OR10 

0.  7/ 

O.  04 

0.71 

0.04 

147 

08  30 

0.75 

0 . 09 

0.7  5 

•-)  v« 

0930 

0.97 

0.04 

0.  97 

0.04 

09  30 

1.13 

n .  08 

1.13 

r. 

1030 

1.20 

0.04 

1.20 

0.04 

_ 

_ _ 

_ 

.  _  _ _ _ 

1030 

1 .07 

0.07 

1.07 

r>  ■ 

J 


22  r> 


Table  H-8.  Hourl v  Average  NO  (t>  -  (NO  ) 

x  x 


Day 

Hour 

NO 

X 

(ppb) 

a  (NO  ) 

X 

(ppb) 

NO  * 

X 

(ppb) 

o(NO  *) 

X 

(ppb) 

145 

1200 

8.56 

0.11 

8.56 

0.11 

noo 

8,81 

0.07 

8.81 

0.07 

1400 

8.98 

0.05 

8.98 

0.05 

1500 

9.46 

0.10 

9.46 

0.10 

1600 

10.00 

0.07 

10.00 

0.07 

1700 

8.54 

0.06 

8.54 

0.06 

1800 

7.45 

0.04 

7.45 

0.04 

1900 

8.42 

0.06 

8.42 

0.06 

2000 

17.79 

0.07 

17.79 

0.07 

2100 

28.84 

0.03 

28.84 

0.03 

2200 

26.79 

0.07 

26.79 

0.07 

146 

1650 

5.71 

0.13 

5.71 

0.13 

1750 

4.21 

0.07 

4.21 

0.07 

1850 

6.98 

0.10 

6.98 

0.10 

147 

0830 

6.24 

0.07 

6.24 

0.07 

0930 

6.48 

0.06 

6.48 

0.06 

1030 

4.77 

0.11 

4.77 

0.11 

Table 

H-9 . 

Standard  Deviation  of 
Wind  Speed  ($  =  ow) 

the  Vertical 

Day 

Hour 

S 

(sec) 

°w 

(m/s) 

o(ow) 

(m/s) 

('w* 

(m/s) 

0  (Ow*) 

(m/s) 

144 

1200 

2.4 

0.535 

0.055 

0.528 

0.055 

1300 

2.4 

0.429 

0.042 

0.423 

0.041 

1430 

3.3 

0.372 

0.025 

0. 369 

0.025 

145 

2305 

3.0 

0.096 

0.010 

0.095 

0.010 

146 

1130 

3.0 

0.722 

0.052 

0.717 

0.052 

Table 

H-10 . 

Average 

Vertical  Wind 

Speed  (<t 

=  W) 

Day 

Hour 

S 

(sec) 

W 

(m/s) 

o(w) 

(m/s) 

W* 

(m/s) 

a  (k*) 
(m/s ) 

144 

1200 

2.4 

0.166 

0.064 

0.166 

0.064 

1300 

2.4 

0.184 

0.058 

0.184 

0.058 

1430 

3.3 

0.209 

0.044 

0.209 

0.044 

145 

2305 

3.0 

0.075 

0.010 

0.075 

0.010 

146 

1130 

3.0 

0.096 

0.071 

0.096 

0.071 

Table  11-11.  Calculation  of  p  for  Sampling 
Times  of  3  min  and  60  min 


Day 

Hour 

(3  win) 

«>tl  (60  min) 

P 

145 

1200 

26. 903 

54.266 

0.  .'  14 

1300 

23.690 

50.309 

0.2  51 

1400 

24.422 

51.092 

0.746 

1500 

1  5 . 59  2 

25.309 

0.  162 

1600 

1 7  .09  7 

29.487 

0. 182 

1700 

18.013 

30 . 89 1 

.) .  ISO 

1 800 

9  .061 

1 8 . 104 

0.2  35 

1900 

5.49  3 

11.359 

0.243 

2000 

2.817 

.5.054 

0. 19  3 

2100 

2.121 

6.407 

0 . 369 

2200 

8.343 

54.046 

0.624 

146 

1650 

l 1 .063 

19.454 

0.188 

1750 

7.029 

? 3.261 

0.  359 

1850 

4.950 

14,954 

0.  369 

147 

0830 

16.623 

22.184 

0.096 

09  30 

27.929 

69 . b6S 

0.305 

1030 

18.854 

42.079 

0.268 

Mean  Value  for  p 

-  0.265 

Standard  Deviation  •  0.116 

Table  H-12. 

.  Calculation 

of  p  for 

Table  H-13. 

Cal cul at 

ion  of  p 

f  U  l* 

Sampling  Times  of 

Sampling 

T i me s  of 

10  min  and 

60  min 

i  min  and  10  min 

Day  Hour 

(10  min)  og 

S 

a 

►— 

3^ 

P 

Day  Hour  0 

r»  (  3  min) 

tj  (10  win) 

P 

14  5  1200 

40.723 

54.266 

0. 160 

145  1200 

26.903 

40.773 

j.  3-4 

I  300 

35.202 

50.  309 

0.199 

non 

2  3.  (<40 

35.202 

2.323 

1401) 

39 . 580 

51.092 

0.142 

1400 

24.42? 

39  .  580 

•'  .vi 

1500 

20.367 

23.309 

0.121 

1500 

1 5.592 

7 0 . 367 

0 .22  2 

1600 

20.3S6 

29  .487 

0.201 

1600 

17.097 

20. 586 

0.154 

1700 

26.337 

30.891 

0.089 

l  700 

18.033 

76.  137 

0.  315 

1800 

12.383 

18.304 

0.209 

1800 

9  .061 

12.585 

0.273 

1900 

6.82  5 

11.359 

0.  284 

19(10 

5.49  3 

6 .82  5 

0.180 

2000 

1.  378 

5.054 

0.225 

2000 

2.837 

3.  378 

0,145 

2100 

3.6V8 

6.407 

0.319 

2100 

2.121 

3.618 

0.444 

2  700 

19.3/1 

54 .046 

0.567 

2200 

H  .  34  7 

19.571 

0.708 

146  16 SO 

16.905 

19.454 

0.078 

146  1650 

1 1 .067 

16,905 

0.342 

1  7  SO 

17.  IV 

23.261 

0.353 

17  50 

7  .929 

12 . 357 

0. 369 

1  HSO 

7.388 

14  .954 

0.394 

1850 

4 . 9  50 

7.3*8 

»un 

IV  OHIO 

19 . 768 

22.184 

0.064 

l  4  7  OH  )<) 

1 6.623 

) 9 , 76R 

■).  144 

09  10 

44.445 

69.665 

0.251 

09  40 

77.«»?9 

44.445 

n .  >k», 

10  30 

29  .405 

42.079 

0.200 

10  30 

18.854 

79 .40  5 

f..  369 

Mean  Value  for 

p  -  0.727 

Mean  Value  for  p 

-  n . 372 

Standard  Deviation  -  0.175 

Standard  Deviation  *  n. 1*4 

227 


ililc  il-14.  Comparison  of  the  Actual 
<>  to  the  Value  Calculated 
from  Kq.  (H— 11)  using  Tg=3  min 


OgWeg-) 


Day 

Hour 

Ac  tual 

Calculated 

145 

1200 

54.3 

55.3 

1300 

50.3 

51.8 

1400 

51.1 

52.6 

1500 

25.3 

25.4 

1600 

29.5 

29.5 

1700 

30.9 

31.1 

1800 

18.3 

18.5 

1900 

11  .4 

11.5 

2000 

5.1 

5.1 

2100 

6.4 

6.4 

2200 

54 .0 

54.2 

146 

1650 

19.5 

19.5 

17  50 

23.3 

23.3 

18  50 

15.0 

15.0 

147 

0830 

22.2 

22.2 

0930 

69.7 

77.9 

1030 

42.1 

42.2 

Table  H-15.  Compur  iser.  at  the  Ae  t 

a  ..  to  the  Value  Cal  eulat  «■ 
0  Q 

from  Kq.  (H— 11'  usine.  t  =10  min 


o^tdeg.) 


Day 

Hour 

Actual 

Calculal ed 

145 

1200 

54.3 

54 .9 

1300 

50.3 

50. 1 

1400 

51.1 

49.7 

1500 

25.3 

25.3 

1600 

29.5 

29.  •- 

1700 

10.4 

31.7 

1800 

18.3 

18. 3 

1900 

11.4 

11.6 

2000 

5.1 

5.1 

2100 

6.4 

6.4 

2200 

54 .0 

54 .2 

146 

1650 

19.5 

19.5 

1/50 

23.3 

2  3.3 

1850 

15.0 

14.9 

147 

0830 

22.2 

22.1 

0930 

69.7 

75.7 

1030 

42.1 

42.4 

Table  H-16.  Mean  Square  Along  Wind  Eddy  Velocity  (f  -  £**) 


Hour 

(■2/«2) 

oCu72) 

(■2/»2) 

<-2/.2) 

oCu12*) 

(*2/»2) 

14S 

1200 

0.959 

0.125 

0.898 

0.109 

1300 

0.778 

0.097 

0.741 

0.090 

1400 

1.440 

0.286 

1.367 

0.276 

1500 

0.934 

0.135 

0.849 

0.105 

1600 

0.729 

0.102 

0.667 

0.081 

1700 

0.576 

0.068 

0.531 

0.051 

1800 

0.401 

0.046 

0.377 

0.040 

1900 

0.095 

0.031 

0.078 

0.025 

2000 

0.026 

0.007 

0.021 

0.004 

2100 

0.021 

0.004 

0.019 

0.004 

2200 

0.017 

0.003 

0.015 

0.003 

U6 

1650 

0.776 

0.134 

0.697 

0.108 

1750 

0.356 

0.066 

0.305 

0.043 

1850 

0.093 

0.015 

0.082 

0.011 

U7 

0830 

0.  342 

0.065 

0.302 

0.051 

0930 

0.454 

0.063 

0.412 

0.046 

1030 

0.700 

0.087 

0.640 

0.063 
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'able 

H- 1 7 .  Mean 

Square  Cross 

Wind 

Eddy  Velocity  ($  ■ 

Day 

Hour 

0(^1 

\jT  2  i. 

oCT7"2*) 

(»2/»2> 

(.n2/n2) 

(■ 2/s2) 

145 

1200 

0.764 

0.263 

0.823 

0.256 

1300 

0.703 

0.168 

0.  711 

0.168 

14  00 

1.757 

0.208 

1.830 

0.195 

1500 

1.026 

0.092 

1.012 

0.091 

1600 

0.899 

0.102 

0.853 

0.090 

1700 

0.963 

0.180 

0.991 

0.178 

1800 

0.308 

0.025 

0.298 

0.023 

1900 

0.088 

0.061 

0.081 

0.060 

2000 

0.016 

0.005 

0.013 

0.004 

2100 

0.007 

0.002 

0.006 

0.001 

2200 

0.032 

0.007 

0.035 

0.007 

146 

1650 

0.840 

0.124 

0.840 

0.124 

1750 

0.310 

0.042 

0.286 

0.035 

1850 

0.092 

0.018 

0.080 

0.013 

147 

0830 

0.256 

0.036 

0.242 

0.033 

0930 

0.345 

0.061 

0.398 

0.030 

1030 

0.604 

0.071 

0.602 

0.071 

Table  H-18.  Mean  Square  Along  Wind  Eddy  Velocity  (4  =  u'~) 
Computed  Using  a  Fixed  Mean 


Day 

Hour 

u '  * 

((n?/s;’) 

n(u'7) 

(in' /s'  ) 

(oF/h2) 

Ku77* 

( n/  /  F.  •' ) 

1A5 

1200 

1.122 

0.  312 

3.119 

0.312 

1300 

2.668 

0.202 

2.671 

0.20  2 

1400 

4 .6  U 

0.38/ 

4.626 

0.  3H7 

1500 

l.79» 

0.142 

1.791 

0.142 

1600 

2.262 

0.207 

2.256 

0.207 

1700 

1.2118 

0.136 

1.816 

0.136 

1800 

0.691 

0.036 

0.689 

0.0  Jh 

1900 

0.159 

0.017 

0.158 

0.017 

2000 

0.117 

0.012 

0.137 

0.012 

2100 

0.109 

0.007 

0.109 

0.007 

2200 

0.626 

0.022 

0.675 

0.022 

146 

1650 

1.619 

0.249 

1.907 

0.248 

1750 

0.910 

0.058 

0.919 

0.058 

1850 

0.  134 

0.031 

0.3V1 

0.031 

147 

08  10 

0 . 584 

0.062 

0.582 

0.062 

0930 

2 . 000 

0.149 

1.997 

o.uy 

1010 

2.223 

0. 196 

2.217 

0.196 
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H 

Table 

H-19 . 

Mean  Square  Cross  Wind 

Eddy  Velocity 

($  =  v77) 

Computed  Using 

a  Fixed 

Mean 

Day 

Hour 

v77 

oG1?) 

v,z* 

o  (v77*) 

(m2/s2) 

(m2/s2) 

(m2/ s2) 

(m'/s2) 

i 

145 

1200 

2.066 

0.203 

2.065 

0.203 

i 

i 

i 

i 

1300 

2.003 

0.223 

1.999 

0.223 

1400 

2.686 

0.239 

2.683 

0.239 

1500 

1.197 

0.159 

1.965 

0.159 

1600 

2.233 

0.128 

2.235 

0.128 

1 

1700 

2.393 

0.326 

2.392 

0.326 

1800 

1.184 

0.065 

1.178 

0.065 

1900 

0.275 

0.057 

0.273 

0.057 

2000 

0.045 

0.008 

0.045 

0.008 

2100 

0.073 

0.005 

0.073 

0.005 

- 

1 

2200 

1.225 

0.048 

1.224 

0.048 

146 

1650 

2.211 

0.139 

2.210 

0.139 

1750 

2.315 

0.149 

2.31-'. 

0.149 

f 

1850 

0.633 

0.041 

0.633 

C.041 

147 

0830 

0.383 

0.039 

0.384 

0.039 

0930 

1.428 

0.080 

1.430 

0.080 

- 

1030 

2.622 

0.180 

2.625 

0.180 

_ 

2  30 

A 

Al'I'hNI)  i  X  I 


l'KKI.IMI  NARY  TIMh  .T.KIKS  ANAJ.YMt, 


1.  INTRODUCTION 

It  is  often  possible  to  describe  complex  physical  phenomena  and  mech¬ 
anisms  better  in  the  frequency  domain  than  in  the  time  domain.  The  time- 
dependent  nature  of  air  pollution  and  meteorological  data  and  the  likelihood 
of  serial  correlation  in  it  makes  obvious  the  need  for  analyzing  such  data  in 
terms  of  its  frequency  composition. 

Many  statistical  analysis  methods  that  operate  on  the  data  in  the 
time  domain  make  no  use  of  the  time  dependency  therein.  For  example,  in  the 
estimation  of  the  frequency  or  cumulative  frequency  distributions,  the  data 
are  sorted  into  order,  after  which  the  time  dependency  information  is  no 
longer  available.  An  important  property  of  time  series  analysis,  on  the  other 
hand,  is  that  much  of  the  time  dependency  in  the  data  is  preserved. 

The  use  of  time  series  in  the-  analysis  of  air  pollution  monitoring 
has  been  suggested  or  carried  out  by  several  authors  including:  Harlow 

and  S i ngpurwa 1  la , ^ ^  Saltzman,^-^  Marcus,^-^  and  Tiao  and  Hammi ng , to 
mention  only  a  few. 

Data  reduction,  harmonic  analysis,  and  filtering  of  Williams  air  quality 
data  was  performed  previously  on  raw  data  from  the  HP9825.I-®  Spectral 
analysis  was  instrumental  in  identifying  certain  instrumentation  problems 
that  were  exhibited  by  a  sawtooth  pattern  in  the  raw  data.  These  problems, 
which  were  diagnosed  early  in  the  experimental  program,  were  subsequently 
c  orrected  . 

An  extensive  treatment  of  spectral  analysis  of  time  series  has  been 
given  by  Koopmans ; ^“5  time  series  data  analysis  and  theory  are  discussed 
extensively  by  Br il 1 inger ; and  a  bibliography  of  time  series  analysis 
papers  has  been  given  by  Wold.I-^ 

The  principal  goal  of  spectrum  analysis  is  to  decompose  the  power  of 
the  given  time  series  into  its  harmonic  components,  that  is,  to  estimate  the 
power  spectrum  from  the  available  data.  The  estimated  spectrum  can  then  be 
used  to  gain  information  about  the  mechanism  that  generated  the  data. 
The  stat  istical  theory  of  spectral  anlaysis  has  been  based  on  the  hypo¬ 
theses  that  the  underlying  process  is  stationary  and  Gaussian,  that  the 
process  mean  is  zero,  and  that  the  spectrum  is  continuous.  Some  of  the 
preprocessing  operations  performed  on  the  data  before  spectrum  analysis  is 
done  (for  example,  removal  of  the  DC  component  of  the  series)  are  intended  to 
bring  the  data  into  reasonable  conformity  with  these  hypotheses. 

Two  models  of  spectral  analysis  have  been  used  here,  univariate  and 
bivariate.  Univariate  spectral  analysis  involves  the  computation  of  the  power 
spectrum  or,  as  it  is  also  called,  the  spectral  density.  The  power  spectrum 
is  the  Fourier  transform  of  the  autocorrelation  function.  The  Fourier  trans¬ 
form  of  a  sequence  of  autocorrelations  contains  precisely  the  same  infor¬ 
mation  found  in  the  original  autocorrelations  their,  re  1  ves .  The  intensity  of 
the  spectrum  at  any  given  frequency  provides  an  estimate  of  the  power  or 


energy  at  that  frequency.  The  accumulated  spectrum  is  essentially  the  first 
integral  of  the  spectral  density  and  can  be  used  to  indicate  the  range  of 
frequencies  that  contain  a  substantial  proportion  of  the  total  power  of  the 
time  series. 

Similarities  between  the  spectra  of  two  series,  such  as  peaks  at  similar 
frequencies,  may  raise  the  possibility  that  the  series  are  related.  To 
investigate  such  possibilities,  the  cross  spectrum  of  two  series  m.<y  be 
computed.  This  computation  is  an  extension  of  the  definition  of  the  spectrum 
and  is  usually  estimated  by  smoothing  the  cross  per iodogram ;  this  is  performed 
in  the  bivariate  spectral  analysis  of  two  series.  Bivariate  spectral  analysis 
allows  for  the  decomposition  of  two  series  into  various  spectral  parameters, 
equivalent  to  descriptive  statistics.  These  include:  the  estimated  spectral 
density  of  each  series;  the  phase  function  that  indicates  the  relative  shift 
of  harmonic  components  of  the  two  time  series  at  the  same  frequency;  the 
transfer  function  that  indicates  the  linear  relationship  between  the  two 
series  in  the  same  way  that  regression  coefficients  relate  the  linear  rela¬ 
tionship  between  two  variables;  and  the  coherence  function  that  provides  an 
estimate  of  the  frequency-dependent  correlation  between  two  series.  It  is 
interesting  to  note  that  if  the  phase  angle  is  a  rapidly  varying  function  of 
frequency,  the  estimated  coherence  can  be  biased  downward  to  an  extent  that  a 
strong  coherence  will  be  masked. 

In  summary,  spectral  analysis  provides  a  means  of  analyzing  serial 
correlations  in  the  data,  as  well  as  a  method  of  examining  frequency  dependent 
correlations.  It  is  indispensable  when  searching  for  periodicity  and  gross 
variations  of  power  with  frequency  in  the  data. 

Spectral  analyses  were  performed  on  a  continuous  period  of  one-hour 
averages,  covering  approximately  one  month,  starting  August  1,  1976.  Missing 
data  and  outliers  identified  by  an  interface  program  to  the  AQAM  data  tape 
were  assigned  values  via  a  simple  linear  interpolation  of  adjacent  values. 
Univariate  and  bivariate  spectral  analyses  were  performed  using  methods 
outlined  by  Koopmans.I-^  Before  spectrum  analysis  was  applied  to  the  data, 
the  mean  (DC  component)  was  removed  from  the  series  to  bring  the  data  into 
conformity  with  stat ionar ity . 

Fi.gure  1-1*  shows  the  time  series  for  observed  CO  at  station  2.  While 
this  series  appears  to  exhibit  a  pronounced  24-hour  period,  one  finds  upon 
close  examination  that  the  daily  period  has  slight  relative  shifts.  Where  the 
periodicity  is  in  fact  a  combination  of  several  other  periodic  phenomena, 
closer  investigation  is  required.  It  is  felt  that  the  roughly  24-hour  periods 
are  actually  beats,  and  unless  the  series  is  filtered,  these  will  not  exhibit 
themselves  in  the  spectrum  as  at  first  anticipated. 

Figure  1-2  shows  the  spectrum  for  this  time  series;  the  peaks  at  16-17 
hours,  8  hours,  and  6  hours,  are  of  interest.  The  largest  peak  at  approxi¬ 
mately  16  hours  is  an  artifact  of  the  sunrise-sunset  period.  (Note  that 
this  data  was  drawn  from  summer  months.)  Periodicity  reflected  by  peaks  at 
higher  frequencies  appears  to  be  an  artifact  of  source  emissions  activity 
during  the  day.  It  is  felt  that  the  higher  frequencies  at  8  and  6  hours  may 

•Figures  appear  consecutively  at  the  end  of  the  appendix. 
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I)i-  related  to  aircraft  activity.  'flu?  x-ax  i  h  of  t  f  i  c  spec  train  has  been  n-plot- 
1  ed  in  terms  of  period  but  wan  original  Ly  computed  as  a  frequency  in  terms 
of  cycles  in  units  of  1/(2*M*DT),  where  DT=1  hoar,  and  M  is  tin-  truncation 
parameter  that  has  a  value  of  LOO.  Hie  truncation  parameter,  also  referred 
to  as  the  number  of  lags,  controls  the  smoothing  o f  the  per  iodog ram  and 
has  been  chosen  approximately  equal  to  10":  of  the  data  length,. 

Figure  1-3  shows  the  accumulated  spectrum  for  the  CO  series.  It  is 
interesting  to  note  that  the  higher  frequencies,  above  the  frequencies  related 
to  diurnal  fluctuation,  appear  to  contribute  significantly  to  the  total  power 
of  the  observed  CO  series. 

Figure  1-4  shows  the  time  series  for  observed  N0X  at  station  2  during 
this  same  period.  That  the  structure  of  this  series  is  very  similar  to  the  CO 
series  is  confirmed  by  the  locations  of  the  peaks  in  the  spectral  density  shown 
in  Fig.  1-5.  Figure  1-6  shows  the  coherence  between  the  CO  and  N0X  series. 
It  is  interesting  to  note  that  while  the  peak  coherence  is  associated  with  the 
diurnal  cycle,  the  coherence  remains  strong  at  the  higher  frequencies,  sug¬ 
gesting  that  these  two  pollutants  may  be  associated  with  the  same  set  of 
sources . 

Figure  1-7  shows  the  time  series  of  aircraft  CO  emissions  on  ground 
level  lines.  This  emission  rate  is  computed  in  AQAM  on  the  basis  of  the 
hour-by-hour  records  of  aircraft  activity.  The  spectrum  of  this  series  shows 
the  strong  diurnally  induced  peak  at  16  hours  along  with  small  enhancements  at 
8  and  6  hours  where  significant  peaks  were  observed  in  the  observed  CO 
concentration  spectrum  (Fig.  1-2).  The  coherence  between  the  spectra  of 
observed  CO  concentrations  and  aircraft  CO  emission  rate  also  shows  peaks 
corresponding  to  these  periodicities;  however,  the  wildly  oscillating  magni¬ 
tude  of  the  coherence  and  accompanying  phase  (not  shown)  are  suggestive  of 
computational  instabilities  possibly  being  driven  by  the  strong  low-frequenv 
peaking  in  both  spectra. 


2.  CONCLUSIONS  AND  RECOMMENDATIONS 

While  the  background  concentrations  due  to  strong  natural  cycles  an 
significant,  harmonics  at  higher  frequencies  thought  to  be  associated  with 
emissions  source  activity  have  been  identified.  The  presence  of  large  peaks 
in  the  spectrum,  have  the  effect  of  badly  biasing  estimates  in  the  low  end 
of  the  spectrum  due  to  the  inevitable  side  lobe  distortion  of  the  spectral 
estimates  (filter  leakage).  A  quasi  difference  filter  can  be  used  to  balance 
the  spectrum  of  a  time  series  with  a  large  peak  at  Low  frequencies.  The 
balancing  process  is  called  prewhitening  by  blackmail  and  Tukey  (  1  959). 

A  prewhitening  filter  can  be  used  to  help  smooth  the  larger  natural  peaks, 
improving  the  resolution  of  the  higher  frequency  spectra. 

Spectral  analysis  can  be  used  to  study  the  effect  of  filters  used  on 
the  air  quality  time  series  and  also  in  the  design  of  such  filters.  Kilters 
are  used  for  purposeful  modification  of  a  discrete  time  series.  It  is 
important  to  understand  the  operation  of  the  filter,  and  to  know  how  to  select 
the  "parameters"  of  the  filter  to  achieve  specific  objectives  of  data 
modi f  icat ion . 


3  1 ! 


The  spectral  analyses  presented  here  illustrate  the  benefits  of  analyzing 
time-dependent  air  quality  data  in  the  frequency  domain.  Estimation  of 
correlation  vis  a  vis  the  coherence  function  between  two  series  provides 
important  insights  into  complex  physical  systems  generating  the  data,  and 
the  ability  of  the  AQAM  to  reproduce  features  prevalent  in  the  observations. 
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